data(vaccine.trial)
print(with(vaccine.trial,
sensitivityJR(z=treatment,s=hiv.outcome,y=logVL,
beta0=c(-1,-.75,-.5,-.25,0,.25,.5,.75,1),
beta1=c(-1,-.75,-.5,-.25,0,.25,.5,.75,1),
phi=c(0.95,0.90), selection="infected",
groupings=c("placebo","vaccine"),
N.boot=50)
))
Run the code above in your browser using DataLab