# NOT RUN {
# A variation on example(rstar.ci):
# }
# NOT RUN {
obj <- rstar.ci(quinedata, thetainit=c(coef(quine.nb1),quine.nb1$theta), floglik=logLikNbin,
datagen=genDataNbin, fscore=gradLikNbin, fpsi=function(theta) theta[2], R=1000,
psidesc="Coefficient of EthN", npoints=5, away=0.1)
plot(obj)
# }
Run the code above in your browser using DataLab