0

上手、

データに対して区分的回帰を実行して、各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

さて、誰かがそれをどのように扱うか教えてくれれば

前もって感謝します。

英語が下手でごめんなさい。

4

0 に答える 0