data(deermice)
Originally,
the columns of this data frame
represent the sex (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
deermice
has 38 deermice whereas
Perom
had 36 deermice
(Perom
has been withdrawn.)
In deermice
the two semi-adults have been classified as adults.
The sex
variable has 1 for female, and 0 for male.
Otis, D. L. et al. (1978) Statistical inference from capture data on closed animal populations, Wildlife Monographs, 62, 3--135.
posbernoulli.b
,
posbernoulli.t
,
fill1
.head(deermice)
fit1 <- vglm(cbind(y1, y2, y3, y4, y5, y6) ~ sex + age,
posbernoulli.t(parallel.t = TRUE), data = deermice, trace = TRUE)
coef(fit1)
coef(fit1, matrix = TRUE)
Run the code above in your browser using DataLab