116
117
bg.data = data.frame(y = 0:11, wts = c(227,123,72,42,21,31,11,14,6,4,7,28))
117
fit = vglm(y ~ 1, fam=betageometric, bg.data, weight=wts, trace=TRUE)
118
fitg = vglm(y ~ 1, fam= geometric, bg.data, weight=wts, trace=TRUE)
118
fit = vglm(y ~ 1, betageometric, bg.data, weight = wts, trace = TRUE)
119
fitg = vglm(y ~ 1, geometric, bg.data, weight = wts, trace = TRUE)
124
125
# Very strong evidence of a beta-geometric: