Я использую смешанную модель в R
( lme4
) для анализа некоторых данных повторных измерений. У меня есть переменная реакции (содержание волокна в кале) и 3 фиксированных эффекта (масса тела и т. Д.). В моем исследовании всего 6 участников, по 16 повторных измерений для каждого (хотя у двух только 12 повторений). Субъектами являются ящерицы, которым давали разные комбинации пищи в разных «обработках».
Мой вопрос: могу ли я использовать идентификатор субъекта в качестве случайного эффекта?
Я знаю, что это обычный курс действий в моделях продольных смешанных эффектов, чтобы принять во внимание случайную выборку природы субъектов и тот факт, что наблюдения внутри субъектов будут более тесно коррелировать, чем наблюдения между субъектами. Но трактовка идентификатора субъекта как случайного эффекта предполагает оценку среднего значения и дисперсии для этой переменной.
Поскольку у меня есть только 6 предметов (6 уровней этого фактора), достаточно ли этого, чтобы получить точную характеристику среднего значения и дисперсии?
Помогает ли в этом отношении тот факт, что у меня есть довольно много повторных измерений для каждого субъекта (я не понимаю, как это имеет значение)?
Наконец, если я не могу использовать идентификатор субъекта в качестве случайного эффекта, позволит ли его включение в качестве фиксированного эффекта контролировать тот факт, что я повторил измерения?
Изменить: я просто хотел бы уточнить, что когда я говорю «могу ли я» использовать идентификатор субъекта в качестве случайного эффекта, я имею в виду «это хорошая идея для». Я знаю, что могу подобрать модель с коэффициентом всего 2 уровня, но наверняка это будет неоправданно? Я спрашиваю, в какой момент становится разумным думать о том, чтобы рассматривать предметы как случайные эффекты? Кажется, литература советует, что 5-6 уровней - это нижняя граница. Мне кажется, что оценки среднего значения и дисперсии случайного эффекта не будут очень точными, пока не будет более 15 уровней факторов.