我正在尝试在 lmer 模型中执行苏格兰结测试,我想要 Prov 均值排名! 我应该在其中放入什么=??
lmer1 <- with(inDB, lmer(stVOL8 ~ Rep + (1|Prov), data=inDB))
sk2 <- SK(lmer1, which='Prov')
summary(sk2)
Error in XXlist[[ii]] <- XX : attempt to select less than one element in integerOneIndex
运行代码对代表进行排名有效!
lmer1 <- with(inDB, lmer(stVOL8 ~ Rep + (1|Prov), data=inDB))
sk2 <- SK(lmer1, which='Rep')
summary(sk2)
Goups of means at sig.level = 0.05
Means G1
9 99.91 a
5 99.88 a
8 99.86 a
4 99.85 a
2 99.84 a
1 99.82 a
3 99.72 a
7 99.61 a
6 99.42 a
我猜你正在使用
SK()
来自 ScottKnott 包。
您的问题是您正在尝试对模型的“随机效应”组件运行测试。混合模型框架的基本限制之一是不可能对随机效应水平之间的差异进行推断;您可以测试有关效应分布的假设(例如,测试 RE 的方差为零的零假设),但无法测试有关分组变量水平之间效应差异的假设。 要获得实现目标的更多帮助,您可能需要在
交叉验证