如何使用“ how”功能进行不平衡的重复设计

Question

我有一组对照样地和经过处理的样地，这些样地是多年来进行采样的。我在vegan包中运行了prc函数，并希望执行置换测试以检查多年来控制图和处理图的差异是否很大。由于我的数据不平衡，因此我无法使用函数。我的代码如下：

library(vegan)
year=as.factor(c(rep(1995,8),rep(1999,8),rep(2001,8),rep(2013,4),rep(1995,4),
       rep(1999,4),rep(2001,4),rep(2013,4)))
treatment=as.factor(c(rep("control",28),rep("treated",16)))

我已经写了这个，但是我确定这是错误的，因为这里缺少治疗方法：

h1 <- how(within = Within(type = "series", mirror = F),
          blocks = year, nperm = 999
)

任何建议，我们将不胜感激。

我有一组对照样地和经过处理的样地，这些样地是多年来进行采样的。我在纯素食包中运行prc函数，并希望执行置换测试以检查是否对控制进行处理...

Answer 1

在原假设下，control或treated组的样本

是]可交换的，因此您不希望在置换设计中使用它们；您真的想对它们进行置换以为测试统计信息生成基于置换的空分布。