BayesMFSurv (version 0.1.0)

mfsurv.summary: mfsurv.summary()

Description

Returns a summary of a mfsurv object via summary.mcmc.

Usage

mfsurv.summary(object, parameter = c("betas", "gammas", "lambda"))

Arguments

object

an object of class mfsurv, the output of mfsurv.

parameter

one of three parameters of the mfsurv output. Indicate either "betas", "gammas" or "lambda".

Value

list. Empirical mean, standard deviation and quantiles for each variable.

Examples

# NOT RUN {
set.seed(95)
bgl <- Buhaugetal_2009_JCR
bgl <- subset(bgl, coupx == 0)
bgl <- na.omit(bgl)
Y   <- bgl$Y
X   <- as.matrix(cbind(1, bgl[,1:7]))
C   <- bgl$C
Z1  <- matrix(1, nrow = nrow(bgl))
Y0  <- bgl$Y0
model1 <- mfsurv(Y ~ X | C ~ Z1, Y0 = Y0,
                N = 50,
                burn = 20,
                thin = 15,
                w = c(0.1, .1, .1),
                m = 5,
                form = "Weibull",
                na.action = 'na.omit')

mfsurv.summary(model1, "betas")
# }