#' @name print.impacts.SEMCMC
#' @rdname print.impacts.SEMCMC
#' @title Display summary of impacts from fitted with SEMCMC.
#'
#' @description This function will print summary of the impacts (direct, indirect and total) from
#' a SEMCMC object. Output is similar to spdep::impacts
#' @param x A SEMCMC.impacts object.
#' @param ... Extra argument to compute the impacts (not used).
#' @keywords spatial models
#' @export
#'
#' @importFrom stats sd
#'
#' @examples
#' data(columbus, package = "spdep")
#'
#' W <- spdep::nb2mat(col.gal.nb, style = "W")
#' m.form <- CRIME ~ INC + HOVAL
#' slm.mcmc <- SEMCMC(m.form, data = columbus, W = W, model = "slm")
#' impacts(slm.mcmc, W)
print.impacts.SEMCMC <- function(x, ...) {
#Lots of checks here
#Summary statistics
post.values <- lapply(x, function(X) {
apply(X, 2, function(Y) { c(mean(Y), sd(Y))})
})
#Posterior means
post.mean <- t(do.call(rbind, lapply(post.values, function(X){X[1,]})))
post.sd <- t(do.call(rbind, lapply(post.values, function(X){X[2,]})))
cat("Impact measures (posterior mean):\n")
print(post.mean)
cat("\nImpact measures (posterior st. dev.):\n")
print(post.sd)
return(list(post.mean = post.mean, post.sd = post.sd))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.