【发布时间】:2014-10-19 13:22:55
【问题描述】:
我有一个大 (>100,000) 单列浮点时间序列数据。我想找到数据中关于时间的结构变化(在我的案例索引中)。为了做到这一点,我正在使用 R strucchange 包。但我遇到了错误。我希望用我的数据集生成以下类型的图
这是我正在使用的代码:
data = read.csv("tst1.csv", header = FALSE)
library(strucchange)
## fit, visualize and test OLS-based CUSUM test
## with a Cramer-von Mises functional
ocus <- efp(data ~ 1, type = "OLS-CUSUM")
sctest(ocus, functional = "meanL2")
## estimate breakpoints
bp <- breakpoints(data ~1)
summary(bp)
plot.ts(data)
lines(bp)
lines(fitted(bp, type = "mean"), col = 2)
我收到以下错误消息:
Loading required package: zoo
Attaching package: ‘zoo’
The following objects are masked from ‘package:base’:
as.Date, as.Date.numeric
Loading required package: sandwich
Error in model.frame.default(formula, data = data) :
invalid type (list) for variable 'data'
Calls: efp -> model.frame -> model.frame.default
Execution halted
我的数据 tst1.csv 看起来像:
0.6
0.78
0.54
0.96
0.43
1.43
1.03
0.83
0.68
3.23
3.09
2.87
...
非常感谢任何帮助。
谢谢
【问题讨论】:
-
您在哪里确切收到错误消息?我会假设打电话给
efp()。str(data)的输出是什么? -
@StephanKolassa,当我使用 Rscript test.R 在终端上运行代码时,我收到了此错误消息。我是 R 的新手。我不知道如何获取 str(data)。但是,当我在 R 控制台上输入数据时,我得到一个数据框。谢谢
标签: r machine-learning time-series sampling