如何处理 R 中逻辑回归模型的 glm.fit 错误“NA/NaN/Inf”？

2024-05-21 22:45:08 103

当我们为逻辑回归模型创建一般线性模型时，我们需要将分布族指定为二项式。当我们不指定分布系列时，会出现错误“NA/NaN/Inf”。因此，在创建逻辑回归模型时，需要在glm函数中使用family="binomial"。

示例1

以下代码段创建了一个示例数据框-

iv1<-rpois(20,5)
iv2<-rpois(20,2)
iv3<-rpois(20,5)
DV1<-sample(0:1,20,replace=TRUE)
df1<-data.frame(iv1,iv2,iv3,DV1)
df1

创建以下数据框-

  iv1  iv2 iv3 DV1
1   5   2   6  0
2   3   1   3  1
3   3   4   8  1
4   5   3   3  1
5   8   2   6  1
6   3   1   4  0
7   6   1   8  1
8   3   1   7  0
9   9   2   6  0
10  7   2   4  0
11  6   4   5  1
12 12   2   4  1
13  6   2   2  0
14  5   1   3  0
15  4   1  10  0
16  3   3   4  0
17  4   1   6  1
18  9   3   4  1
19  7   1   3  1
20  4   3   4  0

要为df1中的数据创建逻辑回归模型，请将以下代码添加到上述代码段中-

Model_1<-glm(factor(DV1)~iv1+iv2+iv3,data=df1)

glm.fit(x=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,-

'y'中的NA/NaN/Inf

此外：警告消息-

1:In:'-'对因子没有意义Ops.factor(y,mu)

2:在:'-'对因子没有意义Ops.factor(eta,offset)

3:In:'-'对因子没有意义Ops.factor(y,mu)

要使用分布族为二项式为df1中的数据创建逻辑回归模型，请将以下代码添加到上述代码段中-

iv1<-rpois(20,5)
iv2<-rpois(20,2)
iv3<-rpois(20,5)
DV1<-sample(0:1,20,replace=TRUE)
df1<-data.frame(iv1,iv2,iv3,DV1)
Model_1<-glm(factor(DV1)~iv1+iv2+iv3,data=df1,family="binomial")
summary(Model_1)

输出

如果您将上述所有代码作为单个程序执行，它会生成以下输出-

Call:
glm(formula = factor(DV1) ~ iv1 + iv2 + iv3, family = "binomial",
data = df1)

Deviance Residuals:
   Min    1Q      Median    3Q     Max
-1.61472 -1.05484 -0.07657 1.07422 1.71351

Coefficients:
             Estimate  Std. Error z value Pr(>|z|)
(Intercept) -2.59874   2.15616  -1.205   0.228
iv1          0.26684   0.22055   1.210   0.226
iv2          0.38736   0.47527   0.815   0.415
iv3          0.06822   0.23316   0.293   0.770

(Dispersion parameter for binomial family taken to be 1)

Null deviance: 27.726 on 19 degrees of freedom
Residual deviance: 25.223 on 16 degrees of freedom
AIC: 33.223

Number of Fisher Scoring iterations: 4

示例2