Я пытаюсь проанализировать свои данные с помощью моделей роста . Люди были разделены на две группы: одна группа прошла обучение, а другая - нет. Моя зависимая переменная - MAS, мои данные - MAS_long, и у меня есть три момента измерения. Я хочу посмотреть только на мою группу вмешательства и на то, как каждый человек изменился со временем.
Все работало хорошо, когда я просматривал все данные, не отделяя свою группу вмешательства. Затем я построил синтаксис, основываясь только на данных моей группы вмешательства. Я центрировал свою переменную MAS, и код больше не работает . Это работает хорошо, когда я не центрирую MAS, и это работает с центрированной зависимой переменной MAT. Но MAS не работает ...
Вот как я построил свою модель роста (на основе Энди Филда - 2012):
Intercept_MAS <- gls(MAS ~ 1, data=MAS_long, method = "ML", na.action = na.exclude)
randomIntercept_MAS <-lme(MAS ~ 1, data = MAS_long, random = ~1|Nr, method = "ML", na.action = na.exclude)
time_randomIntercept_MAS<- update(randomIntercept_MAS, .~. + time) #I added time as fixed effect
summary(time_randomIntercept_MAS)
, пока здесь все не идет хорошо. Затем я добавляю случайный наклон:
time_randomSlope_MAS <- update(time_randomIntercept_MAS, random = ~time|Nr, control = list(opt = 'optim'))
И появляется эта ошибка:
Ошибка в lme.formula (исправлено = MAS ~ время, данные = MAS_long,: проблема Optim,код ошибки конвергенции = 1
message =
К сожалению, после «message =» нет сообщения. У кого-нибудь есть идея, что R хочет мне сказать?