summary.mvna: Summary method for objects of class 'mvna'

Description Usage Arguments Value Author(s) See Also Examples

Description

Summary method for mvna objects. The function returns a list containing the cumulative transition hazards, variance and other informations.

Usage

1
2
3
4
5
6
7
## 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, ...)

Arguments

object

An object of class mvna.

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 summary.mvna.

Value

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 var.type argument. Default will be var.aalen.

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.

Author(s)

Arthur Allignol, arthur.allignol@uni-ulm.de

See Also

mvna

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
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

aallignol/mvna documentation built on May 10, 2019, 3:13 a.m.