progdismodel {changeLOS}R Documentation

progressive disability model

Description

Change in LOS and impact of an intermediate event on mortality can also be investigated in a so-called progressive disability model. This multi-state model can be described and the Aalen-Johansen estimator for transition probabilities can be computed. Further will be computed and plotted the mortality, attributable mortality and the population attributable mortality.

Usage

progdismodel(data, plot=FALSE, file1, file2, file3, max.time,
             lwd=2, cex=1.2, lty1=1, lty2=2, lty3=3,
             color1=4, color2=1, color3=2)

Arguments

data data.frame of the form data.frame(id, j.01, j.02, j.03, j.12, j.13, cens):
id:
id (patient id, admision id)
j.01:
observed time for jump from 0 to 1
j.02:
observed time for jump from 0 to 2
j.03:
observed time for jump from 0 to 3
j.12:
observed time for jump from 1 to 2
j.13:
observed time for jump from 1 to 3
cens:
censoring time (either in initial or intermediate state)
plot TRUE/FALSE
file1 ps file name for the first plot (mortality)
file2 ps file name for the second plot(attributable mortality)
file3 ps file name for the third plot (population attributable mortality)
max.time the last time point of interest
lwd the line widths
cex character expansion factor
lty1 the line type of line 1 of the first plot
lty2 the line type of line 2 of the second plot
lty3 the line type of line 3 of the third plot
color1 the color of line 1 of the first plot
color2 the color of line 2 of the second plot
color3 the color of line 3 of the third plot

Details

The first plot shows P( death, t ), P(death | risk factor absent,t) and P(death | risk factor present,t): mortality.

The second plot shows P(death | risk factor present,t) - P(death | risk factor absent,t): attributable mortality.

The third plot shows (P(death,t) - P(death | risk factor absent,t))/P(death,t): population attributable mortality.

Value

A list of

trans an object of class trans
aj an object of class aj
times.par
PAR population attributable mortality
AR attributable mortality
death P(death,t)
death.given.rfa P(death | risk factor absent,t)
death.given.rfp P(death | risk factor present,t)
tab table with the values at the last time point. if max.time is not missing, this will be the last time point.

References

Matthias Wangler mw@imbi.uni-freiburg.de

Examples

## Not run: 
 data(los.data)
 res <- progdismodel(los.data)
## End(Not run)

[Package changeLOS version 1.1.8 Index]