私はR次のような非線形モデルを扱ってきました。
Y = Alpha1*time + Alpha2*sin(2*pi*time/Alpha3) + Alpha4*(-1)^time
そして、ベルヌーイ変数がAlpha1*timeand 切片に影響するかどうかを並べます。そのようなベルヌーイ変数は次のようになります。
varia<-rep("BEFORE","AFTER"),each=30)
持つ
Y<- -2.5+rnorm(60)+2*sin(2*pi*time/8)+2.5*(-1)^time
time<-seq(1,60)
がこのモデルに適合することがわかりましたが、このカテゴリ変数の効果がなければ、パッケージは非線形モデルの項ではなく変数に対する効果を推定できることnlsもわかりました。nlme私の質問は次のとおりです。これらのパッケージのどれが私を助けることができますか? そして、このベルヌーイ変数をコードに追加するにはどうすればよいでしょうか? よろしくお願いします。