代碼如下：setwd("目錄")
expFile="rocSigExp.txt"
clinicalFile="clinical.txt"
exp=read.table(expFile,sep="\t",header=T,check.names=F,row.names=1)
cli=read.table(clinicalFile,sep="\t",header=T,check.names=F,row.names=1)
exp=exp[,3:ncol(exp)]
samSample=intersect(row.names(exp),row.names(cli))
exp=exp[samSample,]
cli=cli[samSample,]
pFilter=0.05
outTab=c()
outTab=rbind(outTab,c("id",colnames(cli),"SigNum"))
colnames(outTab)=c("id",colnames(cli),"SigNum")
for(i in colnames(exp)){
clinicalPvalVector=c()
sigSum=0
for(clinical in colnames(cli)){
tab1=table(cli[,clinical])
labelNum=length(tab1)
rt1=cbind(expression=exp[,i],clinical=cli[,clinical])
if(labelNum==2){
cliTest<-wilcox.test(expression ~ clinical, data=https://bbs.csdn.net/topics/rt1)
}else{
cliTest<-kruskal.test(expression ~ clinical, data = rt1)}
pValue=https://bbs.csdn.net/topics/cliTest$p.value
clinicalPvalVector=c(clinicalPvalVector,pValue)
if(pValue<pFilter){
sigSum=sigSum+1
}
}
geneClinical=c(i,clinicalPvalVector,sigSum)
outTab=rbind(outTab,geneClinical)
}
write.table(outTab,file="clinicalCor.xls",sep="\t",col.names=F,row.names=F,quote=F)
運行的時候出現報錯
uj5u.com熱心網友回復：
有沒有大神來解答一下呀轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/284155.html