如何使用dlply将package strucchange的断点函数应用于数据子集?[R]

时间:2012-07-31 08:59:32

标签: r plyr piecewise

好吧,

我尝试对我的数据进行分段回归,以找到每个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

为此,我使用了package 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

好吧,如果有人能告诉我如何处理它

提前谢谢。

抱歉英语不好。

0 个答案:

没有答案