假设每个病人的基因在网络中的度是一个csv文件,将这些csv合并成一个矩阵,程序如下
mt是自己定义的一个空矩阵,行数等于所有基因的数目,列数等于csv文件的数目
for(dif in 1:ncol(mt))
{sss<-paste("csv文件的地址",dif,sep="/")
Mutationname<-read.csv(paste(sss,"csv",sep="."))
for(i in 1:nrow(Mutationname))
{
if(Mutationname[i,2]==0)
{
next
}
else
{
a<-as.character(Mutationname[i,1])。
mt[a,dif]<-Mutationname[i,2]
}
}
print(dif)
}