r - 为什么在 R 中分解时间序列的 stl 函数不起作用？

Question

我收集了这些数据（80 行），并试图用它stl()来分解这个时间序列：

    V1
1  6069
2  5989
3  5882
4  5875
5  6263
6  6168
7  6209
8  5171
9  5938
10 6296
11 6314
12 6003
13 5923
14 4389
15 4985
16 6273
17 5438
18 5126
19 6571
20 6783
21 7093
22 5010
23 7683
24 7233
25 6758
26 7356
27 7437
28 7520
29 7444
30 7466
31 7673
32 7945
33 8782
34 8245
35 8369
36 7750
37 7649
38 7187
39 6949
40 7520
41 6912
42 7983
43 7406
44 8361
45 6500
46 8224
47 8985
48 7663
49 8785
50 6251
51 8453
52 8515
53 7698
54 8833
55 8991
56 8583
57 8865
58 8689
59 6276
60 8141
61 8060
62 8966
63 8233
64 6526
65 7561
66 8160
67 7622
68 7976
69 7619
70 7810
71 7674
72 6784
73 4910
74 6999
75 6665
76 6604
77 5767
78 6691
79 6163
80 2258

为了分解时间序列，我尝试使用stl()应该与对象形式类一起使用的ts：

myts <- read.table(file="clipboard", header=FALSE)
myts <- ts(myts, frequency=4)
stl(myts, "periodic")

但无论我做什么，我都会收到此错误（例如：frequency=1在调用ts函数时设置或仅ts使用比我的数据没有更多参数的情况：

Error in stl(myts, "periodic") : only univariate series are allowed

我不知道我做错了什么，因为我在网上找到的所有示例都具有相同的结构。我将不胜感激任何帮助。

Answer 1

ts假设您的数据是矩阵是一个问题。您需要指定要使用的变量。这应该有效：

x <-
  "V1
1  6069
2  5989
3  5882
4  5875
5  6263
6  6168
7  6209
8  5171
9  5938
10 6296
11 6314
12 6003
13 5923
14 4389
15 4985"
y <- read.table(text = x, header = TRUE)
myts <- ts(y, frequency = 4)
stl(myts[, 1], "periodic")

问题在这里讨论。