提取百分比解释了ssa结果的差异（Rssa）

Question

我正在使用Rssa软件包来分解时间序列，女巫工作得很好，除了我不能从每个特征向量得到解释方差的百分比（如果这些是解释这个的正确词）。但是，这些百分比在我可以用这个包绘制的图表的顶部注明。 让我举个例子：

d=rnorm(200,10,3)
plot(d,type="l")
ssa=ssa(d, L = 100,digits=0)
plot(ssa,type="vector")   #the percentage I want is in the title of each individual graph

# to reconstruct the trend and the residuals
res <- reconstruct(ssa, groups = list(1))
trend <- res$F1

如何在矢量中获得这些百分比？特别是因为我想循环多个系列。 谢谢！

Answer 1

似乎组件系列的加权范数代码隐藏在包中。 我从Rssa:::.plot.ssa.vectors.1d.ssa中提取代码并将其包装成一个小函数：

component_wnorm <- 
function(x) {
  idx <- seq_len(min(nsigma(x), 10))
  x <- ssa
  total <- wnorm(x)^2
  round(100*x$sigma[idx]^2 / total, digits = 2)
}

component_wnorm(ssa)
[1] 92.02  0.35  0.34  0.27  0.27  0.25  0.22  0.20  0.20  0.18

Answer 2

最新版本的Rssa具有功能contributions。 因此，您可以使用

> s <- ssa(d, L=100)
> c <- contributions(s)*100
> print(c[1:10], digits = 2)
 [1] 92.41  0.28  0.26  0.26  0.26  0.23  0.23  0.21  0.20  0.20