Description Usage Arguments Value See Also Examples
Plot estimated baseline transition intensities from an object of class
idm
optionally with confidence limits.
1 2 3 4 |
x |
a |
conf.int |
If TRUE show confidence limits |
citype |
Type of confidence limits, can be "shadow" or "bars" |
add |
If TRUE add to existing plot |
axes |
If TRUE axes are drawn |
col |
Color of the lines |
lwd |
Width of the lines |
lty |
Type of the lines |
xlim |
Limits for x-axis |
ylim |
Limits for y-axis |
xlab |
Label for x-axis |
ylab |
Label for y-axis |
legend |
If TRUE a legend is drawn, which can be further controlled via |
transition |
Choose one of the transition intensities: |
... |
Passed to |
Print a plot of the baseline transition intensities of an illness-death model estimated using a Weibull approach.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | library(lava)
library(prodlim)
m <- idmModel(scale.lifetime=1/10,scale.illtime=1/8)
distribution(m,"X") <- binomial.lvm()
regression(m,latent.lifetime~X) <- 0.7
set.seed(30)
d <- sim(m,100)
fit.weib <- idm(formula02=Hist(observed.lifetime,event=seen.exit)~1,
formula01=Hist(time=list(L,R),event=seen.ill)~1,data=d,conf.int=FALSE)
plot(fit.weib)
## Not run:
## FIXME: the limits for the 01 transition are a bit wide!?
## with bootstrap confidence limits
fit.weib <- idm(formula02=Hist(observed.lifetime,event=seen.exit)~1,
formula01=Hist(time=list(L,R),event=seen.ill)~1,data=d,conf.int=TRUE)
plot(fit.weib)
## End(Not run)
|
Loading required package: prodlim
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.