代码如下:
gene_14=read.csv('gene_14.csv',header=F)
gene_results_p=matrix(rep(666,1428),ncol=14)
gene_results_OR=matrix(rep(666,1428),ncol=14)
gene_names=colnames(gene_14)
for (i in c(1:length(gene_14[1,]))){
+ for (n in c(i:length(gene_14[1,]))){
+ a=fisher.test(gene_14[,i],gene_14[,n])
+ gene_results_p[i,n]=p.adjust(a$p.value,method='BH');gene_results_OR[i,n]=a$estimate
+ }
+ }
Error in gene_results_OR[i, n] <- a$estimate :
被替换的项目不是替换值长度的倍数
gene_con=list("OR VALUE",gene_results_OR,"P VALUE",gene_results_p)
write.csv(gene_con,"gene_con.csv")
按上述代码运行始终报错,请大神帮忙看看代码哪里错了
gene_14.csv是14个基因(14列),218例患者(218行)的表格,内容为0(未突变)或1(突变)