如何搭配使用PERL

Question

我有這樣

pop A B C D E 
P1 T/T C/C C/C T/T C/C 
P2 A/A G/G C/C T/T C/C 
1 A/A G/G C/C T/T C/C 
2 A/A G/G C/C T/T C/C 
3 A/T A/C A/G A/T A/C 
4 T/A T/G T/C T/A T/G 
5 G/A G/T G/C G/A G/T 
6 C/A C/T C/G C/A C/T 
pop A B C D E 
P1 T/T C/C C/C T/T C/C 
P2 A/A G/G C/C T/T C/C 
1 A/A G/G C/C T/T C/C 
2 A/A G/G C/C T/T C/C 
3 A/T A/C A/G A/T A/C 
4 T/A T/G T/C T/A T/G 
5 G/A G/T G/C G/A G/T 
6 C/A C/T C/G C/A C/T 

我想要做我的數據

1. 我想將所有A/A到A的folowwing事情數據列和計算比例， T/T至T,C/C至C, G/G至G,Z/Z至-和-/-至-級其餘的人物，如A/T，G/T，C/G，T/C到H

2. 現在我想通過比較P1與P2知道從A到E狀態，如果P1=P2則狀態從A到E是單或P1任何一個或P2包含Z/Z或-/-然後從狀態到A是E其他單從狀態到A是E聚

3. 我想在pop柱P2在彈出的列匹配1 A至E，如果1 pop列pop列匹配p2，它的地位是聚才，我想以其他方式給予1正因爲如此，如果它是單聲道我不想做任何事情。

4. 現在我將計算# 1 S和# H的

5. 最後，我將計算%sim這個公式=((#1*2+#H)/((#1+#H)*2))*100。

6. 我要重複相同的步驟第二套父母P1和P2

進出料會是這樣

pop A B C D E A B C D E   
P1 POLY POLY MONO MONO MONO POLY POLY MONO MONO MONO #1's #H's %sim 
P2 A G C T C        
1 A G C T C 1 1 C T C 2 0 100 
2 A G C T C 1 1 C T C 2 0 100 
3 A G C - C H H H H H 0 5 100 
4 H H H H H H H H H H 0 5 50 
5 H H H H H H H H H H 0 5 50 
6 H H H H H H H H H H 0 5 50 
7 H H H H H H H H H H 0 5 50 

現在我成功了，直到第1步，但不知道如何進一步進行，這是我直到現在的代碼

#!/usr/bin/perl -w 
use strict; 
open(FILE, "<input.txt") || die "File not found"; 
my @lines = <FILE>; 
my @newlines; 
foreach(@lines) { 
    $_ =~ s/AA/A/g; 
    $_ =~ s/TT/T/g; 
    $_ =~ s/GG/G/g; 
    $_ =~ s/CC/C/g; 
    $_ =~ s/AT/H/g; 
    $_ =~ s/AG/H/g; 
    $_ =~ s/AC/H/g; 
    $_ =~ s/TA/H/g; 
    $_ =~ s/TG/H/g; 
    $_ =~ s/TC/H/g; 
    $_ =~ s/GA/H/g; 
    $_ =~ s/GT/H/g; 
    $_ =~ s/GC/H/g; 
    $_ =~ s/CA/H/g; 
    $_ =~ s/CT/H/g; 
    $_ =~ s/CG/H/g; 
    $_ =~ s/ZZ/-/g; 

    push(@newlines,$_); 
} 
open(FILE, ">input1.txt") || die "File not found"; 
print FILE @newlines; 
close(FILE); 

Answer 1

這會替換並創建前兩行。我不明白說明書的其餘部分（其中的7號線是從哪裏來的？）

#!/usr/bin/perl 
use warnings; 
use strict; 
use feature qw{ say }; 

use Data::Dumper; 


*ARGV = *DATA{IO} unless @ARGV; 

my (@parents, @rows); 

sub { 
    my $header = <>; 
    push @parents, map [ split ' ', <> ], 1, 2; 
    push @rows, map [ split ' ', <> ], 1 .. 6; 
}->() for 1, 2; 

for (map @$_, @parents, @rows) { 
    s= ([ACTG])/\1 =$1=x; 
    s= ([-Z]) /\1 =-=x; 
    s= .  /. =H=x; 
} 

say join "\t", 'pop', ('A' .. 'E') x 2; 

print 'P1'; 
for my $parent (0, 1) { 
    print join "\t", q(), map { 
     my $p1 = $parents[ $parent * 2 ][$_]; 
     my $p2 = $parents[ 1 + $parent * 2 ][$_]; 
     ($p1 eq $p2 or '-' eq $p1 or '-' eq $p2) ? 'mono' : 'poly'; 
    } 1 .. 5; 
} 
print "\n"; 


__DATA__ 
pop A B C D E 
P1 T/T C/C C/C T/T C/C 
P2 A/A G/G C/C T/T C/C 
1 A/A G/G C/C T/T C/C 
2 A/A G/G C/C T/T C/C 
3 A/T A/C A/G A/T A/C 
4 T/A T/G T/C T/A T/G 
5 G/A G/T G/C G/A G/T 
6 C/A C/T C/G C/A C/T 
pop A B C D E 
P1 T/T C/C C/C T/T C/C 
P2 A/A G/G C/C T/T C/C 
1 A/A G/G C/C T/T C/C 
2 A/A G/G C/C T/T C/C 
3 A/T A/C A/G A/T A/C 
4 T/A T/G T/C T/A T/G 
5 G/A G/T G/C G/A G/T 
6 C/A C/T C/G C/A C/T