# NOT RUN {
petersen(LRH01, M=120)
abundance(LRH01, n=120, model="H")
abundance(LRH01, n=120, model="I")
abundance(LRH01, n=120, model="T", B=2)
# The threshold model with B >= 3 is equivalent
# to the homogeneity model for this data set
# because max(LRH01$xi)=2
# }
Run the code above in your browser using DataLab