Foro - Perl en Español

¡Hola perleros!
Tengo dos tablas de datos que quiero cruzar. En la más pequeña (llamémosla "pequeña") tengo posiciones cromosómicas en esta forma:



Así para cada cromosoma, cientos de entradas.

Y quiero cruzar esta tabla con otra principal (llamémosla "general") enorme con casi todas las posiciones en el genoma que tiene tres columnas (cromosoma, coordenada y un valor numérico), ordenado por cromosoma y posición.



Lo que quiero hacer es buscar las posiciones de "pequeña" en "general" y sacar a otro archivo la posición y el valor numérico para dicha posición. Es decir quiero filtrar el archivo "general" en función de las posiciones en "pequeña".

¿Cómo puedo manejar estos datos? Me imagino que el archivo "general" no lo puedo cargar en memoria por lo que éste lo tendría que abrir línea a linea. Pero entonces no sé cómo hacer la búsqueda de cada entrada de "pequeña". Como los archivos están ordenados por cromosoma, si quiero buscar el valor numérico para, por ejemplo, chr18 12224555 tendría que leer prácticamente todo "general" hasta llegar a la sección de chr18.

Espero haberlo explicado de manera entendible. Gracias por vuestra ayuda.

D.

Yo creo que sí puedes leer "general" a memoria y convertirlo en un hash, salvo que me digas que tiene varios gigabytes de tamaño o que debes realizar esta operación varias veces al día, y tarda mucho tiempo en cada ocasión, para leerlo. En esos casos, hay que pensar algo distinto. Mientras tanto, sigue una solución normal.

El hash sería, poniendo de valor de la clave, el valor de la primera columna y el valor de la segunda columna de general. Y como valor del hash Algo así:

$cromosoma{"$primera,$segunda"} = $tercera;

Luego, solo te queda abrir "pequeña", línea a línea, y sacar el valor que buscas:

say $cromosoma{"$primera,$segunda"};

Gracias explorer. Pues el "general" ocupa 700Mb, pero tampoco se va a ejecutar muchas veces el programa, la verdad.

Voy a por ello.

Gracias,

D.

hola de nuevo explorer. He intentado hacer algo como lo que me has sugerido, pero creo que tengo algún error de concepto y lío en el uso de las variables, porque al ejecutar el programa me da errores de declaración de variables:



aquí debajo puedes ver mi código. Pienso que al declarar la variable $cromosoma en el ámbito general, podría usarla dentro del bucle, pero parece que no la he declarado bien?

Perl dice que no sabe de dónde viene $cromosoma. Es decir, que no ha sido declarada antes de ser usada por primera vez.

Debes declararla con un my() antes de usarla.

Y, sobre todo, darle una valor cada vez.

Por otra parte, me parece que la línea 34 está mal... hay un operador de asignación, dentro de un if().

Gracias explorer, he modificado el if() de la línea 34 para que la comparación sea de strings, pero sigo sin entender por qué tengo que declarar $cromosoma, si la he declarado globálmente en la línea 13, ¿o no?. Quiero usar la variable $cromosoma que viene del archivo "general" (ARGV[0]) y compararla con la variable $chr declarada dentro del while (<IN2>). Si vuelvo a declararla con un my() dentro del while, ¿no se comportaría como una variable local?

En la línea 13 estás declarando (y definiendo) la variable $cromosoma como local dentro del primer while() (fíjate en las llaves, que son las que delimitan el contexto en que es conocida la variable).

En cambio, en la línea 34, el programa se encuentra con una variable desconocida, $cromosoma (porque era local en el while anterior, no global), y además, no tiene ningún valor (no se lo has puesto).

Prueba a poner esta línea:



que, aunque no es lo que te he respondido en mi respuesta, es así como estás construyéndolo en la primera estructura.

Esta es mi versión: