Question

我有一个关于时间序列填空的问题，使用R包{Rssa}。示例系列的长度为334，包含几个NA。

> my.L = length(b1)%/%2
> s1.gap <- ssa(b1,L=my.L)
Error in ssa(b1, L = my.L) : 
  Nothing to decompose: the given field shape is empty
In addition: Warning messages:
1: In eval(expr, envir, enclos) :
  Some field elements were not covered by shaped window. 292 elements will be ommited
2: In eval(expr, envir, enclos) :
  Nothing to decompose: the given field shape is empty

我想用igapfill {Rssa}填补我的NAs。由于igapfill想要一个SSA形状的物体，我必须先做一个ssa，然后出现问题：

> s2.gap <- ssa(b1,L=23)
Warning message:
In eval(expr, envir, enclos) :
  Some field elements were not covered by shaped window. 269 elements will be ommited

> summary(s2.gap) #for L=23

Call:
ssa(x = b1, L = 23)

Series length: 334, Window length: 23,  SVD method: eigen
Special triples:  0

Computed:
Eigenvalues: 1, Eigenvectors: 1,    Factor vectors: 0

Precached: 0 elementary series (0 MiB)

Overall memory consumption (estimate): 0.008934 MiB

帮助建议（并且默认情况下）窗口长度L大约是数据长度的一半，这对我来说不起作用。只有在用L到23之后我才会收到此错误消息，但警告仍然是：

> summary(s2.gap)  #for L=15

Call:
ssa(x = b1, L = 15)

Series length: 334, Window length: 15,  SVD method: eigen
Special triples:  0

Computed:
Eigenvalues: 15,    Eigenvectors: 15,   Factor vectors: 0

Precached: 0 elementary series (0 MiB)

Overall memory consumption (estimate): 0.01062 MiB

但是我不会得到超过1个特征向量，直到我进一步减少（仍然是错误消息）。

MaterialList

我不太明白我的问题是什么，谷歌等也没有帮助我。我感谢任何建议： - ）

Answer 1

（可能我的回答为时已晚。希望它对您有所帮助。）

当然，Rssa与NA合作。但是，在ssa之前调用igapfill在Rssa中确实是不清楚的。最近一本书Singular Spectrum Analysis with R的引用： “请注意，迭代方法通常使用参数值 force.decompose = FALSE的初步调用中的ssa为避免分解，由于间隙位置，这可能是不可能的。“

也就是说，在调用s1.gap <- ssa(b1, L=my.L, force.decompose = FALSE)之前使用igapfill。所引用的书中的companion web-site可能对其他技巧很有帮助。

ssa（）时间序列gapfilling R

1 个答案: