版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

xiezhancnu

木蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 5291.5
紅花: 1
帖子: 147
在線: 84.2小時(shí)
蟲(chóng)號(hào): 1387285
注冊(cè): 2011-09-02
專業(yè): 生物信息學(xué)

[求助] perl程序求大神加以修改~~

現(xiàn)在有個(gè)文件夾TOTAL，該文件夾內(nèi)有兩個(gè)文件分別是a.txt和b.txt（見(jiàn)附件），我現(xiàn)在想生成這樣的結(jié)果：
a.txt b.txt（這一行可沒(méi)有）
Os02t0249600-01 6236.420821 19861.20539
Os02t0242600-01 6550.062
Os03t0427300-01 7423.620595 26454.94742
Os01t0762500-00 5592.607026 31007.32272
第二行空格的意思是因?yàn)镺s02t0242600-01在文件a.txt內(nèi)沒(méi)有對(duì)應(yīng)的相，所以就空著！我現(xiàn)在寫(xiě)了一個(gè)程序，如下：
@gbk = glob('TOTAL/*');
for($n=0;$n<=@gbk;$n++){
open(input,"$gbk[$n]"

;
while(my $line={
   my @tmp=split("\t",$line);
   my $key = $tmp[0];
   my $value = $tmp[4];
   if (exists $hash{$key}){
      $hash{$key} = $hash{$key}."\t".$value;
   }else{
      $hash{$key} = $value;
   }
}
}

delete $hash{GeneID};
while(my($key,$value)=each %hash){
print "$key\t$value\n";
}
生成的結(jié)果是這樣的：
Os02t0249600-01 6236.420821 19861.20539
Os02t0242600-01 6550.062
Os03t0427300-01 7423.620595 26454.94742
Os01t0762500-00 5592.607026 31007.32272
也就是無(wú)法生成空格，第二行的6550.062應(yīng)該在19861.20539這一列，請(qǐng)大神在我的程序的基礎(chǔ)上加以修改或者另寫(xiě)程序以實(shí)現(xiàn)我想要的結(jié)果，謝謝了~~~

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : a.txt

2012-07-27 00:05:34, 740 bytes

附件 2 : b.txt

2012-07-27 00:05:48, 562 bytes

» 猜你喜歡

材料/農(nóng)業(yè)專業(yè)，07/08開(kāi)頭均可，過(guò)線就行已經(jīng)有4人回復(fù)
327求調(diào)劑已經(jīng)有5人回復(fù)
361求調(diào)劑已經(jīng)有3人回復(fù)
一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑已經(jīng)有5人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有6人回復(fù)
材料專業(yè)求調(diào)劑已經(jīng)有10人回復(fù)
341求調(diào)劑(一志愿湖南大學(xué)070300) 已經(jīng)有5人回復(fù)
0854 考研調(diào)劑招生了！AI 方向已經(jīng)有15人回復(fù)
384求調(diào)劑已經(jīng)有6人回復(fù)
生物學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

磁性粒子團(tuán)聚成這個(gè)樣子了，求大神指點(diǎn)啊~~ 已經(jīng)有37人回復(fù)
急急急，，急求大神幫忙翻譯一下，跪謝�。�！已經(jīng)有2人回復(fù)
請(qǐng)幫忙編寫(xiě)一個(gè)perl程序處理文本已經(jīng)有26人回復(fù)
求大神幫我看下我這個(gè)方案已經(jīng)有10人回復(fù)
這個(gè)反應(yīng)怎么操作，求做過(guò)此反應(yīng)的大神指點(diǎn) 已經(jīng)有11人回復(fù)
Perl 生物信息學(xué)編程求助，運(yùn)行出現(xiàn)問(wèn)題，煩請(qǐng)幫忙查找一下錯(cuò)誤，謝謝！已經(jīng)有3人回復(fù)
求教一個(gè)origin中的問(wèn)題，求大神指點(diǎn) 已經(jīng)有4人回復(fù)
Perl 編程請(qǐng)大神們幫忙�。�！生物信息的牛人請(qǐng)進(jìn) 已經(jīng)有11人回復(fù)
matlab作業(yè)哪位大神麻煩給做一下,十萬(wàn)火急,周五之間要交已經(jīng)有3人回復(fù)
程序語(yǔ)言版擬開(kāi)新子版“Shell&Perl&Python” 已經(jīng)有10人回復(fù)
【求助】求高人幫我修改程序已經(jīng)有5人回復(fù)

1樓 2012-07-27 00:16:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wizardfan

至尊木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 599 (博士)
貴賓: 1.818
金幣: 24632.2
散金: 197
紅花: 48
沙發(fā): 2
帖子: 2254
在線: 400.7小時(shí)
蟲(chóng)號(hào): 1879241
注冊(cè): 2012-07-05
性別: GG
專業(yè): 生物信息學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xzhdty: 金幣+1, 謝謝參與 2012-07-28 00:01:13
xiezhancnu: 金幣+5 2012-08-11 09:35:09

use strict;
#use Data:

umper;
opendir DIR, "total";
my @files = readdir DIR;
closedir DIR;
my %hash;
#store all data into a big hash
foreach my $file(@files){
next unless ($file=~/\.txt$/);#only deal with text file
&dealOneFile($file);
}
#print Dumper(%hash);

#find the union set of all keys from all files
my %key;
my @doneFiles = keys %hash;
@doneFiles = sort {$a cmp $b} @doneFiles;#sort the file namess
foreach my $tmp(@doneFiles){
my %oneFile = %{$hash{$tmp}};
# print Dumper(%oneFile);
foreach (keys %oneFile){
$key{$_}++
}
}
#print Dumper(%key);

#print header
$"="\t";
print "file\t@doneFiles\n";
#print content
foreach my $key(keys %key){
print "$key";
foreach my $file(@doneFiles){
print "\t";
print "$hash{$file}{$key}" if (exists $hash{$file}{$key});
}
print "\n";
}

sub dealOneFile{
my $filename = "$_[0]";
open IN,"total/$filename";
; #remove header
while(my $line={
chomp($line);
my ($key,$a,$b,$c,$value) = split("\t",$line);
$hash{$_[0]}{$key}=$value;
}
}

贊一下

回復(fù)此樓

2樓2012-07-27 03:12:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhusheng303

木蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3137.3
帖子: 68
在線: 45.6小時(shí)
蟲(chóng)號(hào): 880479
注冊(cè): 2009-10-22
性別: GG
專業(yè): 生物信息學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

有以個(gè)問(wèn)題要問(wèn)你？
你在處理這個(gè)事，最終是不是只有這兩文件？
如果是，解決辦法如下，
while(my $line={
   my @tmp=split("\t",$line);
   my $key = $tmp[0];
   my $value = $tmp[4];
   if (exists $hash{$key}){
      $hash{$key} = $hash{$key}."\t".$value;
   }else{
      $hash{$key} = “". "\t" .$value;  ####或者寫(xiě)成 “$hash{$key} = “-". "\t" .$value”,用“-”表示缺失數(shù)據(jù)
   }

贊一下

回復(fù)此樓

積極面對(duì)

3樓2012-07-27 08:58:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhusheng303

木蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3137.3
帖子: 68
在線: 45.6小時(shí)
蟲(chóng)號(hào): 880479
注冊(cè): 2009-10-22
性別: GG
專業(yè): 生物信息學(xué)

【答案】應(yīng)助回帖

如果超過(guò)3個(gè)文件，解決方法如下，
@gbk = glob('TOTAL/*');
my $hash; ##### 構(gòu)建一個(gè)hash變量
for($n=0;$n<=@gbk;$n++){
      open(input,"$gbk[$n]";
      while(my $line={
               my @tmp=split("\t",$line);
               my $key = $tmp[0];
               my $value = $tmp[4];
               next if（$key eq "GeneID"

;  #### 不記錄“GeneID”這一行
               $hash->{$key}{$n} = $value;  #### 每一個(gè)文件一個(gè)值
   }
}
my $fileNum=$#gbk;
foreach my $key (keys %{$hash}){
         print "$key";
         foreach my $i (0..$fileNum){
                     if(exists $hash->{$key}{$i}){
                              print "\t$$hash->{$key}{$i}";
                     }else{
                              print "\t";
                     }
         }
         print "\n";

}

贊一下

回復(fù)此樓

積極面對(duì)

4樓2012-07-27 09:26:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiezhancnu

木蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 5291.5
紅花: 1
帖子: 147
在線: 84.2小時(shí)
蟲(chóng)號(hào): 1387285
注冊(cè): 2011-09-02
專業(yè): 生物信息學(xué)

引用回帖:

4樓: Originally posted by zhusheng303 at 2012-07-27 09:26:47
如果超過(guò)3個(gè)文件，解決方法如下，
@gbk = glob('TOTAL/*');
my $hash; ##### 構(gòu)建一個(gè)hash變量
for($n=0;$n<=@gbk;$n++){
      open(input,"$gbk";
      while(my $line=< ...

跑不通啊，你試試�。�

回復(fù)此樓

5樓2012-07-27 10:18:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhusheng303

木蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3137.3
帖子: 68
在線: 45.6小時(shí)
蟲(chóng)號(hào): 880479
注冊(cè): 2009-10-22
性別: GG
專業(yè): 生物信息學(xué)

【答案】應(yīng)助回帖

★
xzhdty: 金幣+1, 謝謝 2012-07-28 00:00:57

我寫(xiě)了一個(gè)完整的perl腳本，具體見(jiàn)“parseFile.pl”（http://good.gd/2155638.htm），測(cè)試的結(jié)果見(jiàn)“結(jié)果.bmp”（http://good.gd/2155640.htm）
"第二行空格的意思是因?yàn)镺s02t0242600-01在文件a.txt內(nèi)沒(méi)有對(duì)應(yīng)的相",你說(shuō)a文件中沒(méi)有“Os02t0242600-01”這行，但是，其實(shí)是有的，見(jiàn)圖“文件a.bmp”（http://good.gd/2155643.htm）。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

xiezhancnu

積極面對(duì)

6樓2012-07-27 13:38:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wizardfan

至尊木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 599 (博士)
貴賓: 1.818
金幣: 24632.2
散金: 197
紅花: 48
沙發(fā): 2
帖子: 2254
在線: 400.7小時(shí)
蟲(chóng)號(hào): 1879241
注冊(cè): 2012-07-05
性別: GG
專業(yè): 生物信息學(xué)

【答案】應(yīng)助回帖

引用回帖:

5樓: Originally posted by xiezhancnu at 2012-07-27 10:18:20
跑不通啊，你試試！！...

Have you tried my codes which can also allow multiple file processing? The commented-out print statements allow you to check the middle result.

贊一下

回復(fù)此樓

7樓2012-07-27 23:28:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiezhancnu

木蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 5291.5
紅花: 1
帖子: 147
在線: 84.2小時(shí)
蟲(chóng)號(hào): 1387285
注冊(cè): 2011-09-02
專業(yè): 生物信息學(xué)

引用回帖:

3樓: Originally posted by zhusheng303 at 2012-07-27 08:58:57
有以個(gè)問(wèn)題要問(wèn)你？
你在處理這個(gè)事，最終是不是只有這兩文件？
如果是，解決辦法如下，
while(my $line=<input>{
   my @tmp=split("\t",$line);
   my $key = $tmp;
   my $valu ...

有很多文件的！

回復(fù)此樓

8樓2012-07-28 22:09:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiezhancnu

木蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 5291.5
紅花: 1
帖子: 147
在線: 84.2小時(shí)
蟲(chóng)號(hào): 1387285
注冊(cè): 2011-09-02
專業(yè): 生物信息學(xué)

送鮮花一朵

引用回帖:

6樓: Originally posted by zhusheng303 at 2012-07-27 13:38:37
我寫(xiě)了一個(gè)完整的perl腳本，具體見(jiàn)“parseFile.pl”（http://good.gd/2155638.htm），測(cè)試的結(jié)果見(jiàn)“結(jié)果.bmp”（http://good.gd/2155640.htm）
"第二行空格的意思是因?yàn)镺s02t0242600-01在文件a.txt內(nèi)沒(méi)有對(duì)應(yīng) ...

你好，呵呵，我之前把文件弄錯(cuò)了，呵呵，a文件，b文件的名字弄顛倒了，呵呵，另外a、b文件只是測(cè)試文件，總共有七個(gè)文件的，呵呵，我還沒(méi)測(cè)試一下您的程序是否能跑通，送您鮮花一朵吧，呵呵，這樣的話生成結(jié)果是這樣的：
Os02t0249600-01       6236.420821       19861.20539
Os02t0242600-01                                  6550.062
Os03t0427300-01       7423.620595       26454.94742
Os01t0762500-00       5592.607026       31007.32272
這樣也許您就不覺(jué)得奇怪了吧！

贊一下

回復(fù)此樓

9樓2012-07-28 22:28:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 xiezhancnu 的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 265求調(diào)劑 +10	梁梁校校 2026-03-17	10/500	2026-03-23 21:17 by 一切OK
[考研] 化學(xué)308分求調(diào)劑 +3	你好明天你好 2026-03-23	3/150	2026-03-23 20:11 by macy2011
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +3	SYA！ 2026-03-23	3/150	2026-03-23 19:09 by macy2011
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-17	10/500	2026-03-23 13:05 by 醉在風(fēng)里
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +8	材料學(xué)257求調(diào)劑 2026-03-23	9/450	2026-03-23 13:01 by ztnimte
[考研] 276求調(diào)劑 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在風(fēng)里
[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋?xiě)?/a>
[考研] 尋找調(diào)劑 +4	倔強(qiáng)芒? 2026-03-21	4/200	2026-03-22 16:14 by 木托莫露露
[考研] 311求調(diào)劑 +3	26研0 2026-03-20	3/150	2026-03-22 14:46 by ColorlessPI
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 268求調(diào)劑 +9	簡(jiǎn)單點(diǎn)0 2026-03-17	9/450	2026-03-21 15:37 by lature00
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 求調(diào)劑 +3	白QF 2026-03-21	3/150	2026-03-21 13:12 by zhukairuo
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽(yáng)w
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx

24小時(shí)熱門版塊排行榜

[求助] perl程序求大神加以修改~~

» 本帖附件資源列表

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

» 本帖已獲得的紅花（最新10朵）

【答案】應(yīng)助回帖

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: