## Not run:
# set.seed(1)
# sim <- fsvsim(n = 500, series = 3, factors = 1) # simulate
# res <- fsvsample(sim$y, factors = 1, keeptime = "all") # estimate
# covs <- covmat(res) # extract
#
# # Trace plot of determinant of posterior covariance matrix
# # at time t = n = 500:
# detdraws <- apply(covs[,,,500], 3, det)
# ts.plot(detdraws)
# abline(h = mean(detdraws), col = 2) # posterior mean
# abline(h = median(detdraws), col = 4) # posterior median
# abline(h = det(covmat(sim)[,,500]), col = 3) # implied by DGP
#
# # Trace plot of draws from posterior covariance of Sim1 and Sim2 at
# # time t = n = 500:
# ts.plot(covs[1,2,,500])
# abline(h = covmat(sim)[1,2,500], col = 3) # "true" value
#
# # Smoothed kernel density estimate:
# plot(density(covs[1,2,,500], adjust = 2))
#
# # Summary statistics:
# summary(covs[1,2,,500])
# ## End(Not run)
Run the code above in your browser using DataLab