data(catb)
desc<-descriptive(catb,dfreq=TRUE,dtype="nbcap",t=11)
plot(desc)
# The heterogeneity graph is mildly convex, mostly because of
# 2 sites where the catbirds have been seen six times.
closedp.0(catb,dfreq=TRUE,dtype="nbcap",t=11)
# The residuals are OK. The occupancy rate (N/50) estimates
# varies between 40.6\% for M0 and 49.4\% for Mh Gamma.
# One could also try fitting M0 removing the 2 unusual sites:
closedp.0(catb,dfreq=TRUE,dtype="nbcap",t=11,t0=5)
# M0 fits very well and the occupancy rate estimate is 42.6\%.
Run the code above in your browser using DataLab