當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

xiezhancnu

木蟲 (小有名氣)

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

[求助] perl程序求助，請(qǐng)大神們賜教~~~

我現(xiàn)在有十個(gè)基因，附件里的兩個(gè)文件分別是這十個(gè)基因任意兩個(gè)之間的相關(guān)系數(shù)和p-alue的值，附件應(yīng)該是很容易看懂的，呵呵~~，現(xiàn)在呢，我想改變一下輸出的格式，將兩個(gè)文件整合一下，輸出成這種格式：
基因1 基因2 相關(guān)系數(shù) p-value （這一行可有可無）
DR_2577 DR_A0328 0.81186 0.0000000171
……
輸出結(jié)果不算可有可無的第一行的話，應(yīng)該總共有45行，就是這么一個(gè)程序~~
我覺得應(yīng)該是說清楚了，呵呵~~若是哪里說的不夠清楚的話，請(qǐng)大神們主動(dòng)提出來，呵呵~~謝謝啦~~~

回復(fù)此樓

» 本帖附件資源列表

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

2012-08-10 21:25:34, 11.5 K

附件 2 : cor.xls

2012-08-10 21:25:40, 12 K

» 猜你喜歡

305分求調(diào)劑（食品工程）已經(jīng)有7人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有5人回復(fù)
299求調(diào)劑已經(jīng)有6人回復(fù)
材料292調(diào)劑已經(jīng)有5人回復(fù)
085404電子信息284分求調(diào)劑已經(jīng)有4人回復(fù)
收08調(diào)劑生已經(jīng)有3人回復(fù)
081700 調(diào)劑 267分已經(jīng)有9人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
086003食品工程求調(diào)劑已經(jīng)有4人回復(fù)
0856材料專碩353求調(diào)劑已經(jīng)有5人回復(fù)

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

請(qǐng)教碳材料的氮?dú)馕矫摳角€分析！請(qǐng)各位大神不吝賜教！~~ 已經(jīng)有5人回復(fù)
請(qǐng)幫忙編寫一個(gè)perl程序處理文本已經(jīng)有26人回復(fù)
Perl 生物信息學(xué)編程求助，運(yùn)行出現(xiàn)問題，煩請(qǐng)幫忙查找一下錯(cuò)誤，謝謝！已經(jīng)有3人回復(fù)
【緊急求助】offer鑒定及相關(guān)問題~~請(qǐng)有經(jīng)驗(yàn)的蟲友不吝賜教~~ 已經(jīng)有10人回復(fù)
Perl 編程請(qǐng)大神們幫忙！��！生物信息的牛人請(qǐng)進(jìn) 已經(jīng)有11人回復(fù)
【求助/交流】求構(gòu)載體高手賜教構(gòu)建亞細(xì)胞定位載體pCAMBIA1302 已經(jīng)有4人回復(fù)
【求助】新手上路請(qǐng)不吝賜教（熒光法測(cè)臨界膠束濃度原理）已經(jīng)有4人回復(fù)
【求助】扣式電池組裝的幾個(gè)基本問題，請(qǐng)高手賜教已經(jīng)有8人回復(fù)
【求助】關(guān)于催化燃燒，請(qǐng)蟲友賜教。已經(jīng)有7人回復(fù)
【求助/交流】求教一個(gè)關(guān)于酶方面的問題，請(qǐng)高手不吝賜教已經(jīng)有10人回復(fù)
【求助】求大俠賜教：非水溶液滴定已經(jīng)有3人回復(fù)
【求助】熟悉XRD、XRF的大俠來賜教已經(jīng)有6人回復(fù)

1樓 2012-08-10 21:39:31

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

jackie1179

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1929.2
紅花: 1
帖子: 460
在線: 47.2小時(shí)
蟲號(hào): 1119284
注冊(cè): 2010-10-11
專業(yè): 生物信息學(xué)

★
xzhdty: 金幣+1, 謝謝 2012-08-13 22:22:35

借鑒了wizardfan在你新帖中的split方法，對(duì)原程序稍加修改，或許可以滿足你的要求。

CODE:

#!/usr/bin/perl 

use warnings;

$cor_file = "cor.csv";

$pvalue_file = "pvalue.csv";

$out = "out";

undef %cor;

undef %pvalue;

open(COR,"$cor_file") || die "Cannot open this file$!";

$tmp= ;

chomp($tmp);

$line =~ s/^\t//;

@gene = split(/\t/,$tmp);

print @gene,"\n";

$cnt = 0;

while() {

  $cnt++;

  chomp;

  @array = split(/\t/,$_);

  for($j = $cnt; $j < 10 ; $j++ ) {

    $cor{$gene[$cnt - 1],$gene[$j]} = $array[$j + 1];

  }

}

close COR;

foreach (sort keys %cor) {

  print $_,"=>","$cor{$_}\n";

}

open(PVALUE,"$pvalue_file") || die "Cannot open this file$!";

;

$cnt = 0;

while() {

  $cnt++;

  chomp;

  @array = split(/\t/,$_);

  for($j = $cnt; $j < 10; $j++) {

    $pvalue{$gene[$cnt - 1],$gene[$j]} = $array[$j + 1];

  }

}

close PVALUE;

open(OUT,">$out") || die "Cannot write to this file$!";

print OUT "gene1\tgene2\tcorrelation efficient\t p-value\n";

for($i = 0; $i < 10; $i++) {

  for($j = 0; $j < 10; $j++) {

    if(exists $cor{$gene[$i],$gene[$j]}) {

      printf OUT "%-8s\t%-8s\t%8s\t%8s\n",$gene[$i],$gene[$j],$cor{$gene[$i],$gene[$j]},$pvalue{$gene[$i],$gene[$j]};

    }

  }

}

close OUT;

贊一下

回復(fù)此樓

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

xiezhancnu

9樓2012-08-13 10:01:38

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

查看全部 10 個(gè)回答

wizardfan

至尊木蟲 (著名寫手)

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

本來是想幫你的，但是以前幫過你以后，一點(diǎn)回音都沒有http://www.gaoyang168.com/bbs/viewthread.php?tid=4773466&fpage=1，我也建議別人不要幫助沒有基本禮儀的人。

贊一下

回復(fù)此樓

2樓2012-08-11 01:58:24

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

xiezhancnu

木蟲 (小有名氣)

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

引用回帖:

2樓: Originally posted by wizardfan at 2012-08-11 01:58:24
本來是想幫你的，但是以前幫過你以后，一點(diǎn)回音都沒有http://www.gaoyang168.com/bbs/viewthread.php?tid=4773466&fpage=1，我也建議別人不要幫助沒有基本禮儀的人。

不是的，我暈~~，現(xiàn)在放假了，那幾天家里有事回去了，我家農(nóng)村的，家里又上不了網(wǎng)！在家?guī)Я耸畞硖�，就把那個(gè)帖子給忘了，平時(shí)小木蟲上的又不多，請(qǐng)您不要生氣了哈~~現(xiàn)在把金幣給您~~暈~~·

贊一下

回復(fù)此樓

3樓2012-08-11 09:30:25

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

jackie1179

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1929.2
紅花: 1
帖子: 460
在線: 47.2小時(shí)
蟲號(hào): 1119284
注冊(cè): 2010-10-11
專業(yè): 生物信息學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xiezhancnu: 金幣+10 2012-08-11 09:45:21
xzhdty: 金幣+1, 謝謝參與 2012-08-13 22:21:08

你的原始文件要保存為文本文檔才好操作（用制表符分隔），因?yàn)橛胮erl操作excel表格我也不會(huì)

CODE:

#!/usr/bin/perl 

use warnings;

$cor_file = "cor.csv";

$pvalue_file = "pvalue.csv";

$out = "out";

undef %cor;

undef %pvalue;

open(COR,"$cor_file") || die "Cannot open this file$!";

 =~ /(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)/;

for($i = 1; $i <= 10; $i++) {

  push @gene,$$i;

}

print @gene,"\n";

$cnt = 0;

while() {

  $cnt++;

  $_ =~ /^\S+\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)/;

  for($j = $cnt; $j < 10 ; $j++ ) {

    $cor{$gene[$cnt - 1],$gene[$j]} = ${$j + 1};

  }

}

close COR;

foreach (sort keys %cor) {

  print $_,"=>","$cor{$_}\n";

}

open(PVALUE,"$pvalue_file") || die "Cannot open this file$!";

;

$cnt = 0;

while() {

  $cnt++;

  $_ =~ /^\S+\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)/;

  for($j = $cnt; $j < 10; $j++) {

    $pvalue{$gene[$cnt - 1],$gene[$j]} = ${$j + 1};

  }

}

close PVALUE;

open(OUT,">$out") || die "Cannot write to this file$!";

print OUT "gene1\tgene2\tcorrelation efficient\t p-value\n";

for($i = 0; $i < 10; $i++) {

  for($j = 0; $j < 10; $j++) {

    if(exists $cor{$gene[$i],$gene[$j]}) {

      printf OUT "%-8s\t%-8s\t%8s\t%8s\n",$gene[$i],$gene[$j],$cor{$gene[$i],$gene[$j]},$pvalue{$gene[$i],$gene[$j]};

    }

  }

}

close OUT;

贊一下

回復(fù)此樓

4樓2012-08-11 09:32:35

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

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

[考研] 07化學(xué)280分求調(diào)劑 +4	722865 2026-03-23	4/200	2026-03-24 00:01 by chixmc
[考研] 085600材料與化工調(diào)劑 +7	A-哆啦Z夢(mèng) 2026-03-23	12/600	2026-03-23 23:16 by 星空星月
[考研] 291求調(diào)劑 +8	hhhhxn.. 2026-03-23	8/400	2026-03-23 23:15 by peike
[考研] 0703化學(xué)求調(diào)劑 +4	奶油草莓. 2026-03-22	5/250	2026-03-23 19:37 by pswait
[考研] 333求調(diào)劑 +6	87639 2026-03-21	10/500	2026-03-23 10:41 by Iveryant
[考研] 323求調(diào)劑 +6	洼小桶 2026-03-18	6/300	2026-03-23 00:29 by king123！
[考研] 352求調(diào)劑 +3	大米飯！ 2026-03-22	3/150	2026-03-22 23:28 by king123！
[考研] 317求調(diào)劑 +12	申子申申 2026-03-19	18/900	2026-03-22 22:23 by luoyongfeng
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 求調(diào)劑院校信息 +6	CX 330 2026-03-21	6/300	2026-03-22 15:25 by 無懈可擊111
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	4/200	2026-03-22 15:12 by 朱芷琳
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 317求調(diào)劑 +5	申子申申 2026-03-19	9/450	2026-03-20 22:26 by JourneyLucky
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu