summary.mvna {mvna} | R Documentation |
Summary method for mvna
objects. The function returns a list
containing the cumulative transition hazards, variance and other
informations.
## S3 method for class 'mvna' summary(object, level = 0.95, var.type = c("aalen", "greenwood"), ci.fun = c("log", "linear", "arcsin"), ...) ## S3 method for class 'mvna' print.summary(x, ...)
object |
An object of class |
level |
Level of the pointwise confidence interval. Default is 0.95. |
var.type |
Which of the "aalen" or "greenwood" variance estimator should be displayed and used to compute the pointwise confidence intervals. Default is "aalen". |
ci.fun |
Which transformation to apply to the confidence intervals. One of "linear", "log", "arcsin". Default is "log". |
... |
Further arguments. |
x |
An object of class |
Returns an object of class mvna
which is a list of data frames
named after the possible transitions. Each data frame contains the
following columns:
time |
Event times at which the cumulative hazards are estimated. |
na |
Estimated cumulative transition hazards. |
var.aalen or var.greenwood |
Variance estimates. The name depends
on the |
lower |
Lower bound of the pointwise confidence interval. |
upper |
Upper bound. |
n.risk |
Number of individuals at risk of experiencing an event just before t. |
n.event |
Number of transitions at time t. |
Arthur Allignol, arthur.allignol@gmail.com
data(sir.adm) ## data set transformation data(sir.adm) id <- sir.adm$id from <- sir.adm$pneu to <- ifelse(sir.adm$status==0,"cens",sir.adm$status+1) times <- sir.adm$time dat.sir <- data.frame(id,from,to,time=times) ## Possible transitions tra <- matrix(ncol=4,nrow=4,FALSE) tra[1:2,3:4] <- TRUE na.pneu <- mvna(dat.sir,c("0","1","2","3"), tra,"cens") summ.na.pneu <- summary(na.pneu) ## cumulative hazard for 0 -> 2 transition: summ.na.pneu$"0 2"$na