Это дополнительный вопрос к Какая модель для моих данных? (тестирование различий наклона для трех групп).
Решение сработало для примерного набора данных (см.связанный вопрос), но я столкнулся с ошибкой при использовании другого набора данных:
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
Теперь я хотел бы знать, почему эта ошибка происходит.Поскольку решение работает для другого набора данных, оно может иметь какое-то отношение к вышеуказанным данным?