r - 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。