EAlasso (version 0.1.0)

summary.MHLS: Summarizing Metropolis-Hastings sampler outputs

Description

Summary method for class "MHLS".

Usage

# S3 method for MHLS
summary(object, ...)

Arguments

object

an object of class "MHLS", which is a result of MHLS.

...

additional arguments affecting the summary produced.

Value

mean, median, standard deviation, 2.5% quantile and 97.5% quantile for each beta and its subgradient.

Details

This function provides a summary of each sampled beta and subgradient.

Examples

# NOT RUN {
#' set.seed(123)
n <- 10
p <- 5
X <- matrix(rnorm(n * p), n)
Y <- X %*% rep(1, p) + rnorm(n)
sigma2 <- 1
lbd <- .37
weights <- rep(1, p)
LassoResult <- Lasso.MHLS(X = X, Y = Y, lbd = lbd, type = "lasso", weights = weights)
B0 <- LassoResult$B0
S0 <- LassoResult$S0
summary(MHLS(X = X, PE = rep(0, p), sig2 = sigma2, lbd = lbd,
     weights = weights, B0 = B0, S0 = S0, niter = 50, burnin = 0,
     type = "coeff"))
# }