版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調劑公告

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

xiezhancnu

木蟲 (小有名氣)

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

[求助] perl程序求助，請大神們賜教~~~

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

回復此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內容負責。
本內容由用戶自主發(fā)布，如果其內容涉及到知識產(chǎn)權問題，其責任在于用戶本人，如對版權有異議，請聯(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

» 猜你喜歡

085601材料工程找調劑已經(jīng)有11人回復
總分293求調劑已經(jīng)有9人回復
334分一志愿武理材料求調劑已經(jīng)有5人回復
327求調劑已經(jīng)有6人回復
329求調劑已經(jīng)有10人回復
275求調劑已經(jīng)有15人回復
291求調劑已經(jīng)有6人回復
309求調劑已經(jīng)有7人回復
297求調劑已經(jīng)有11人回復
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有6人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

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

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

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

jackie1179

木蟲 (正式寫手)

應助: 11 (小學生)
金幣: 1929.2
紅花: 1
帖子: 460
在線: 47.2小時
蟲號: 1119284
注冊: 2010-10-11
專業(yè): 生物信息學

【答案】應助回帖

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

你的原始文件要保存為文本文檔才好操作（用制表符分隔），因為用perl操作excel表格我也不會

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;

贊一下

回復此樓

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

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

查看全部 10 個回答

wizardfan

至尊木蟲 (著名寫手)

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

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

贊一下

回復此樓

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

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

xiezhancnu

木蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

xiezhancnu

木蟲 (小有名氣)

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

引用回帖:

4樓: Originally posted by jackie1179 at 2012-08-11 09:32:35
你的原始文件要保存為文本文檔才好操作（用制表符分隔），因為用perl操作excel表格我也不會

#!/usr/bin/perl

use warnings;

$cor_file = "cor.csv";
$pvalue_file = "pvalue.csv";
...

你好，剛剛看了一下您的程序，發(fā)現(xiàn)您的程序采用了正則表達式捕獲變量的方法，但是實際要處理的基因并不是10個，是3000多個，這樣的話，您的程序就完全不是用了，/^\S+\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)\t(\S+)/總不能延長到3000多個吧，不知道您可不可以寫一個適用于任意數(shù)目基因的程序？謝謝了~~

贊一下

回復此樓

5樓2012-08-11 15:00:36

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

查看全部 10 個回答

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

[考研] 291求調劑 +5	Y-cap 2026-03-29	6/300	2026-03-29 13:18 by mumin1990
[考研] 求調劑，一志愿南京航空航天大學，080500材料科學與工程學碩，總分289分 +7	@taotao 2026-03-29	7/350	2026-03-29 12:03 by longlotian
[考研] 0703化學調劑，求導師收 +9	天天好運來上岸?/a> 2026-03-24	10/500	2026-03-28 22:17 by chemzp
[考研] 315求調劑 +4	akie... 2026-03-28	5/250	2026-03-28 21:05 by zhq0425
[考研] 一志愿北京工業(yè)大學，324分求調劑 +4	零八# 2026-03-28	4/200	2026-03-28 15:01 by 17865157980
[考研] 求調劑 +3	QiMing7 2026-03-25	4/200	2026-03-28 14:30 by QiMing7
[考研] 一志愿南昌大學324求調劑 +7	hanamiko 2026-03-27	7/350	2026-03-28 09:56 by 李上岸0921
[考研] 085600，材料與化工321分調劑 +4	大饞小子 2026-03-27	6/300	2026-03-27 14:11 by 松花缸1201
[考研] 333求調劑 +3	question挽風 2026-03-23	3/150	2026-03-27 11:29 by 不吃魚的貓
[考研] 材料求調劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 0703化學338求調劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 359求調劑 +4	王了個楠 2026-03-25	4/200	2026-03-27 08:43 by 不吃魚的貓
[考研] 求調劑，一志愿南京航空航天大學大學，080500材料科學與工程學碩 +4	@taotao 2026-03-26	5/250	2026-03-27 08:10 by hypershenger
[考研] 081200-11408-276學碩求調劑 +4	崔wj 2026-03-26	4/200	2026-03-27 08:04 by chemisry
[考研] 材料專碩 335 分求調劑 +4	拒絕冷暴力 2026-03-25	4/200	2026-03-25 18:45 by haxia
[考研] 0854電子信息求調劑 +7	α____ 2026-03-22	9/450	2026-03-25 13:37 by α____
[考研] 0854人工智能方向招收調劑 +4	章小魚567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 材料調劑 +3	iwinso 2026-03-23	3/150	2026-03-25 11:29 by greychen00
[考研] 一志愿河北工業(yè)大學0817化工278分求調劑 +7	jhybd 2026-03-23	12/600	2026-03-24 09:03 by jhybd
[考研] 284求調劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

xiezhancnu

[求助] perl程序求助，請大神們賜教~~~

» 本帖附件資源列表

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

jackie1179

【答案】應助回帖

wizardfan

xiezhancnu

xiezhancnu

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] perl程序求助，請大神們賜教~~~

» 本主題相關價值貼推薦，對您同樣有幫助: