progdismodel {changeLOS} | R Documentation |
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.
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)
data |
data.frame of the form data.frame(id, j.01, j.02, j.03, j.12, j.13, cens):
|
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 |
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.
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. |
Matthias Wangler mw@imbi.uni-freiburg.de
## Not run: data(los.data) res <- progdismodel(los.data) ## End(Not run)