Foro - Perl en Español

por **beafm** » 2011-01-09 10:02 @460

Hola de nuevo:

Estoy atascadísima en el mismo sitio. A ver, tengo lo siguiente en mi programa:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

#PRIMER MARCO DE LECTURA
 
     my (@posinicio);
     my (@posfinal);
     my ($posicion) = 0;
 
     my (@genes) = seleccionar_secuencias( $posicion, $DNA, \@posinicio, \@posfinal, $longitudmin );
 
print "@genes\n";
 
 
#SEGUNDO MARCO DE LECTURA. 
 
     my (@posinicio1);
     my (@posfinal1);
     $posicion = 1;
   
     my (@genes1) = seleccionar_secuencias( $posicion, $DNA, \@posinicio1, \@posfinal1, $longitudmin );
 
print "@genes1 \n";
Coloreado en 0.003 segundos,  usando GeSHi 1.0.8.4

La subrutina es:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

sub seleccionar_secuencias {
 
  my($posicion, $secuencia, $posinicio, $posfinal, $longitudmin) = @_;
  
  my($codoninicio);
  my($codon);
  my($i);
  my(@genes);
 
for ($posicion; $posicion < length ($secuencia) - 2 ; $posicion = $posicion + 3) {
 
          $codoninicio = substr ( $secuencia, $posicion, 3 );
 
                        my($pos2) = $posicion;
                        my($pos3) = $posicion;
                        my($gen);
 
              if ($codoninicio =~ /ATG/i) {
                        
                        do{
                            $codon = substr ( $secuencia, $pos2, 3 );
                   
                            $gen .= $codon;
 
                            $final = $pos2 + 1;
 
                            $pos2 = $pos2 + 3;
 
                          }until ( $codon =~ /TA[AG]/i or $codon =~ /TGA/i);
 
                          if ( $longitudmin < length ($gen) ) { 
 
                                    $genes[$i] = $gen;
        
                                    $$posinicio[$i] = $pos3 + 1;
        
                                    $$posfinal[$i] = $final;
 
                                    $i = $i + 1;
                
                          }
                                                  
                }
          }
             
  return @genes;
 
}
Coloreado en 0.001 segundos,  usando GeSHi 1.0.8.4

Una parte de $DNA sería lo siguiente:

Sintáxis: [ Descargar ] [ Ocultar ]

Using text Syntax Highlighting

GATCTCGCAACTCTTAATAAAGCTAATTCCTGCTAATTCGCTATACACTAAATCTTAATAAACTTGCTTC

ATTATTTTCTCGATAAATATTTTATTTATAACAGTTTGTTAAATTGACAGGGTCCTTGCTAAGTTTCTAT

ATACAACATAAACATACTATACGAAATAACATTTTTCTCAAGGTATGTAGAGTATAAATGCTCCAACGTG

GCACAATTTCATTTCGTTTAGTTTACTAAAAACCCATATTCGCGTCGCACTCACACAGAATCTAGAGTTT

AGTACATTTTACCGGCTCAACAATTAACCATTCAAATTGTATTTCAAGTATGAGTAACAATCGAATAAAA

TATTGTATTCTTGCTTTCCTATAAAAGTATTTTGATAATCGTAAAGCTGGGAATACCATTTGCCTGTACT

TTGGCAGTTATGGTGCTGCCATGTTGCTCCTTTTATAGTTATATCCTGTATGTTAAAAACTTTTATGGAA

CTATCGTACTTCTGTTTAACAGCATCATAAGCCTTAAATGCTAAATTTCGATTATTAACTGGTATTCCTG

AACAATCACGTATACCGTTTTTGAAGTCTTGGACCATGAAATCGTCGCCTGAGCGTCTATCAAACGCACC

TAATTGCATCAGCTGAGTCATGTTCATATTATGTCTTGAACTGTCATCGAATTCAATACTCATAGTAGTC

AACGCAAAGAGAGATACTGAACATATTGGCACATCATCTAAATTGCCATTTCATAATGAAAAGCAACCCT

TCTGACATTCGTTTATTTAAAGAAACACTCGTCAATCAAACGTGTGTAATCACCATCATTAACTTCGACT

GTATTAAATAATACAAGTCTACCCAAACTCAGTTATCGTTCCGTTAGTCACATTACTATATGCATTAAAA

CAATCAACTACTGCCTATCTGAACTCAGACGAATTGGTATTATAAGCATCAATAGTTAGGACACACATTA

TGAGTAATAAAACTAGAGTAATCATTGCTATATCGAAGAGTATAACTTATTTGGGTTGAAAACTATCACT

AGACGTAAAATCTAATATACTAAACGAGTTACGAGGTTCAGTTGTAGTTCTTTATTTATCATCTTTTCTA

ATGATTGCTATCATCGGGTATTCTATGTAAGTAGACAAATTTGTAATACATTTACGTTATTCATTTCCTT

TCCTTGAAGGTGTTGTAATATTCAATAGATCCCTTGAAAATTACATCAGCTCTATAAAGGTTGGCTCATT

GCTAAAACTCAAGTGTTGCTCAAATGGTTTAATAATGTCAGCTACGCAGTTTTGGTATCTGTTTAACGAT

ACGTATAAGTACGATGTGTTTTCCTTCTGACCTAGTATATTACAGTACACAGGTAATATAAGTAGCACCT

TCTATAGATATTTATTAAACTGTTACCTATCACTGTAATAGCTTATTACAGAACTGAGGAACGAGGTTAG

TTCTTCTCAATTTATGGCATTTAGAATTATTTTTATTATTACGTTTGATTAGTTTTGTTCACAACTAGCA

ATATGTCTAAATAACTTATCGTGTAAAGTTGATAGAGCAAATAGAACACGTAGCTCAGAAAACACAATTC

CTCTAATTTAATCCTCCCAAGTCCACGACATACCCACTACAGGTTGGTCTCTCGACTGATGGATAGAGAG

TCTATCAATCTTGTATTCTTAGCGGTGAACCACCACCATCCAAAAAACTAGTGAACTACTTATATCTTAG

TAATGAGTCGAGTAACAAACTTTTTACGATAATGATACTATGAAGAAGCATCTACCAGTTGTGCATATGT

CAAAATAGCGTTATAATTTTATAACATCACACAAATATATTAATCGAATAGTCCTAAGTCACACCGTAGT

TGACAAAGCGGCGTTGCTAGTTCCATAAGTAATATTTGTACGAGCATATACGATTAATTGATTCGTAAAA

TATGATCTTTTTATGCTAATACTATTATTCATTACTGTCTACCATATTTTGTAGTTCATCATTATTTTAT

TTATACAATTACTAGCTTTGATTTTTGAATTGTTATTTGGTTATTGATAGCAATTGATATTAAACGCTTA

TTCTTTAATATTATTGCACATACATCGTATTTACTGAATATGCCGGATAAAAACAACAGACTTGTGTCCCColoreado en 0.000 segundos,  usando GeSHi 1.0.8.4

Cuando hago print @genes me sale todo bien, pero luego para print @genes1 no me sale nada. Me sale esto:

Sintáxis: [ Descargar ] [ Ocultar ]

Using bash Syntax Highlighting

tuxdistro@linux-swp1:~/programa_de_informatica>perl migetorf.pl ficheroFASTA ficheroSALIDA 300Coloreado en 0.002 segundos,  usando GeSHi 1.0.8.4

y no me sale nada más... Creo que el resto del programa lo tengo bien hecho, pero si esa subrutina no me funciona, no me va nada y no puedo avanzar en nada más. Y es que si pongo $posicion = 2; que también lo tengo que hacer, sí me da, pero con el 1 no. He probado con el depurador, y para $posicion = 1; no me imprime nada en $codon, y en $gen siempre tiene el mismo gen, pero en @genes no hay nada de nada.

Agradecería mucho si alguien me puede decir qué hago mal, porque he probado a cambiarlo de mil maneras, y nada.

Hola de nuevo:

Ahora mismo no tengo dudas sobre el programa que estoy haciendo, pero si le doy a que me imprima algo por pantalla, me dice "out of memory". Tengo todos los pasos por referencia que he podido, y también he intentado recoger lo máximo posible en subrutinas. Tengo solo una subrutina más complicada, que es la de extracción de los ORF.

¿Si le doy a que imprima en un fichero también se queda sin memoria? ¿Hay alguna forma de hacer que funcione?

Otra cosa. He visto que he recibido un mensaje privado, pero no puedo verlo. ¿Por qué?

Gracias

por **explorer** » 2011-01-12 07:53 @370

¿Cuánta memoria tiene el sistema en donde estás? Si estás en Linux, ejecuta el comando free.

Es muy raro de que un programa de este tipo agote la memoria... podría tratarse de un bucle infinito, o, efectivamente, que estés duplicando la información.

Si estás llamando a las subrutinas con un argumento que sea la secuencia de ADN, y esa secuencia tiene cien mil bases, pues toda esa cantidad de memoria es la que se reserva en cada paso (aunque luego la debería de soltar, desde luego).

Para leer los mensajes privados tienes que entrar en la sección "Bandeja de Entrada".

por **beafm** » 2011-01-13 01:57 @122

Ok, ya me funciona el programa, creo que tenía una función que no se me ejecutaba bien.

Al final he hecho lo que me dijo usted de hacer un array de array y luego aplicarles el sort() para ordenarlos. ¿Existe alguna manera de ordenar los elementos de mayor a menor? Es que el sort() lo hace de menor a mayor, pero una de las veces que lo uso tiene que ser en el orden contrario.

Gracias.

por **explorer** » 2011-01-13 05:57 @290

Si, por ejemplo, tenemos

sort { $a <=> $b } @array

para ordenar numéricamente el @array, de menor a mayor, con

sort { $b <=> $a } @array

lo ordena de mayor a menor.

Más información en perldoc -f sort.

por **beafm** » 2011-01-13 15:32 @688

Ok. Ya lo último que pregunto, porque ya lo tengo hecho.

Tengo que imprimir una cabecera y la proteína que sea dividida en tramos de 60 aminoácidos. Yo lo hago, pero me mete un espacio que no debe salir. No sé por qué me mete ese espacio.

El código es éste:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

imprimir_secuencia ($posicionesordenadas1[$uu]->[2], 60);
 
sub imprimir_secuencia {
 
    my($cadena, $longitud) = @_;
 
    for ( my($pos); $pos < length($cadena) ; $pos = $pos + $longitud ) {
        print SALIDA substr($cadena, $pos, $longitud), "\n";
    }
}
Coloreado en 0.001 segundos,  usando GeSHi 1.0.8.4

Lo que me sale a mí es:

Sintáxis: [ Descargar ] [ Ocultar ]

Using text Syntax Highlighting

>1-300000_1  [5451 - 5954] Schizosaccharomyces pombe 972h- chromosome II, complete sequence

 MSNPESLKKQVEPPGYNELFMVEDVCNVDLEQGLDLCKPEKVNKQSQRSRQSRQSLFTNT

IKPQKDKMNIKTNKIKEFLNDLFTEFSKFHNSYYPDGRISTRSNFRWPLLIIWSIIIVFA

VDKKFEVQKFLSIWINENRFYSEIWVPIAIYVCLLVLMLLSLICMFPLColoreado en 0.000 segundos,  usando GeSHi 1.0.8.4

Y lo que me tiene que dar es:

Sintáxis: [ Descargar ] [ Ocultar ]

Using text Syntax Highlighting

>1-300000_1 [5451 - 5954] Schizosaccharomyces pombe 972h- chromosome II, complete sequence

MSNPESLKKQVEPPGYNELFMVEDVCNVDLEQGLDLCKPEKVNKQSQRSRQSRQSLFTNT

IKPQKDKMNIKTNKIKEFLNDLFTEFSKFHNSYYPDGRISTRSNFRWPLLIIWSIIIVFA

VDKKFEVQKFLSIWINENRFYSEIWVPIAIYVCLLVLMLLSLICMFPLColoreado en 0.000 segundos,  usando GeSHi 1.0.8.4

Es idéntico, menos en el primer espacio de la proteína.

por **explorer** » 2011-01-13 16:07 @713

Yo no veo problemas en el código, así que sospecho que el espacio en blanco está en la secuencia leída.

Haz la prueba de poner un

print "[", substr($posicionesordenadas1[$uu]->[2], 0, 1), "]\n";

Y así sabrás si hay o no, realmente, un espacio en blanco en la primera posición.

por **beafm** » 2011-01-15 09:59 @458

He hecho lo que me ha dicho, y me sale [M] lo que quiere decir que no hay espacios... además, al comenzar el programa, cuando saqué la secuencia ya me aseguré de quitarle todo lo que tuviera extraño, incluidos los espacios...

Así que el espacio a lo mejor es cosa de la pantalla, no lo sé...

Previamente a hacer ese print() tengo este otro, no sé si el problema podría estar ahí:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

print SALIDA ">$cabecera[0]_$numero1  [$posicionesordenadas1[$uu]->[0] - $posicionesordenadas1[$uu]->[1]] (REVERSE SENSE) $cabecera[1]\n ";Coloreado en 0.001 segundos,  usando GeSHi 1.0.8.4

Yo diría que no hay problema, porque \n me salta de línea, no me mete espacios...

No sé qué puede ser, la verdad.

Foro - Perl en Español

Más problemas con los ORF

Atasco... no avanzo al variar la posición

Publicidad

Re: Más problemas con los ORF

Re: Más problemas con los ORF

Re: Más problemas con los ORF

Re: Más problemas con los ORF

Re: Más problemas con los ORF

Re: Más problemas con los ORF

Re: Más problemas con los ORF

Re: Más problemas con los ORF

Re: Más problemas con los ORF

¿Quién está conectado?