There's a bug when we try to specify formula_time = ~1 in fitSmoothHazard.fit. This would happen, for example, if we assumed an exponential hazard.
library(casebase)
#> See example usage at http://sahirbhatnagar.com/casebase/N<-1000; p<-30nzc<-0.33*px<-matrix(rnorm(N*p), N, p)
dimnames(x)[[2]] <- paste0("x", seq_len(p))
beta<- rnorm(nzc)
fx<-x[, seq(nzc)] %*%
Code and results of Section 4 of the paper "Fine-Gray subdistribution hazard models to simultaneously estimate the absolute risk of different event types: cumulative total failure probability may exceed 1", by Peter Austin, Ewout Steyerberg & Hein Putter
Code repository for the manuscript 'External validation of the performance of competing risks prediction models: a guide through modern methods' (in preparation)
Code and supplementary materials for the manuscript "Multiple imputation for cause-specific Cox models: assessing methods for estimation and prediction"
There's a bug when we try to specify
formula_time = ~1infitSmoothHazard.fit. This would happen, for example, if we assumed an exponential hazard.