R Commander中GLM建模报错求助：二项变量建模失败

问题描述 #

我正在完成一项质量控制相关作业，需通过广义线性模型（GLM）分析哪些定性变量会影响响应变量。我使用R Commander工具建模，将取值为S和N的二项定性变量作为响应变量，其余变量也均为定性变量，点击确认后出现报错。手动编写脚本：

GLM.1 <- glm(A.demanda.foi.atendida.e.foi.no.prazo. ~ UF, family=binomial(logit), data=Tabela_Procon_última)

同样失败，收到以下错误/警告信息：

RcmdrMsg: [5] ERROR: y values ​​must be 0 <= y <= 1

RcmdrMsg: [6] ERROR: Error evaluating argument 'object' in method selection for function 'summary'

Warning in sprintf(gettextRcmdr("the model %s is no longer available"), : one argument not used by format 'The model is no longer available'

RcmdrMsg: [7] ERROR: Templates are no longer available

RcmdrMsg: [8] ERROR: object 'GLM.1' not found

原因分析 #

• 核心错误：响应变量格式不匹配：glm的二项分布模型要求响应变量是0/1数值型，或是能被识别为二项分类的因子型。你的响应变量是取值为"S"和"N"的字符型（或未正确转换的因子），R无法将其映射到0-1范围，导致第一个报错，后续所有错误都是模型创建失败引发的连锁反应。
• 连锁报错逻辑：由于GLM.1对象未成功生成，后续调用summary等函数时找不到目标对象，进而触发R Commander的一系列关联报错。

解决方案 #

• 转换响应变量为0/1数值型：
  使用ifelse函数将"S"映射为1，"N"映射为0：

  Tabela_Procon_última$resposta_binaria <- ifelse(Tabela_Procon_última$A.demanda.foi.atendida.e.foi.no.prazo. == "S", 1, 0)
  

  用转换后的变量重新建模：

  GLM.1 <- glm(resposta_binaria ~ UF, family=binomial(logit), data=Tabela_Procon_última)
  

• 或转换响应变量为二项因子：
  将响应变量转换为因子并指定参考水平，glm会自动适配二项模型要求：

  Tabela_Procon_última$A.demanda.foi.atendida.e.foi.no.prazo. <- factor(Tabela_Procon_última$A.demanda.foi.atendida.e.foi.no.prazo., levels = c("N", "S"))
  

  之后直接运行原建模脚本即可。
• 验证模型有效性：
  模型创建后，先输入GLM.1确认对象存在，再调用summary(GLM.1)查看模型结果。

内容的提问来源于stack exchange，提问作者Flora Rodrigues