Мне отчаянно нужна помощь. У меня есть переменная непрерывного отклика, фиксированный эффект лечения (переменная фактора с тремя уровнями), фиксированный эффект времени (переменная фактора с девятью фиксированными моментами времени) и эффект субъекта, который является случайным эффектом (у меня 34 субъекта, в основном каждый субъект пробует каждую обработку в каждый момент времени). Я пытаюсь приспособить смешанную модель с неструктурированной ковариацией, используя следующую команду, и я продолжаю получать ошибку ниже. Кроме того, для того, чтобы R вернул ошибку, требуется EVER. Проблема, кажется, в аспекте corSymm. Я никогда не делал этого раньше, и R все еще довольно труден для меня, поэтому любая помощь будет принята с благодарностью. Спасибо!
cs1 <- corSymm(form = ~ 1 | Subject)
cs1 <- Initialize(cs1, data = Mydata)
model <- gls(Response ~ Time+Treatment+Treatment*Time, correlation = cs1, data = Mydata)
anova(model)
Ошибка в gls (Ответ ~ Время + Лечение + Лечение * Время,
корреляция = cs1,: предел оценки функции достигнут без
схождение (9)