Как разрешить фактор-специфическую дисперсию случайного эффекта в lme

Question

Я предполагаю, что отклонения случайных эффектов в моей модели смешанного эффекта будут разными для разных уровней фиксированного коэффициента 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.

Answer 1

Это Rpub обсуждает, как сделать конкретные факторы-отклонения в обоих LME и lmer:

http://rpubs.com/bbolker/6298