подготовить более 50 компонентов в пакете RSSA в R - PullRequest
Новая
       25

подготовить более 50 компонентов в пакете RSSA в R

0 голосов
/ 28 июня 2018 
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')

но это не сработало!

Пример: аналогично этому случаю доступ к собственным значениям в пакете RSSA в R

1 Ответ

0 голосов
/ 28 июня 2018

При просмотре страницы справки для ?ssa мы видим параметр с именем neig, который задокументирован как;

целое число, количество желательных собственных пятен. Если «NULL», то будет использоваться правильное значение по умолчанию, см. «Подробности»

Используя это как именованный параметр: 

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

enter image description here

И 

> 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 «нормальным» числом для использования, хотя, глядя на значения результатов от neig = 50 и neig-200, я не вижу заметной точки отсечения для 50-го собственного значения. Но ... они должны установить это в коде, который я показал вам, как получить доступ.)

...