| 查看: 589 | 回复: 2 | ||
xiezhancnu木虫 (小有名气)
|
[求助]
另开新帖求perl程序高手指导~~
|
|
http://muchong.com/bbs/viewthread.php?tid=4816445 就是这个，我需要一个同时处理3000多个基因的程序，具体内容请点击链接进入另一个帖子，我急用的，呵呵~~~求求求高手赐教， PS：其实说真的，我perl程序的引用真心不会用，每次一牵扯这种二维数组类型的数据就写不出来，很悲催~~~.,,真心求帮助~~先谢过~~呵呵~~~ |
回复此楼» 猜你喜欢
一志愿华东理工大学081700，初试分数271
已经有5人回复
-
081700 调剂 267分
已经有7人回复
-
材料调剂
已经有3人回复
-
一志愿河北工业大学0817化工278分求调剂
已经有10人回复
-
085600材料与化工调剂
已经有12人回复
-
材料292调剂
已经有4人回复
-
291求调剂
已经有8人回复
-
材料专硕找调剂
已经有3人回复
-
材料专业求调剂
已经有11人回复
-
一志愿国科过程所081700，274求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- perl程序继续中，脸红贴上自写程序，求拍砖。
已经有4人回复
- 请帮忙编写一个perl程序处理文本
已经有26人回复
- Perl 生物信息学编程求助，运行出现问题，烦请帮忙查找一下错误，谢谢！
已经有3人回复
- 求高手指点witting反应
已经有14人回复
- 程序语言版拟开新子版“Shell&Perl&Python”
已经有10人回复
- 【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
- 【求助】一个没想清楚的问题—纳米粒子和量子点的区别
已经有21人回复
- 【求助】电化学阻抗谱低频感抗弧不是半圆趋势 该怎样拟合
已经有22人回复
- 【求助】求助高手帮忙
已经有8人回复
- 【求助】有很多错误的程序，找不出来了
已经有18人回复
- 【求助】请问下，河北省每年公务员招考时几月份？
已经有8人回复
- 【求助】Linux 安装
已经有17人回复
- 【求助】帮我编个程啊~C语言，C++，Matlab，VB...都行~
已经有14人回复
- 【求助】大家能不能推荐几个国外做基因转染（药物控释）的牛人？（呵呵，回帖就送~）
已经有39人回复
wizardfan
至尊木虫 (著名写手)
- 应助: 599 (博士)
- 贵宾: 1.818
- 金币: 24632.2
- 散金: 197
- 红花: 48
- 沙发: 2
- 帖子: 2254
- 在线: 400.7小时
- 虫号: 1879241
- 注册: 2012-07-05
- 性别: GG
- 专业: 生物信息学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiezhancnu: 金币+30 2012-08-14 13:11:13
感谢参与，应助指数 +1
xiezhancnu: 金币+30 2012-08-14 13:11:13
|
use strict; use Data:  umper;#create index for names to avoid the situation that two files have different order of gene names my %idx; my %idxReverse; open IN,"pvalue.csv"; my $line = chomp($line); my @arr = split(",",$line); my $len = scalar @arr; for(my $i=1;$i<$len;$i++){ $idx{$arr[$i]}=$i; $idxReverse{$i}=$arr[$i]; } close IN; #populate the data structure my %hash;#values my @tags;#the elements to be displayed in the header line after gene1 and gene2 &dealOneFile("cor.csv","correlation"  ;&dealOneFile("pvalue.csv","pvalue" ;#print the result #header print "gene1,gene2"; foreach my $tag(@tags){ print ",$tag"; } print "\n"; #values foreach my $aa(sort {$a<=>$b} keys %hash){ my %tmp = %{$hash{$aa}}; foreach my $bb(sort {$a<=>$b} keys %tmp){ print "$idxReverse{$aa},$idxReverse{$bb}"; foreach my $tag(@tags){ print ",$hash{$aa}{$bb}{$tag}"; } print "\n"; } } sub dealOneFile{ my $filename = $_[0]; my $tag = $_[1]; push(@tags,$tag); open IN,"$filename"; my $line = chomp($line); my @header = split(",",$line); shift @header; my $count = 0; while($line= chomp($line); $count++; my @arr = split(",",$line); my $current = shift(@arr); my $currentIdx = $idx{$current}; for(my $i=$count;$i<$len;$i++){ my $second = $header[$i]; my $secondIdx = $idx{$second}; if ($currentIdx<$secondIdx){ $hash{$currentIdx}{$secondIdx}{$tag}=$arr[$i]; }else{ $hash{$secondIdx}{$currentIdx}{$tag}=$arr[$i]; } } } } 不知道什么原因，你老是把文件名搞错，自己认真点吧。Perl是很基础的生物信息学工具，好好弄懂，别老是一不会就来提问了。自己的本领才是最关键的。 |
2楼2012-08-13 06:02:26
wizardfan
至尊木虫 (著名写手)
- 应助: 599 (博士)
- 贵宾: 1.818
- 金币: 24632.2
- 散金: 197
- 红花: 48
- 沙发: 2
- 帖子: 2254
- 在线: 400.7小时
- 虫号: 1879241
- 注册: 2012-07-05
- 性别: GG
- 专业: 生物信息学
★
xzhdty: 金币+1, 谢谢 2012-08-13 22:23:33
xzhdty: 金币+1, 谢谢 2012-08-13 22:23:33
|
这个程序可以处理任意多个基因，任意多个数据文件（比如还有个文件叫similarity),只要内存够。但是有一个要求，输入格式必须是CSV，comma separated values。 例如 ,DR_2577,DR_A0328,DR_2258,DR_0109,DR_0465,DR_0539,DR_1356,DR_2207,DR_2121,DR_1962 DR_2577,NA,1.71E-08,0.000439,0.026534,1.15E-06,0.712607,0.052184,7.77E-17,2.91E-10,4.74E-17 DR_A0328,NA,NA,0.011085,0.067408,4.92E-07,0.535037,0.019445,4.43E-09,1.27E-07,1.07E-07 DR_2258,NA,NA,NA,0.824648,0.00147,0.054368,0.016167,0.003185,0.000694,0.00236 DR_0109,NA,NA,NA,NA,0.09646,0.161788,0.696339,0.004726,0.007971,0.040283 DR_0465,NA,NA,NA,NA,NA,0.189011,0.012672,1.10E-06,4.67E-06,0.000102 DR_0539,NA,NA,NA,NA,NA,NA,0.034369,0.57409,0.306483,0.573561 DR_1356,NA,NA,NA,NA,NA,NA,NA,0.217671,0.024297,0.381359 DR_2207,NA,NA,NA,NA,NA,NA,NA,NA,6.84E-10,5.93E-13 DR_2121,NA,NA,NA,NA,NA,NA,NA,NA,NA,1.38E-09 DR_1962,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA |
3楼2012-08-13 06:05:12
