Foro - Perl en Español

por **alexclipse** » 2013-12-26 07:13 @342

Hola, estoy intentando crear un programa que me cuente el número de '1' y '0' dentro de un archivo. Los '1' representan el número de metilaciones de una cadena de ADN. Hay dos enzimas de restricción, H y M, y por el momento quiero contar el número de metilaciones en cadenas cortadas por la enzima H. Para ello he hecho este programa:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

#!/usr/bin/perl
 
$metilaciones = 'C:\BIF\practica5\practicaperl\metilacion.txt';    # ruta del archivo
open( METILACIONES, $metilaciones );                               # abrir archivo y asignar filehandle
@lineas             = <METILACIONES>;
$numerolineas       = @lineas;                                     # contar número líneas del archivo
$contador           = 0;                                           # contador del número de H
$contador2          = 0;                                           # contador de líneas
$cuentametilaciones = 0;                                           # contador de 1
while ( $contador2 < $numerolineas ) {                             # mientras haya líneas
    if ( $lineas[$contador2] =~ 'H' ) {                            # si la línea tiene H
        $contador = $contador + 1;                                 # sumo uno al contador de H
        if ( $lineas[$contador2] =~ '1' ) {                        # si la línea con H tiene 1
            $cuentametilaciones++;                                 # sumo uno al contador de metilaciones
        }
    }
    $contador2++;                                                  # siguiente línea
}
print "$contador\n";                                               # imprimir número de H
print "$cuentametilaciones\n";                                     # imprimir número de metilaciones
close METILACIONES;                                                # cerrar archivo
Coloreado en 0.003 segundos,  usando GeSHi 1.0.8.4

El problema es que lo que estoy haciendo es sumar 1 a las metilaciones cada vez que la línea tiene '1', no contando el número de '1'. ¿Cómo hago para contar los '1'? Muchas gracias.

#!/usr/bin/perl
use v5.10;
use autodie;                    # «Es mejor morir que regresar con deshonor» --proverbio Klingon
 
my $haches = 0;
my $unos   = 0;
 
open my $METILACIONES, '<', 'kk.txt';
 
while (my $linea = <$METILACIONES>) {
    next if -1 == index $linea, 'H';    # pasamos a la siguiente si la $linea no tiene 'H'
 
    $haches++;
 
    $unos += $linea =~ tr/1/1/;         # contamos los '1' que hay, y los acumulamos
}
 
say "H: $haches. 1: $unos.";
 
close $METILACIONES;
Coloreado en 0.001 segundos,  usando GeSHi 1.0.8.4

por **alexclipse** » 2013-12-27 04:35 @232

Muchas gracias por la ayuda, aunque conseguí solucionarlo de otro modo.

Foro - Perl en Español

Contar caracteres de un archivo

Contar caracteres de un archivo

Publicidad

Re: Contar caracteres de un archivo

Re: Contar caracteres de un archivo

¿Quién está conectado?