我需要使用R Studio为以下数据拟合线性乘线性模型。我已经安装了独立模型，但不确定如何通过线性模型生成线性模型。 （“分类数据分析简介”第二版的表7.28）
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(link=log), data=relig.TBC)
q=summary(relig.TBC.glm)
q