- Code: Select all
buckrd <- convert.inp("d:/RD_buckA.inp",
covariates = c("spec2", "sex", "trt2"))
buckrd$spec2 <- factor(buckrd$spec2,
levels = c(0, 1),
labels = c("msp","pman"))
buckrd$sex <- factor(buckrd$sex,
levels = c(0, 1),
labels = c("female", "male"))
buckrd$trt2 <- factor(buckrd$trt2,
levels = c(0, 1),
labels = c("placebo", "vaccine"))
#Robust
buckrd.process=process.data(data=buckrd,model="Robust", time.intervals=c(0,0,0,0,0.867,0,0,0,0,0.767,0,0,0),
groups=c("sex","spec2","trt2"))
#add 2 age bins
buckrd.ddl=make.design.data(buckrd.process)
buckrd.ddl=add.design.data(buckrd.process,buckrd.ddl,parameter="S",type="age",
bins=c(0,.5,3),name="tages")
buckrd.models=function()
{
#survival
S.tages.adtrt2=list(formula=~tages+trt2)
S.spec2.trt2=list(formula=~spec2*trt2)
S.trt2=list(formula=~trt2)
S.spec2.adtages=list(formula=~spec2+tages)
S.spec2=list(formula=~spec2)
S.tages=list(formula=~tages)
S.null=list(formula=~1)
#Gamma
#random
GammaDoublePrime.random=list(formula=~1,share=TRUE) # gamma' = gamma''
#random time
GammaDoublePrime.randomt=list(formula=~time,share=TRUE)
#no emigration
GammaDoublePrime.dot=list(formula=~1,share=TRUE,fixed=0) # gamma'=gamma''=0
p.null=list(formula=~1)
p.sex=list(formula=~sex, share=TRUE)
p.time=list(formula=~time, share=TRUE)
p.tages=list(formula=~tages, share=TRUE)
p.spec2=list(formula=~spec2, share=TRUE)
p.trt2=list(forumala=~trt2, share=TRUE)
p.session=list(forumala=~session, share=TRUE)
cml=create.model.list("Robust")
results=mark.wrapper(cml,data=buckrd.process,ddl=buckrd.ddl)
return(results)
}
buckrd.results=buckrd.models()
Thank you