可以学习我们的课程,高级分析内容里面有详细微生物共存网络分析代码讲解和示例数据;https://bdtcd.xetslk.com/s/qZRiF
我在R中运行如下代码:
Abu=read.table('family_nor.out',header=T)
Abu<-as.matrix(Abu)
#Min support filtering for each taxa
Abu[Abu<0.01]<-0
#Filter taxa by abundance sum
Abu<-Abu[which(rowSums(Abu)>10),]
#Filter taxa by Occupancy
table<-Abu
table[table>0]<-1
table1<-Abu[which(rowSums(table)>=100),]
Abu<-table1
data<-co_occurrence_network(Abu,0.5,0.001)
在这最后一步总是出现报错: Error in rcorr(t(matrix), type = "spearman") : must have >1 column
别人用相同的代码都可以正常运行,但是我却不行,我猜想可能是数据格式的问题,但是修改了很多遍依然没能解决这个问题,希望老师可以不吝赐教,谢谢
可以学习我们的课程,高级分析内容里面有详细微生物共存网络分析代码讲解和示例数据;https://bdtcd.xetslk.com/s/qZRiF
如果觉得我的回答对您有用,请随意打赏。你的支持将鼓励我继续创作!