如何对一组特定的因子类型变量值运行逻辑回归答案

【问题标题】：How can I run logistic regression on a specific group of factor type variable values如何对一组特定的因子类型变量值运行逻辑回归
【发布时间】：2014-04-24 15:02:25
【问题描述】：

我想对分类变量的特定组（范围 0f 值）进行逻辑回归。我执行了以下步骤： 1.我将变量分组：

cut_Var3 <- cut(dat$Var3,breaks=c(0,3,6,9))

table(cut_Var3) 的结果给了我这个输出（cut_Var3 变成了一个因素）：

# (0,3] (3,6] (6,9] 
#    5     4     4

我想对其他变量进行逻辑回归，但仅针对 (3,6) 的水平。所以我将能够对第二组的 4 个观察值进行回归。 2. 我尝试编写这行代码（以及其他变体）：

ff <- glm( TargetVar ~ relevel(cut_Var3,3:6), data = dat)

但没有运气。

我应该怎么做才能正常运行？

附上一个示例数据集：

dat <- read.table(text = " TargetVar  Var1    Var2       Var3
 0        0        0         7
 0        0        1         1
 0        1        0         3
 0        1        1         7
 1        0        0         5
 1        0        1         1
 1        1        0         0
 1        1        1         6
 0        0        0         8
 0        0        1         5
 1        1        1         4
 0        0        1         2
 1        0        0         9
 1        1        1         2  ", header = TRUE)

【问题讨论】：

标签： r logistic-regression

【解决方案1】：

对于relevel，您需要完全按照因子中出现的方式指定级别标签：

glm( TargetVar ~ relevel(cut_Var3,"(3,6]"), data = dat)

Call:  glm(formula = TargetVar ~ relevel(cut_Var3, "(3,6]"), data = dat)

Coefficients:
                    (Intercept)  relevel(cut_Var3, "(3,6]")(0,3]  
                           0.75                            -0.35  
relevel(cut_Var3, "(3,6]")(6,9]  
                          -0.50  

Degrees of Freedom: 12 Total (i.e. Null);  10 Residual
  (1 observation deleted due to missingness)
Null Deviance:      3.231 
Residual Deviance: 2.7  AIC: 24.46

【讨论】：