在 R 中绘制 RSSA 包中的 50 多个组件

Question

require(Rssa)
t=ssa(co2,200) #200 here should be the number of components
plot(t)       # this way it plots only the first 50 not 200!

以上代码仅生成前 50 个组件的图形。我需要绘制超过 50 个组件

我试过了

plot(t$sigma[1:200],type='l',log='y')

但是没用！

示例：类似于此案例

Answer 1

查看 ?ssa 的帮助页面，我们看到一个名为 neig 的参数，记录为；

integer, number of desired eigentriples. If 'NULL', then sane default value will be used, see 'Details'

将其用作命名参数：

 t=ssa(co2, neig=200) 
 plot(t)

并且：

> t$sigma
  [1] 78886.190749   329.031810   327.198387   184.659743    88.695271    88.191805    52.380502
  [8]    40.527875    31.329930    29.409384    27.157698    22.334446    17.237926    14.175096
 [15]    14.111402    12.976716    12.943775    12.216524    11.830642    11.614243    11.226010
 [22]    10.457529    10.435998  snipped the remaining 200 numbers.

（显然，包作者不认为 200 是要使用的 "sane" 数字，尽管查看 neig=50 和 neig-200 的结果值我没有看到可辨别的切割点第 50 个特征值。但是......他们必须在我向您展示了如何访问的代码中设置它。）

在 R 中绘制 RSSA 包中的 50 多个组件

plot more than 50 components in RSSA package in R

r

svd