1
Rの非線形混合効果回帰
驚いたことに、Googleを使用して次の質問に対する答えを見つけることができませんでした。 私はいくつかの個人からのいくつかの生物学的データを持っていますが、それはおおよそS字状の成長挙動を示しています。したがって、標準のロジスティック成長を使用してモデル化したい P(t) = k*p0*exp(r*t) / (k+p0*(exp(r*t)-1)) p0はt = 0での開始値、kはt-> infinityでの漸近極限、rは成長速度です。私が見る限り、nlsを使用してこれを簡単にモデル化することができます(私の理解が不足しています:時間とデータをスケーリングすることにより、標準のロジット回帰を使用して類似のものをモデル化できないのはなぜですか?編集:ニック、どうやら、例えばプロポーションですが、めったにhttp://www.stata-journal.com/article.html?article=st0147この接線に関する次の質問は、モデルが異常値を処理できるかどうかです> 1)。 ここで、3つのパラメーターk、p0、およびrに対する固定(主にカテゴリー)およびランダム(個別のID、場合によってはスタディID)の効果を許可したいと思います。nlmeはこれを行う最良の方法ですか?SSlogisモデルは、私がやろうとしていることに対して賢明なようです、それは正しいですか?次のいずれかが賢明なモデルから始まりますか?開始値を正しく取得できないようで、update()はランダムな効果に対してのみ機能し、固定された効果に対しては機能しないようです-ヒントはありますか? nlme(y ~ k*p0*exp(r*t) / (k+p0*(exp(r*t)-1)), ## not working at all (bad numerical properties?) data = data, fixed = k + p0 + r ~ var1 + var2, random = k + p0 + r ~ 1|UID, start …