The goal of **calibmsm** is to provide a set of tools for estimating calibration plots when validating an existing (i.e. previously developed) multistate model. While the vignettes and articles for this package focus on assessing the calibration of a multistate model, calibration plots can also be produced for predicted risks from dynamic models, landmark super models, competing risks models and standard single outcome survival models, where predictions can be made at any landmark time.

For a detailed guide on how to use **calibmsm** please see the Overview vignette.

The package can be installed from CRAN as follows:

You can install the development version of calibmsm from GitHub with:

This is a basic example which shows you how to assess the calibration of the transition probabilities at 5 years follow up for individuals out of state `j = 1`

at time `s = 0`

using the BLR-IPCW approach. Please see the *Overview* vignette for examples of how to assess calibration using the pseudo-value and MLR-IPCW approaches, as well as details of the methodology.

The predicted transition probabilities are stored in `tps0`

, the individuals data are stored in `ebmtcal`

, and the data in `msdata`

format are stored in `msebmtcal`

. Calibration curves are estimated using `calib_msm`

using the BLR-IPCW method (`calib.type = 'blr'`

) with inverse probability of censoring weights are calculated based on variables `year`

, `age`

, `prophylaxis`

and donor gender `match`

. The calibration curves are estimated using restricted cubic splines with 3 knots. A 95% confidence interval is calculated using bootstrapping with 200 bootstrap replicates.

```
## Load calibmsm
library(calibmsm)
## Extract relevant predicted risks from tps0
tp.pred <- dplyr::select(dplyr::filter(tps0, j == 1), any_of(paste("pstate", 1:6, sep = "")))
## Calculate observed event probabilities
dat.calib.blr <-
calib_msm(data.mstate = msebmtcal,
data.raw = ebmtcal,
j = 1,
s = 0,
t = 1826,
tp.pred = tp.pred,
calib.type = "blr",
curve.type = "rcs",
rcs.nk = 3,
w.covs = c("year", "agecl", "proph", "match"),
CI = 95,
CI.R.boot = 200)
## Plot calibration plots
plot(dat.calib.blr, combine = TRUE, nrow = 2, ncol = 3)
```

If you encounter a bug, please file an issue with a minimal reproducible example on GitHub.