5

我假设我的混合效应模型中的随机效应方差对于不同水平的固定因子会有所不同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，但是这个函数不允许考虑组内误差的方差异质性（在我的情况下最好考虑）。我的问题是如何使用函数为随机效应合并“Btyp”特定的方差lme？

下面你可以看到它是如何与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   21951

Random effects:
 Groups   Name Variance Std.Dev. Corr 

 Probe    Pind 158.6058 12.5939         
          pF     2.4289  1.5585  -1.000 

 Probe    Bind 134.6383 11.6034         
          pF     2.7619  1.6619  -1.000 

 Probe    Tind 490.6714 22.1511         
          pF    46.3533  6.8083  -1.000 

 Residual      316.9860 17.8041    

Number of obs: 2530, groups: Probe, 45

Pind, Bind,Tind是不同水平的指标变量BTyp。

4

0 回答 0