我有一个基于所有真菌丰度的混沌距离矩阵:
CR10 CR11 CR13 CR14 CR17 CR18 CR19
CR11 0.4531840
CR13 0.4288178 0.4624915
CR14 0.5903908 0.5466617 0.4942469
CR17 0.4784990 0.3387325 0.6136265 0.5779121
CR18 0.7649840 0.7537409 0.7526077 0.5632825 0.4153391
CR19 0.3772907 0.4579895 0.3208187 0.3706775 0.5644193 0.7380274
CR20 0.4598706 0.5529427 0.6424340 0.6690386 0.3855154 0.5509150 0.6406800
和具有相同环境图的33个环境变量的表。
我跑步时:
fungAbundAdonis <- lapply(colnames(env2), function(x) {
form <- as.formula(paste("OTU.table2", x, sep="~"))
z <- adonis(form, data = env2, permutations=999)
return(data.frame(env = rownames(z$aov.tab), Rsq = z$aov.tab$R2,P = z$aov.tab$P))}
)
我收到此错误:
adonis错误(form,data = env2,排列= 999): 公式的右侧没有可用的术语。
我不明白为什么,因为当我将相同的脚本与图的距离矩阵从1到9以及12 15和16以及这些图的环境表一起使用时,它工作正常。有人知道错误的根源是什么吗? 谢谢!
答案 0 :(得分:1)
您的问题没有可重复的例子,我不得不猜测。但是,如果变量在右侧是 constant ,我可以重现您的错误消息。当您子集env2
且在该选定子集中一个变量只有一个值时可能会发生这种情况。 (这仅涉及vegan 2.5-x或发行版本:vegan 2.6-0不会显示错误消息。)