| read sample | | samp<−c(rep(0, 8), rep(1, 21), rep(2, 27), rep(3, 39), rep(4, 29), rep(5, 28), rep(6, 16), | | rep(7, 15), rep(8, 6), rep(9, 5), rep(10, 2), rep(11, 1), rep(14, 3)) | | MLE estimate | | K<−14 | | hatrhoMLE<−MLERoMM1K(K, samp) | | Bayesian estimate | | a<−1.0 | | b<−1.0 | | set.seed(13579) | | hatrhoBayes<−EBaRoMM1K(K, samp, a, b) | | Bayesian Jeffreys estimate | | set.seed(13579) | | hatrhoJe<−EJeRoMM1K(K, samp) | | Bootstrap corrected estimate | | set.seed(13579) | | hatrhoBoot<−EBoRoMM1K(K, samp) | | c(hatrhoMLE, hatrhoBayes, hatrhoJe, hatrhoBoot) | | > 0.8405396 0.8403625 0.8403625 0.8405683 |
|
