ggwang 发表于 2012-6-4 18:07 
崇拜楼上的,不过以楼上编程的能力应当会发出更好的文章的,我一篇研究要用到SupF、MeanF以及Lc这几个指标, ...
你是在进行时间序列的断点检测吗?
R-project的 AER,strucchange PARCKAGE可以用来检验时间序列的稳定性。
http://ftp.ctex.org/mirrors/CRAN ... nge/strucchange.pdf
http://cran.r-project.org/web/packages/AER/AER.pdf
以下代码来源于AER.PDF,调用了strucchange的supF函数
library("strucchange")
infl <- usm[, "infl"]
unemp <- usm[, "unemp"]
usm <- ts.intersect(diff(infl), lag(diff(infl), k = -1), lag(diff(infl), k = -2),
lag(diff(infl), k = -3), lag(diff(infl), k = -4), lag(unemp, k = -1),
lag(unemp, k = -2), lag(unemp, k = -3), lag(unemp, k = -4))
colnames(usm) <- c("dinfl", paste("dinfl", 1:4, sep = ""), paste("unemp", 1:4, sep = ""))
usm <- window(usm, start = c(1962, 1), end = c(2004, 4))
fs <- Fstats(dinfl ~ ., data = usm)
sctest(fs, type = "supF")
plot(fs)