我有 2 两个相异矩阵。一个是在 111 个站点之间比较观察到的数据，另一个是使用空模型生成的。
我想使用adnois函数vegan来测试观察到的差异是否与空模型预期的差异显着。然而，阿多尼斯函数只会在公式的左侧采用一个相异矩阵。
有谁知道如何为这个测试建模？
谢谢
问问题
1409 次
1 回答
2
这个问题的答案是：
meanjac <- function(x) mean(vegdist(x, method='jaccard', diag=TRUE))
test <- oecosimu(x, nestfun=meanjac, method="r1", nsimul = 10^3, statistic='adonis')
它将一个函数传递给 oecosimu，以获取 Jaccard 相异矩阵的平均值，然后使用 'r1' 方法通过随机打乱二元社区矩阵但根据观察到的占有率分配物种占有率并比较它来生成空社区矩阵到观察到的相异矩阵。
感谢 Jari 为我指明了正确的方向......
于 2014-03-11T16:39:27.523 回答