Reputation: 634
I'm trying to replicate the eigenvector plotting image from Rssa package
set.seed(1)
x <- cumsum(rnorm(200))
library(Rssa)
s <- ssa(x)
r <- reconstruct(s)
plot(s, t="paired")
But I get something completely different. What am I doing wrong?
par(mfrow = c(4,4), mar = c(2, 2, 2, 1))
n_comp = 5
combn(n_comp, 2) |>
apply(MARGIN = 2, \(i) plot(x= r[[ i[1] ]],
y= r[[ i[2] ]],
type = "l"))
Upvotes: 0
Views: 15