tests/glm.R

In free1way.docreg: Additional Documentation and Regression Tests for `stats::free1way'

library("free1way.docreg")
options(digits = 5)
set.seed(29)

N <- 500
x <- gl(2, N, labels = c("A", "B"))
y <- rlogis(length(x), location = c(0, 2)[x])
y <- factor(y < 0)

ci <- confint(m <- glm(y ~ x, family = binomial()), test = "LRT")
c(coef(m)["xB"], ci["xB",])
ci <- confint(m <- glm(y ~ x, family = binomial()), test = "Rao")
c(coef(m)["xB"], ci["xB",])

ft <- free1way(y = y, groups = x, B = 10000)
coef(ft)
print(ft, test = "LRT")
confint(ft, test = "LRT")
print(ft, test = "Rao")
confint(ft, test = "Rao")
print(ft, test = "Wald")
confint(ft, test = "Wald")
print(ft, test = "Permutation")
confint(ft, test = "Permutation")

y <- relevel(y, levels(y)[2])
ci <- confint(m <- glm(y ~ x, family = binomial(link = "cloglog")), test = "LRT")
c(coef(m)["xB"], ci["xB",])
ci <- confint(m <- glm(y ~ x, family = binomial(link = "cloglog")), test = "Rao")
c(coef(m)["xB"], ci["xB",])

ft <- free1way(y = y, groups = x, link = "loglog", B = 10000)
coef(ft)
print(ft, test = "LRT")
confint(ft, test = "LRT")
print(ft, test = "Rao")
confint(ft, test = "Rao")
print(ft, test = "Wald")
confint(ft, test = "Wald")
print(ft, test = "Permutation")
confint(ft, test = "Permutation")