Foro - Perl en Español

por **angelluigi** » 2016-11-16 08:55 @413

Buenas.

Me están pidiendo que diga las posiciones en un array que presenten el valor mínimo. Ando un poco liado con esto. Si alguien me puede ayudar...

Muchas gracias.

#!/usr/bin/env perl
use v5.20;
 
my @array = map { 10 + int rand 90 } 1 .. 100;
 
#say "@array";
 
my @posiciones_minimo;
my $valor_minimo = 10E10;
 
while(my($i, $valor) = each @array) {
    if ($valor_minimo > $valor) {
        $valor_minimo = $valor;
        @posiciones_minimo = ();
    }
 
    if ($valor_minimo == $valor) {
        push @posiciones_minimo, $i;
    }
}
 
say "Mínimo encontrado: $valor_minimo";
say "Posiciones en las que se encuentra: @posiciones_minimo";
 
__END__
Coloreado en 0.002 segundos,  usando GeSHi 1.0.8.4

Una posible ejecución:

Sintáxis: [ Descargar ] [ Ocultar ]

Using text Syntax Highlighting

$./kk.pl

Mínimo encontrado: 10

Posiciones en las que se encuentra: 43 59 90Coloreado en 0.000 segundos,  usando GeSHi 1.0.8.4

Bueno, al final, no es tan complejo... con apenas 5 líneas principales lo resolvemos.

La explicación del funcionamiento te lo dejo como deberes para casa.

Naturalmente, se puede hacer de otras formas.

por **angelluigi** » 2016-11-16 19:22 @849

Muchas gracias, acabo de empezar el master en bioinformática y todavía ando un poco verde.

Lo he probado y la solución concuerda perfectamente. Más o menos creo que lo he entendido. Lo que pasa es que no he visto el bucle while() ... each(), que has usado. Yo he sacado el valor mínimo con un foreach y pensaba sacar la posición, o sea, el índice, con un bucle for.

Te dejo el código para ver si puedes sacarlo con un bucle for, gracias y perdón por las molestias.

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

my @skew = ();
 
$contador_G = 0;
$contador_C = 0;
$min_skew   = 0;
 
$secuencia = "CCTATCGGTGGATTAGCATGTCCCTGTACGTTTCGCCGCGAACTAGTTCACACGGCTTGATGGCAAATGGTTTTTCCGGCGACCGTAATCGTCCACCGAG";
 
for ( my $i = 0; $i < ( length $secuencia ); $i++ ) {
    $base = substr( $secuencia, $i, 1 );
    if    ( $base =~ /C/ ) { $contador_C = $contador_C + 1; }
    elsif ( $base =~ /G/ ) { $contador_G = $contador_G + 1; }
    my $diferencia_GC = $contador_G - $contador_C;
 
    push @skew, $diferencia_GC;
    foreach (@skew) {
        if ( $_ < $min_skew ) { $min_skew = $_; }
 
    }
 
}
 
$skew = join( '', (@skew) );
print $skew, "\n\n\n";
 
$valor_minimo = $min_skew;
 
my @posiciones_minimo;
 
while ( my ( $i, $valor ) = each @skew ) {
    if ( $valor_minimo > $valor ) {
        $valor_minimo      = $valor;
        @posiciones_minimo = ();
    }
 
    if ( $valor_minimo == $valor ) {
        push @posiciones_minimo, $i;
    }
}
 
print @posiciones_minimo;
 
Coloreado en 0.002 segundos,  usando GeSHi 1.0.8.4

por **explorer** » 2016-11-17 18:28 @811

No necesitas hacer dos bucles, ya que con uno te vale, ya que estás recorriendo toda la secuencia:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

#!/usr/bin/env perl
 
my @skew;
 
my $contador_G;
my $contador_C;
my $min_skew   = length $secuencia;     # hay que inicializarlo a un valor alto o máximo
my @posiciones_minimo;
 
my $secuencia = "CCTATCGGTGGATTAGCATGTCCCTGTACGTTTCGCCGCGAACTAGTTCACACGGCTTGATGGCAAATGGTTTTTCCGGCGACCGTAATCGTCCACCGAG";
 
for ( my $i = 0; $i < length $secuencia; $i++ ) {
 
    $base = substr $secuencia, $i, 1;
 
    if    ( $base eq 'C' ) { $contador_C++; }
    elsif ( $base eq 'G' ) { $contador_G++; }
 
    my $diferencia_GC = $contador_G - $contador_C;
 
    push @skew, $diferencia_GC;
 
    if ($min_skew > $diferencia_GC) {
        $min_skew = $diferencia_GC;     # tenemos un nuevo mínimo
        @posiciones_minimo = ();        # reiniciamos el almacén de posiciones
    } 
 
    if ( $diferencia_GC == $min_skew ) {        # si es un mínimo,
        push @posiciones_minimo, $i;            # recordamos su posición
    }
}
 
print join('|', @skew), "\n\n";         # -1|-2|-2|-2|-2|-3|-2|-1|-1|0|1|1|1|1|1|2|1|1|1|2|2|1|0|-1|-1|0|0|0|-1|0|0|0|0|-1|0|-1|-2|-1|-2|-1|-1|-1|-2|-2|-2|-1|-1|-1|-2|-2|-3|-3|-4|-3|-2|-3|-3|-3|-2|-2|-2|-1|0|-1|-1|-1|-1|-1|0|1|1|1|1|1|1|0|-1|0|1|0|1|1|0|-1|0|0|0|0|0|-1|0|0|-1|-2|-2|-3|-4|-3|-3|-2
 
print "Valor mínimo: $min_skew\n";     # Valor mínimo: -4
 
print "@posiciones_minimo\n";           # 52 96
 
Coloreado en 0.002 segundos,  usando GeSHi 1.0.8.4

por **angelluigi** » 2016-11-17 18:39 @818

Muchas gracias, lo acabo de probar y va perfecto.

por **explorer** » 2016-11-19 20:45 @906

Esta es otra versión, más corta, usando otra técnica.

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

#!/usr/bin/env perl
use v5.14;
 
my $secuencia = 'CCTATCGGTGGATTAGCATGTCCCTGTACGTTTCGCCGCGAACTAGTTCACACGGCTTGATGGCAAATGGTTTTTCCGGCGACCGTAATCGTCCACCGAG';
 
my %visto;
my $contador = 0;
 
while ($secuencia =~ m/(?<letras_repetidas>(?|(?<letra>C)+|G+))/g) {
    my $largo = length $+{letras_repetidas};
 
    $contador += $+{letra} eq 'C' ? -$largo : $largo;
 
    push @{$visto{ $contador }}, pos($secuencia);       # posiciones basadas en 0, pero terminan más allá,
                                                        # así que es como si estuvieran basadas en 1
}
 
#use Data::Dumper;
#say Dumper \%visto;
 
my $minimo = (sort {$a <=> $b} keys %visto)[0];
 
say "Mínimo [$minimo] en las posiciones [@{$visto{$minimo}}]";         # Mínimo [-4] en las posiciones [53 97]
Coloreado en 0.001 segundos,  usando GeSHi 1.0.8.4

La idea se basa en localizar los grupos de letras 'C' y 'G', contando las veces que se repiten, y la posición en la que se encuentran.

por **angelluigi** » 2017-01-08 14:10 @632

Buenas, otra vez.

Repasando otra vez el ejercicio me ha venido una duda, en el penúltimo script que me mandaste.

if ($min_skew > $diferencia_GC) {
$min_skew = $diferencia_GC; # tenemos un nuevo mínimo
@posiciones_minimo = (); # reiniciamos el almacén de posiciones
}

No logro entender lo que haces con @posiciones_minimo

Muchas gracias.

por **explorer** » 2017-01-08 16:57 @748

El objetivo es saber en qué posiciones encontramos el mínimo (valor más pequeño) de la diferencia entre 'C' y 'G'.

Entonces, en el if() nos preguntamos si hemos encontrado un valor más pequeño que el que tenemos hasta ahora ($min_skew). Si la $diferencia_GC es menor que $min_skew, resulta que hemos encontrado un mínimo más pequeño, así que todas las @posiciones_minimo que teníamos antes, ya no nos valen (no nos interesan), y por eso reiniciamos el array.

En la línea siguiente hay otro if(), que se cumplirá en el caso de que hayamos encontrado un mínimo más pequeño: en la línea 24 igualamos el valor de las variables, así que la condición de la línea 28 se cumple, por lo que la 29 se ejecuta siempre que encontremos un nuevo mínimo más pequeño. Y la labor de la línea 29 es almacenar la posición de ese nuevo mínimo (de ese, y de todos los demás que sigan iguales a $min_skew).

Foro - Perl en Español

Posiciones de un array

Posiciones de un array

Publicidad

Re: Posiciones de un array

Re: Posiciones de un array

Re: Posiciones de un array

Re: Posiciones de un array

Re: Posiciones de un array

Re: Posiciones de un array

Re: Posiciones de un array

¿Quién está conectado?