我正在学习如何做PERMANOVA
,当我从这个例子开始时,这真的是一个挑战。以下代码来自adonis
中的帮助页面。我可以通过这个例子,但无法弄清楚为什么它是incorrect (no strata)
?这是否意味着我必须获取分层数据才能使用adonis
?
data(dune)
data(dune.env)
### Example of use with strata, for nested (e.g., block) designs.
dat <- expand.grid(rep=gl(2,1), NO3=factor(c(0,10)),field=gl(3,1) )
Agropyron <- with(dat, as.numeric(field) + as.numeric(NO3)+2) +rnorm(12)/2
Schizachyrium <- with(dat, as.numeric(field)-as.numeric(NO3)+2) +rnorm(12)/2
Y <- data.frame(Agropyron, Schizachyrium)
mod <- metaMDS(Y)
plot(mod)
### Hulls show treatment
with(dat, ordihull(mod, group=NO3, show="0"))
with(dat, ordihull(mod, group=NO3, show="10", col=3))
### Spider shows fields
with(dat, ordispider(mod, group=field, lty=3, col="red"))
### Correct hypothesis test (with strata)
adonis(Y ~ NO3, data=dat, strata=dat$field, perm=999)
### Incorrect (no strata)
adonis(Y ~ NO3, data=dat, perm=999)
答案 0 :(得分:1)
在PERMANOVA
中,PER
位代表排列测试。排列测试仅在您根据研究的实验设计洗涤真正可交换的样本时才有效。在这个例子中,该示例使用了一个被阻止的研究,因此如果您在的strata
级别内置换样本,并且从不在级别之间置换,则排列测试仅有效strata
。
评论Incorrect (no strata)
所指的是什么; 这些数据的无限制排列测试不正确。