上手、
データに対して区分的回帰を実行して、各IDの最終的なブレークポイントを見つけ、この情報を抽出しようとしています。
例:私の個人データのサブセット
ID time y
7G009 0 9
7G009 108,33 13
7G009 185,69 16
7G009 309,22 20
7G009 515,08 21
7G051 0 10
7G051 108,33 14
7G051 185,69 19
7G051 309,22 23
7G051 515,08 25
8S027 0 8
8S027 108,33 13
8S027 185,69 17
8S027 309,22 22
8S027 515,08 23
その目的のために、私はパッケージstrucchangeとdlplyのブレークポイント関数を使用しました。私のコードは次のとおりです。
breakby=dlply(test,.(subject),summarize,break=breakpoints(y~time,h=2))
しかし、エラーが発生します:
Erreur dans breakpoints.formula(nbf ~ cumggd, h = 2) :
minimum segment size must be greater than the number of regressors
さて、誰かがそれをどのように扱うか教えてくれれば
前もって感謝します。
英語が下手でごめんなさい。