有搞bioperl的吗？求一个计算矩阵的的程序（求代码）！紧急

aids260

能不能写个程序把每列上的数字计算个平均值输出（保留到小数点后2位就可以了了）如图所示

输出格式：


>pro521/BGIBMGA123456-TA
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx

aids260

顶，达人，赶快来！

kingwmj

这个不用bioperl，bioperl也很可能没这个功能，因为太简单了。你把要读取的数据读到二维数组里面，然后求平均值就行了。

aids260

回复 3# kingwmj


    用一般的perl 该可以吧，请指点哦！

kingwmj

回复  kingwmj


    用一般的perl 该可以吧，请指点哦！
aids260 发表于 2011-05-08 09:01

    可以

aids260

回复 5# kingwmj


    麻烦帮给弄个代码哦！ 本身水平有限，又搞论文实在是没有时间啊！

aids260

顶一个！！

yinyuemi

回复 7# aids260

1. awk '!/>/{++p;for(i=1;i<=NF;i++) a[i]+=$i}END{for(i=1;i<NF;i++) printf a[i]/p FS;print a[NF]}' file

复制代码

perl:

1. perl -lane 'if(/^[0-9]/){++$p;for $a (0..$#F) {$sum[$a]+=$F[$a]}}END{for $b (0..$#F){$ave=$sum[$b]/$p;printf "%.2f ","$ave"};print ""}' file
   

复制代码

aids260

回复 8# yinyuemi


    大虾，你的perl line能在windows下用吗？？？

yinyuemi

回复 9# aids260


    不大清楚，如果不行的话，把单引号改成双引号，试试