GLM中的奇怪错误
问题描述:
我正在创建多个泊松回归,每个模型都有不同的结果。在这个与其他人没有区别的特殊结果(所有二进制,无论是0还是1)上,我都会遇到这个我从未见过的奇怪错误。GLM中的奇怪错误
mix_1<- glm(count~ offset(log(time)) + age + gender, data=x[x$diagnosis=="match",], family=poisson)
glm模型没有问题,我没有按预期工作。尝试执行cbind和指数化系数时收到错误。
mix_2 <- cbind(exp(coef(mix_1)), exp(confint(mix_1)))
这是我的错误:
Waiting for profiling to be done...
Error in glm.fit(x = Xi, y = Y, weights = W, etastart = LP, offset = o, :
NA/NaN/Inf in 'x'
warning
Warning messages:
1: glm.fit: fitted rates numerically 0 occurred
我没有想到的是数据的一个子集将为这个问题是必要的,所以我没有提供一个。但如果我真的不需要,只需告诉我,我将创建一个子集。
所有帮助非常感谢!
答
我刚刚遇到同样的问题。在Andrew Gelman的博客here上有关MASS和glm.fit的问题进行了一些讨论,这些问题可能至少略有关联。
该问题似乎与confint()
函数。从confint手册:
存在用于其调用那些包MASS其基于简档 似然性类“GLM”和“NLS” 存根方法。
...
confint(glm.D93) # needs MASS to be present on the system
confint.default(glm.D93) # based on asymptotic normality
我发现,如果我使用confint.default()
功能问题消失。我希望我有更多的答案。