Print method for PMMH output
Usage
# S3 method for class 'pmmh_output'
print(x, ...)
Examples
# Create dummy chains for two parameters across two chains
chain1 <- data.frame(param1 = rnorm(100), param2 = rnorm(100), chain = 1)
chain2 <- data.frame(param1 = rnorm(100), param2 = rnorm(100), chain = 2)
dummy_output <- list(
theta_chain = rbind(chain1, chain2),
diagnostics = list(
ess = c(param1 = 200, param2 = 190),
rhat = c(param1 = 1.01, param2 = 1.00)
)
)
class(dummy_output) <- "pmmh_output"
print(dummy_output)
#> PMMH Results Summary:
#> Parameter Mean SD Median CI Lower.2.5% CI Upper.97.5% ESS Rhat
#> param1 -0.01 0.99 0.05 -1.96 1.69 200 1.01
#> param2 -0.02 0.89 -0.01 -1.59 1.77 190 1.00