例题的话,就是这个啦——
library(survival)
data(bcdata)
bcdata<-data.frame(bcdata)#mustbeadataframe
try<-ICBayes(Surv(L,R,type=’interval2’)~x1,data=bcdata,
model=’case2ph’,status=bcdata[,3],p=1)
#generalform
try2<-ICBayes(model=’case2ph’,L=bcdata[,1],R=bcdata[,2],status=bcdata[,3],
xcov=bcdata[,4],p=1)
ICBayes包,应该就是做这个的~
链接在这里:
http://sph.sc.edu/epid_bios/cai/ICBayes-manual.pdf