#r #glm #linearmodels
Вопрос:
Я пытался подогнать линейную смешанную модель к набору данных со следующей специальной характеристикой:
Существует одна категориальная переменная A, принимающая значения {1,2,…, K}, которая считается фиксированным эффектом для {1,2} и случайным эффектом для {3,…,K}. Как я должен реализовать это в функции lmer внутри R?
Чтобы быть более конкретным, рассмотрим простейшую модель y_i = b_0 A ошибка.
(Близким решением будет оценка двух моделей отдельно для {1,2} и {3,…,K}, но это влияет на оценку дисперсии условий ошибки.)