R strucchange RSS and BIC for one breakpoint

为君一笑 提交于 2019-12-12 21:23:34

问题


The reference manual for the R strucchange package states:

As the maximum of a sequence of F statistics is equivalent to the minimum OLS estimator of the breakpoint in a 2-segment partition it can be extracted by breakpoints from an object of class "Fstats" as computed by Fstats.

Is the breakpoint extracted by the above approach always the same as the breakpoint extracted by applying coeftest to an object of class "breakpointsfull" when the "breaks" option is 1? In other words, is the breakpoint extracted by calling breakpoints on an Fstats object (minimum RSS) always equal to the breakpoint extracted by calling coeftest on a breakpointsfull object (minimum BIC) when breaks = 1?

library(lmtest)
library(strucchange)

data("Nile")

fs.nile <- Fstats(Nile ~ 1)
breakpoints(fs.nile)

bp.nile <- breakpoints(Nile ~ 1)
coeftest(bp.nile, breaks = 1)

回答1:


Essentially yes. The breakpoint associated with the supF statistic (supWald or supLR) and the minimum RSS breakpoint in the linear regression model are identical.

Of course, the trimming/minimal segment size needs to be the same. The default is 15% in both breakpoints() (argument h) and Fstats() (arguments from and to). Also, no alternative covariance matrix estimator (vcov) must be used in Fstats().

However, breakpoints(Fstats(formula, ...)) just extracts a simple object of class "breakpoints" while breakpoints(formula, ...) sets up the whole "breakpointsfull" object from which confidence intervals, coefficient tests, information criteria etc. can be computed.



来源:https://stackoverflow.com/questions/30562109/r-strucchange-rss-and-bic-for-one-breakpoint

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!