Следующие данные взяты из таблицы 7.28 в «Введение в категориальный анализ данных» (2-е издание). Я пытаюсь приспособить модель независимости, и я думаю, что она почти у меня есть, но по какой-то причине факторы TBC1 и религиозные1 отсутствуют в резюме. Как я могу заставить их появиться? Я использую неправильную регрессию? Это сделано в R studio.
religious=c(rep(1, 4), rep(2, 4), rep(3, 4), rep(4, 4), rep(5, 4), rep(6, 4), rep(7, 4), rep(8, 4), rep(9, 4))
TBC=c(rep(c(1, 2, 3, 4), 9))
count=c(49, 49, 19, 9, 31, 27, 11, 11, 46, 55, 25, 8, 34, 37, 19, 7, 21, 22, 14, 16, 26, 36, 16, 16, 8, 16, 15, 11, 32, 65, 57, 61, 4, 17, 16, 20)
relig.TBC=data.frame(religious, TBC, count)
relig.TBC$religious=factor(relig.TBC$religious, levels=sort(unique(relig.TBC$religious)))
relig.TBC$TBC=factor(relig.TBC$TBC, levels=sort(unique(relig.TBC$TBC)))
relig.TBC.glm=glm(count~religious+TBC, family=poisson, data=relig.TBC)
q=summary(relig.TBC.glm)
q