这是Which model for my data?的后续问题(测试三组斜率的差异)。
该解决方案适用于示例性数据集(请参阅链接的问题),但是使用另一个数据集时遇到错误:
Input = ("
Group Time Size
A 1 1.08152
A 2 1.10589
A 3 1.13292
B 1 1.04597
B 2 1.05763
B 3 1.07023
B 4 1.08612
B 5 1.10059
B 6 1.11589
B 7 1.13143
B 8 1.14741
B 9 1.16721
B 10 1.18288
C 1 1.04777
C 2 1.06145
C 3 1.07484
C 4 1.08908
C 5 1.10346
C 6 1.11866
C 7 1.13375
C 8 1.14931
C 9 1.16563
C 10 1.18294
")
dat = read.table(textConnection(Input),header=TRUE)
建立模型时
(m1 <- gls(Size ~ Time * Group, dat, correlation = corAR1(form = ~ Time | Group), weights = varIdent(form = ~ 1 | I(Group == "A"))))
并进行成对测试
pairs(emtrends(m1, ~ Group, var = "Time"))
我遇到此错误:
Error in crossprod(x, y) : requires numeric/complex matrix/vector arguments
现在我想知道为什么会发生此错误。由于该解决方案适用于其他数据集,它可能与上述数据有关?