Предположим, у нас есть следующая линейная модель смешанных эффектов: 
Как мы можем вписать эту вложенную модель в R?
Сейчас я попробовал две вещи:
- Перепишите модель следующим образом:
Затем с помощью функции lmer в пакете lme4 поместите модель смешанного эффекта и укажите Xi как случайный, так и фиксированный эффект, ковариантный как: lmer(y ~ X-1+(0+X|subject))
Но когда я передаю результат BIC и делать выбор модели, он всегда выбирает простейшую модель, что не правильно.
Сначала я попытался регрессировать
y_i на
X_i и обработать
X_i как фиксированный эффект, затем я получу оценку наклона, то есть вектора
phi_i.Затем посмотрите
phi_i как новые наблюдения и снова наберите
C_i, чтобы получить
beta.Но это кажется неправильным, поскольку мы не знаем
C_i в реальной задаче, и похоже, что
C_i и
beta совместно определяют коэффициенты.
Итак, есть ли другие способы подгонкиэтот тип модели в R и где мои ошибки?
Спасибо за любую помощь!