前提:文件名不同、文件的行名相同
结果:行名<-行名,列名<-文件名
#首先新建文件tmp.sh 输入这个代码:
perl -le '{print"gene_$_\t".int(rand(1000))foreach 1..99}'
## 然后用perl脚本调用这个tmp.sh文件:
perl -e 'system(" bash tmp.sh >$_.txt")foreach a..z'
## 首先在GSE48213_RAW目录里面生成tmp.txt文件
awk '{print FILENAME"\t"$0}' * |grep -v EnsEMBL_Gene_ID >tmp.txt
## 然后把tmp.txt导入R语言里面用reshape2处理即可!
setwd("D:/Desktop/bind/")
a=read.table('tmp.txt',sep = '\t',stringsAsFactors = F)
library(reshape2)
fpkm <- dcast(a,formula = V2~V1)
+++++++++++++++++++
paste file1 file2 file3 | awk '{printf $1"\t";for(i=2;i<=NF;i=i+2)printf $i"\t";print $i}'
网友评论