data(BBS2001)
desc<-descriptive(BBS2001,dfreq=TRUE,dtype="nbcap",t=50)
plot(desc)
# Note in this plot the convex shape typical of model Mh.
cp<-closedp.0(BBS2001,dfreq=TRUE,dtype="nbcap",t=50,t0=20)
cp
plot(cp)
# The Gamma estimator has the smallest deviance and the all
# its residuals are small.
# Let's calculate a profile confidence interval for the gamma
# estimator:
closedpCI.0(BBS2001,dfreq=TRUE,dtype="nbcap",t=50,t0=20,
m="Mh",h="Gamma",theta=3.5)
# One can also calculate a profile confidence interval for Chao's
# lower bound:
closedpCI.0(BBS2001,dfreq=TRUE,dtype="nbcap",t=50,t0=20,
m="Mh",h="Chao")
Run the code above in your browser using DataLab