小编Olg*_*kis的帖子

我假设混合效应模型中的随机效应方差对于固定因子的不同水平将有所不同BTyp。

这是我的模特

fm2 <- lme(CA ~ 1 + pF+Tiefe+BTyp+Tiefe:pF+BTyp:pF, data=data2, 
           random = list(~ 1 + pF|Probe))
fm2_Btyphet<-update(fm2, weights=varIdent(form=~1|BTyp))

我Btyp使用lmer函数设法将特定的方差纳入随机效应中，但是此函数不允许考虑组内误差的方差异质性（在我的情况下最好考虑）。我的问题是如何使用lme函数将特定于“ Btyp”的差异纳入随机效应？

在下面您可以看到它如何与lmer功能一起使用。

CA ~ 1 + pF + Tiefe + BTyp + Tiefe:pF + BTyp:pF + 
     (0 + Pind + pF | Probe) + (0 + Bind + pF | Probe) + (0 + Tind + pF | Probe) 


 Data: data2 

 AIC   BIC logLik deviance REMLdev

   21987 22092 -10975    21979 …