Extract summary statistics for the posterior correlation matrix which have been stored during sampling
runningcormat extracts summary statistics from the model-implied
correlation matrix
from an fsvdraws object for one point in time.
runningcormat(x, i, statistic = "mean", type = "cor")
x |
Object of class |
i |
A single point in time. |
statistic |
Indicates which statistic should be extracted. Defaults
to |
type |
Indicates whether covariance ( |
Matrix containing the requested correlation matrix summary statistic.
Other extractors:
cormat.fsvdraws(),
covmat.fsvdraws(),
runningcovmat()
set.seed(1)
sim <- fsvsim(n = 500, series = 3, factors = 1) # simulate
res <- fsvsample(sim$y, factors = 1, runningstore = 6) # estimate
cor100mean <- runningcormat(res, 100) # extract mean at t = 100
cor100sd <- runningcormat(res, 100, statistic = "sd") # extract sd
lower <- cor100mean - 2*cor100sd
upper <- cor100mean + 2*cor100sd
true <- cormat(sim, 100)[,,1] # true value
# Visualize mean +/- 2sd and data generating values
par(mfrow = c(3,3), mar = c(2, 2, 2, 2))
for (i in 1:3) {
for (j in 1:3) {
plot(cor100mean[i,j], ylim = range(lower, upper), pch = 3,
main = paste(i, j, sep = ' vs. '), xlab = '', ylab = '')
lines(c(1,1), c(lower[i,j], upper[i,j]))
points(true[i,j], col = 3, cex = 2)
}
}Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.