| sit.adm {mvna} | R Documentation |
Pneumonia status on admission for intensive care unit (ICU) patients, a random sample from the SIR-3 study.
data(sir.adm)
The data contains 765 rows and 4 variables:
Beyersmann, J., Gastmeier, P., Grundmann, H., Baerwolff, S., Geffers, C., Behnke, M., Rueden, H., and Schumacher, M. Use of multistate models to assess prolongation of intensive care unit stay due to nosocomial infection. Infection Control and Hospital Epidemiology, 27:493-499, 2006.
# 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)
diag(tra) <- TRUE
tra[1:2,3:4] <- TRUE
na.pneu <- mvna(dat.sir,c("0","1","2","3"),
tra,"cens")
xyplot(na.pneu,tr.choice=c("0 2","1 2","0 3","1 3"),
aspect=1,strip=strip.custom(bg="white",
factor.levels=c("No pneumonia on admission -- Discharge",
"Pneumonia on admission -- Discharge",
"No pneumonia on admission -- Death",
"Pneumonia on admission -- Death"),
par.strip.text=list(cex=0.9)),
scales=list(alternating=1),xlab="Days",
ylab="Nelson-Aalen esimates")