Я пытаюсь вписать модель линейного роста (LGM) в R, и я понимаю, что первичными шагами будет согласование нулевой модели со временем в качестве предиктора моей независимой переменной Y (с учетом случайных эффектов) и Нулевая модель не допускает случайных эффектов, затем сравните их и посмотрите, достаточно ли случайный эффект, чтобы оправдать использование модели со случайным перехватом.
Мне удалось согласовать модель со случайным перехватом с помощью функции lmer пакета lme4, но я не могу найти в этом пакете функцию, которая позволила бы мне подобрать модель без случайного перехвата.
Я пытался согласовать модели как со случайным перехватом (функция lme), так и без (функция gls) с пакетом nlme, но ни один из них не работал для меня.
Мой оригинальный код:
library(nlme)
LMModel <- lme(Y~Time, random=~Time| ID, data=dataset,
method="ML")
и, выполнив это, я получил сообщение об ошибке «отсутствуют значения в объекте» (очевидно, ссылаясь на мою переменную Time). Таким образом, я добавил преобразование своего набора данных в матрицу с помощью «matr <- as.matrix (dataset)» и добавил недостающую часть управления данными в мой код, который в итоге стал: </p>
LMModel <- lme(Y~Time, random=~Time| ID, data=dataset,
method="ML", na.action = na.exclude(matr))
Запустив это, я получаю сообщение об ошибке: «не удалось найти функцию« 1 »»
Далее я попытался подобрать модель без случайного эффекта с помощью функции gls nlme и получил точно такую же ошибку.
Я чувствую себя совершенно растерянным, так как не могу понять, что означает эта функция 1. Есть идеи о том, что здесь может происходить?
Заранее большое спасибо за помощь!
Federico