inst/userguide/figures/CS1--plotkemstates.R
In nwfsc-timeseries/MARSS: Multivariate Autoregressive State-Space Modeling
###################################################
### code chunk number 11: plotkemstates
###################################################
par(mfrow = c(3, 3))
for (i in 1:9) {
plot(years, y.tss[i, ], xlab = "", ylab = "Index of log abundance", lwd = 2, bty = "l")
lines(years, x.tss[i, ], type = "l", lwd = 2, lty = 2)
lines(years, kem.states[i, ], type = "l", col = 2, lwd = 1, lty = 1)
title(paste("simulation ", i))
}
# legend("topright", c("Observed","True","KalmanEM estimate"),lty = c(-1, 2, 1), pch = c(1, -1, -1),col=c(1,1,2))
nwfsc-timeseries/MARSS documentation built on June 3, 2023, 1:32 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
###################################################
### code chunk number 11: plotkemstates
###################################################
par(mfrow = c(3, 3))
for (i in 1:9) {
plot(years, y.tss[i, ], xlab = "", ylab = "Index of log abundance", lwd = 2, bty = "l")
lines(years, x.tss[i, ], type = "l", lwd = 2, lty = 2)
lines(years, kem.states[i, ], type = "l", col = 2, lwd = 1, lty = 1)
title(paste("simulation ", i))
}
# legend("topright", c("Observed","True","KalmanEM estimate"),lty = c(-1, 2, 1), pch = c(1, -1, -1),col=c(1,1,2))
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.