运行代码如下
geneTerm=data.frame()
term=read.table(termFile, header=T, sep="\t", check.names=F)
for(molude in unique(term[,"Cluster"])){
term1=term[term[,"Cluster"]==molude, , drop=F]
molGenes=term1[1, "geneID"]
molGenes=unlist(strsplit(molGenes,"/"))
molTerm=term1[1,"Description"]
geneTerm=rbind(geneTerm, cbind(molGenes, molTerm))
}
row.names(geneTerm)=geneTerm[,1]
sameGene=intersect(row.names(geneTerm), row.names(data))
geneTerm=geneTerm[sameGene,2,drop=F]
colnames(geneTerm)=c("Term")
data=data[sameGene,]
结果报错
Error in strsplit(molGenes, "/") : 非字符参数
怎么改一下呢
你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答
本次提问扣除的有问必答次数，将会以问答VIP体验卡（1次有问必答机会、商城购买实体图书享受95折优惠）的形式为您补发到账户。
因为有问必答VIP体验卡有效期仅有1天，您在需要使用的时候【私信】联系我，我会为您补发。