没有线性预测变量的Poisson GLM模型

时间:2019-08-09 19:26:21

标签: r linear-regression

我正在尝试在R中运行代码(这是我的新手),并且给了我一个非常大的数据集,我需要使用它来拟合泊松glm,例如log(mu)= B0 + B1x1。设Yi为对象i的响应计数,xi = 1(黑色),xi = 0(白色)。

可以在www.stat.ufl.edu/~aa/glm/data中找到数据集。

我加载了数据,但是很难理解该模型。 这是我到目前为止的代码,但显然我遗漏了一些东西。

str(hdata)
head(hdata)
attach(hdata)
hfit = glm(count ~ factor(race), family = poisson(link = log))
summary(hfit)

#plot the model
par(mfrow = c(2,2))
plot(hfit)

#overdispersion test
library(AER)
dispersiontest(hfit, trafo =1)

#goodness of fit test
sum(resid(hfit, type="pearson")^2)
#pvalue
1 - pchisq(2279.873, 1306)

我需要此模型的帮助,因为我似乎无法分开每个种族,我认为这是我需要做的。当我运行hfit的摘要时,我得到截距-2.38和factor(race)1的1.73。 AIC为1122。此外,当我进行超分散测试时,我得到c = 0.743,如果模型具有等分散,则c =0。对吗?谢谢

0 个答案:

没有答案