data(multcif) # simulated data
multcif$cause[multcif$cause==0] <- 2
times=seq(0.05,3,by=0.1) # to speed up computations use only these time-points
add<-comp.risk(Surv(time,status>0)~const(X)+cluster(id),data=multcif,
multcif$cause,n.sim=0,times=times)
###
out1<-cor.cif(add,data=multcif,cause1=1,cause2=1,theta=log(2+1))
summary(out1)
pad <- predict(add,X=1,Z=0,se=0,uniform=0)$P1
summary(out1,marg.cif=pad)
Run the code above in your browser using DataLab