以下数据来自“分类数据分析简介”（第2版）中的表7.28。我试图拟合一个独立模型，但我想我几乎拥有了它，但是由于某种原因，摘要中缺少了TBC1和虔诚的因素。如何让他们出现？我是否使用了错误的回归？这是在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