我對編碼相當缺乏經驗,但我經常使用Perl合併文件並在兩個文件之間匹配ID和信息。我剛剛嘗試過使用我之前多次使用過的程序來匹配兩個文件,但這次它不起作用,我不明白爲什麼。合併兩個文件的腳本
下面是代碼:
use strict;
use warnings;
use vars qw($damID $damF $damAHC $prog $hash1 %hash1 $info1 $ID $sire $dam $F $FB $AHC $FA $hash2 %hash2 $info2);
open (FILE1, "<damF.txt") || die "$!\n Couldn't open damF.txt\n";
my $N = 1;
while (<FILE1>){
chomp (my $line=$_);
next if 1..$N==$.;
my ($damID, $damF, $damAHC, $prog) = split (/\t/, $line);
if ($prog){
$hash1 -> {$prog} -> {info1} = "$damID\t$damF\t$damAHC";
}
open (FILE2, "<whole pedigree_F.txt") || die "$!\n whole pedigree_F.txt \n";
open (Output, ">Output.txt")||die "Can't Open Output file";
while (<FILE2>){
chomp (my $line=$_);
next if 1..$N==$.;
my ($ID, $sire, $dam, $F, $FB, $AHC, $FA) = split (/\t/, $line);
if ($ID){
$hash2 -> {$ID} -> {info2} = "$F\t$AHC";
}
if ($ID && ($hash1->{$prog})){
$info1 = $hash1 -> {$prog} -> {info1};
$info2 = $hash2 -> {$ID} -> {info2};
print "$ID\t$info2\t$info1\n";
}
}
}
close(FILE1);
close FILE2;
close Output;
print "Done!\n";
,並從兩個輸入文件格式這些片段:
文件1:
501093 0 0 3162
2958 0 0 3163
1895 0 0 3164
1382 0 0 3165
2869 0 0 3166
2361 0 0 3167
754 0 0 3168
3163 0 0 3169
文件2:
49327 20543 49325 0.077 0.4899 0.808 0.0484
49328 15247 49326 0.0755 0.5232 0.8972 0.0499
49329 27823 49327 0.0834 0.5138 0.8738 0.0541
我想匹配column 4 in file 1
與column 1 in file 2
的值。
然後我還想打印從columns 2 and 3 in file 1
和columns 3 and 5 in file 2
的匹配值。
此外,它可能是值得一提的每個文件上有大約條目。
這是我得到的輸出:
11476 0.0362 0.3237 501093 0 0
11477 0.0673 0.4768 501093 0 0
11478 0.0443 0.2619 501093 0 0
注意,它不是通過我創建的第一個散列循環。
怎麼回事?第二個文件名中的空間是否有可能拋棄?你有什麼錯誤嗎? – Ilion
請添加您想要的輸出 – dawg
第一個散列似乎沒有循環。它幾乎看起來像使用vars qw();沒有在文件上工作 – user7675703