複数の異なるランダム効果をnlmeで指定する

Question

私はいくつかの鯨観光データを分析し、nlmeパッケージに線形混合効果モデルを構築しようとしています。 （私もlme4にこのモデルを実行するに開いています。）

私の変数は以下のとおりです。

• mins：出会いの時間（応答変数）
• Id：個々のクジラID（ランダム効果）
• Vessel：血管Id（ランダム効果）
• Sex：動物の性別
• Length：動物の長さ（Year内にネスト）
• Year
• Month。

だから私の確率変数はIdとVesselであり、私はまた、ネストされたランダム効果としてYearとMonthを持っています。

私は次のように出ている：

form1 <- formula(Min ~ length + Sex+ Encounter) 
model1 <- lme(form1, 
       random = list(Id = ~1, 
          Vessel = ~1, 
          Year=~1, 
          Month = ~1), data=wsdata, method="ML") 

しかし、すべての私のランダム効果はId内にネストになります。

私はネストされたランダム効果としてIdとVessel別々のランダム効果などとYearとMonthを定義することができる方法はありますか？一般的に

Answer 1

は、それはあなたがまだ達成しやすくなり、時間的または空間的自己相関や不均一のためのモデルを（必要ない限り、lme4には（あなたは「別の」によって何を意味するかは、私は思う）交差ランダム効果を指定する方がはるかに簡単ですencounter何であるか

• ：）nlmeで、私はいくつかの他のコメント

  library(lme4) 
  fit <- lmer(mins ~ Length + Sex+ (1|Id) + (1|Vessel) + 
         (1|Year/Month), data=wsdata, REML=FALSE) 
  

  で先に行くでしょうか？それはあなたの公式ではあったが、データセットのあなたの説明にはなかった。

• 遭遇時間（遭遇時間？）は歪んでいる可能性が高い。その場合、それらをログ変換したいかもしれない。