Description Usage Arguments Value References Examples

Compute the (integrated) Brier score to evaluate the model fit for interval-censored survival data.

1 2 3 4 5 6 |

`obj` |
an object of class |

`pred` |
predicted values. This can be a matrix of survival probabilities evaluated
at a sequence of time points for a set of new data, a list of |

`btime` |
a vector of length two indicating the range of times that the scores are computed on.
The default |

`type` |
a character string denoting the type of scores returned. For |

If `type = "IBS"`

, this returns the integrated Brier score.

If `type = "BS"`

, this returns the Brier scores.

S. Tsouprou. Measures of discrimination and predictive accuracy for interval-censored data. Master thesis, Leiden University. https://www.math.leidenuniv.nl/scripties/MasterTsouprou.pdf.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 | ```
### Example with dataset miceData
library(survival)
library(icenReg)
data(miceData)
## For proper evaluation, Inf should be set to be a large number, for example, 9999999.
idx_inf <- (miceData$u == Inf)
miceData$u[idx_inf] <- 9999999.
obj <- Surv(miceData$l, miceData$u, type = "interval2")
## Model fit for an NPMLE survival curve with survfit
pred <- survival::survfit(formula = Surv(l, u, type = "interval2") ~ 1, data = miceData)
# Integrated Brier score up to time = 642
sbrier_IC(obj, pred, btime = c(0, 642), type = "IBS")
## Model fit for a semi-parametric model with icenReg::ic_sp()
pred <- icenReg::ic_sp(formula = Surv(l, u, type = "interval2") ~ 1, data = miceData)
# Integrated Brier score up to the largest endpoints of all censoring intervals in the dataset
sbrier_IC(obj, pred, type = "IBS")
## Model fit for an NPMLE survival curve with icenReg::ic_np()
pred <- icenReg::ic_np(miceData[,c('l', 'u')])
# Brier score computed at every left and right endpoints of all censoring intervals in the dataset
sbrier_IC(obj, pred, type = "BS")
``` |

ICcforest documentation built on Feb. 17, 2020, 9:07 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.