print.MHLS: Print Metropolis-Hastings sampler outputs

In seunghyunmin/EAinference: Estimator Augmentation and Simulation-Based Inference

Description

Print a brief summary of the MH sampler outputs.

Usage

## S3 method for class 'MHLS'
print(x, ...)

Arguments

x	an object of class "MHLS", which is a result of MHLS.
...	additional print arguments.

Details

print.MHLS prints out last 10 iterations and a brief summary of the simulation; number of iterations, number of burn-in periods, PE, PEtype and acceptance rate.

Value

Above results are silently returned.

Examples

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 <- lassoFit(X = X, Y = Y, lbd = lbd, type="lasso", weights = weights)
B0 <- LassoResult$B0
S0 <- LassoResult$S0
Result <- MHLS(X = X, PE = rep(0, p), sig2 = sigma2, lbd = lbd, group = 1:p,
     weights = weights, B0 = B0, S0 = S0, niter = 50, burnin = 0,
     type = "coeff")
print(Result)

seunghyunmin/EAinference documentation built on May 9, 2019, 5:58 p.m.