1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
|
\name{cbpp}
\alias{cbpp}
\docType{data}
\title{Serological Incidence of Contagious Bovine Pleuropneumonia in Ethiopia}
\description{
Contagious bovine pleuropneumonia (CBPP) is a major disease of cattle
in Africa, caused by a mycoplasma. This dataset describes the
serological incidence of CBPP in zebu cattle during a follow-up survey
implemented in 15 commercial herds located in the Boji district
of Ethiopia. The goal of the survey was to study the within-herd spread
of CBPP in newly infected herds. Blood samples were quarterly
collected from all animals of these herds to determine their CBPP
status. These data were used to compute the serological incidence of
CBPP (new cases occurring during a given time period). Some data are
missing (lost to follow-up).
}
\usage{data(cbpp)
}
\format{
A data frame with 56 observations on the following 4 variables.
\describe{
\item{\code{herd}}{A factor identifying the herd (1 to 15).}
\item{\code{incidence}}{The number of new serological cases for a
given herd and time period.}
\item{\code{size}}{A numeric vector describing herd size at the
beginning of a given time period.}
\item{\code{period}}{A factor with levels \code{1} to \code{4}.}
}
}
\details{
Serological status was determined using a competitive enzyme-linked
immuno-sorbent assay (cELISA).
}
\source{
Lesnoff, M., Laval, G., Bonnet, P., Abdicho, S., Workalemahu, A.,
Kifle, D., Peyraud, A., Lancelot, R., Thiaucourt, F. (2004)
Within-herd spread of contagious bovine pleuropneumonia in Ethiopian highlands.
\emph{Preventive Veterinary Medicine} \bold{64}, 27--40.
}
\examples{
## response as a matrix
(m1 <- lmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
family = binomial, data = cbpp))
## response as a vector of probabilities and usage of argument "weights"
m2 <- lmer(incidence / size ~ period + (1 | herd), weights = size,
family = binomial, data = cbpp)
## Confirm that these are equivalent:
stopifnot(all.equal(coef(m1), coef(m2)),
all.equal(ranef(m1), ranef(m2)))
}
\keyword{datasets}
|
