
data(Perom)
m
or f
),
the ages (y
: young, sa
: semi-adult, a
: adult),
the weights in grams, and the
capture histories of 38 individuals over 6 trapping
occasions (1: captured, 0: not captured).
The data set was collected by V. Reid and distributed
with the
Otis, D. L. et al. (1978) Statistical inference from capture data on closed animal populations, Wildlife Monographs, 62, 3--135.
posbernoulli.t
.head(Perom)
fit1 <- vglm(cbind(y1, y2, y3, y4, y5, y6) ~ sex + age,
posbernoulli.t(parallel.t = TRUE), data = Perom, trace = TRUE)
coef(fit1)
coef(fit1, matrix = TRUE)
Run the code above in your browser using DataLab