Я пытаюсь подогнать под смешанную модель с помощью функции lmer
из пакета lme4
. Однако я не понимаю, что следует вводить в параметр start
. Моя цель - использовать простую линейную регрессию, чтобы использовать оцененные там коэффициенты в качестве начальных значений для смешанной модели.
Допустим, моя модель следующая:
linear_model = lm(y ~ x1 + x2 + x3, data = data)
coef = summary(linear_model)$coefficients[- 1, 1] #I remove the intercept
result = lmer(y ~ x1 + x2 + x3 | x1 + x2 + x3, data = data, start = coef)
Этот пример представляет собой упрощенную версию того, что я делаю, поскольку я не смогу поделиться своими данными.
Тогда я получаю ошибку следующего вида:
Error during wrapup: incorrect number of theta components (!=105) #105 is the value I get from the real regression I am trying to fit.
Я пробовал много разных решений, пытаясь предоставить список и назвать эти значения theta
, как я видел на некоторых форумах.
Также код Github проверяет, подходит ли длина, но я не могу найти, к чему он относится:
# Assign the start value to theta
if (is.numeric(start)) {
theta <- start
}
# Check the length of theta
length(theta)!=length(pred$theta)
Однако я не могу найти, где определено pred$theta
, и поэтому не понимаю, откуда взялось это значение 105
.
Любая помощь ?