R の非線形混合モデルを自分のデータに当てはめようとしていますが、人々が使用しているヘルプ ファイルと例を理解するのに苦労しています。私は修士課程の学生で、22 日間にわたって 4 頭の羊の福祉対策として行動の変化をモデル化しようとしています。-4 日目から -1 日目はベースラインとして使用され、0 日目は手術中のため測定値がなく、1 日目から 17 日目は私が興味を持っているものです。さまざまなデータセットをプロットすると、手術後 1 日目から 3 日目にかけてかなり落ち込んでいます。多かれ少なかれベースラインに戻ります。最初に説明変数 (Day) を 2 乗するように言われましたが、lmer() を使用しても結果はまったく変わりません。だから私はnlmerを試して、それがうまくいくかどうかを見たい.
私のデータセットは次のようになり、500 行ほどになります。「羊」は因子として読み取られます。
. Behaviour Sheep Day DaySq Block Observed Frequency
1 Standing 2 -4 16 0 49 71.01449
2 Lying 2 -4 16 0 12 17.39130
3 Eating 2 -4 16 0 36 52.17391
4 Ruminating 2 -4 16 0 16 23.18841
5 Moving in pen 2 -4 16 0 0 0.00000
6 XDisturbance 2 -4 16 0 9 13.04348
私のlmerモデルは次のようになります
lmer(Frequency ~ Behaviour * Day + (1|Sheep), behdata2)
私は似たようなことをしたいのですが、非線形として。それは可能ですか？いくつかのバージョンを試しましたが、これを正しく行う方法がわからず、エラー メッセージが表示され続けます。また、R ヘルプ ファイルの例は、私が持っていない多くの知識を前提としているようです。どんな助けでも大歓迎です！