#Simulated competing risks data

n=1000
time1=rweibull(n,2)
time2=rweibull(n,.5)
censtime=runif(n)*2
obstime=pmin(time1,time2,censtime)
d = 1*(obstime==time1)+2*(obstime==time2)

#Estimation of cumulative cause-specific hazard

library(survival)
ch1 = survfit(Surv(obstime,d==1)~1,type="fh2")
ch2 = survfit(Surv(obstime,d==2)~1,type="fh2")
plot(ch1,fun="cumhaz",mark.time=F,main="Cum.haz. 1")
plot(ch2,fun="cumhaz",mark.time=F,main="Cum.haz. 2")

#Estimation of cumulative incidence functions in 'survival'
ci.surv = survfit(Surv(obstime,d,type="mstate")~1)
plot(ci.surv,main="Cum.incidence")

#Estimation of cumulative incidence functions in 'mstate'

library(mstate)
library(ggplot2)
ci.sim = Cuminc(time=obstime,status=d)
plot(x = ci.sim,use.ggplot = TRUE,conf.type = "log",lty = 1:2,conf.int = 0.95)