Displaying 1 result from an estimated 1 matches for "bacteria1".
Did you mean:
bacteria
2004 Nov 01
1
GLMM
...) <- structure (contr.sdif (3),
dimnames = list (NULL, c ("drug", "enc")))
bacteria.pql <- glmmPQL (y ~ trt + I(week > 2), random = ~ 1 |
ID, family = binomial, data = bacteria)
summary (bacteria.pql)
library (repeated)
y1 <- 1 * (bacteria$y == "y")
bacteria1 <- data.frame (bacteria, y1 = y1)
bacteria.glmm <- glmm (y1 ~ trt + I(week > 2), nest = ID, data
= bacteria1, family = binomial)
summary (bacteria.glmm)
library (glmmML)
bacteria.glmmML <- glmmML (y1 ~ trt + I(week > 2), cluster =
bacteria1$ID, data = bacteria1, family = binomi...