Я создаю линейную смешанную модель, для которой я заранее оценил / определил компоненты дисперсии.Теперь я хотел бы добавить их в функцию, но я не знаю, как это сделать.Используемая в настоящее время функция «lmer» из пакета «lme4» действительно сама оценивает компоненты дисперсии, и я не знаю, как их предварительно определить.
Заранее спасибо!
For Example this data:
library(lme4)
f <- rnorm(10, mean=1.4025, sd=0.101)
m <- rnorm(10, mean=1.9025, sd=0.15)
gain <- c(f,m,f,m)
s <- c("f","f","f","f","f","m","m","m","m","m")
sex <- c(s,s,s,s)
i <- c(1:10)
anim <- as.factor(c(i,i,i,i))
invented <- data.frame(gain,sex,anim)
# model
lmm <- lmer(gain ~ sex + (1|anim), data = d)
summary(lmm)
# How can I add variance components to the model?
# Fixed-effect-variance = 0.0135
# Random-effect-variance = 0.0045