开始使用免费开始使用

拟合麻雀存活概率模型

在本练习中,您将根据麻雀的体貌特征,估计其在严冬风暴中存活的概率。数据集 sparrow 已预加载。要预测的结果是 status("Survived","Perished")。我们将使用的变量有:

  • total_length:从喙尖到尾尖的全长(mm)
  • weight:重量(克)
  • humerus:肱骨长度(将翅膀连接到身体的"上臂骨")(英寸)

请记住,使用 glm()文档)创建逻辑回归模型时,必须显式指定 family = binomial

glm(formula, data = data, family = binomial)

您将调用 summary()broom::glance() 来查看检视逻辑回归模型的不同函数。您将查看的诊断之一是与 \(R^2\) 类似的指标,称为伪 $R^2$。

$$ pseudoR^2 = 1 - \frac{deviance}{null.deviance} $$

可以将 deviance 理解为类似方差的概念:它衡量分类数据的变异程度。伪 \(R^2\) 类似于标准回归中的 $R^2$:\(R^2\) 衡量模型"解释的方差"。而伪 \(R^2\) 则衡量"解释的偏差(deviance)"。

本练习是课程的一部分

R 中的监督学习:回归

查看课程

练习说明

  • 按视频中的建议,使用 TRUEFALSE 作为预测结果。在数据框 sparrow 中创建新列 survived,当 status == "Survived" 时取 TRUE。
  • 创建公式 fmla,将 survived 表达为所选变量的函数。打印该公式。
  • 拟合一个逻辑回归模型来预测麻雀的存活概率。将模型赋给变量 sparrow_model
  • 调用 summary() 查看模型系数、deviance 和 null deviance。
  • 对模型调用 glance(),以数据框形式查看 deviance 和其他诊断信息。将 glance() 的输出赋给变量 perf
  • 计算伪 $R^2$。

交互式实操练习

通过完成这段示例代码来试试这个练习。

# sparrow is available
summary(sparrow)

# Create the survived column
sparrow$survived <- ___

# Create the formula
(fmla <- _____)

# Fit the logistic regression model
sparrow_model <- ___

# Call summary
___

# Call glance
(perf <- ___)

# Calculate pseudo-R-squared
(pseudoR2 <- ___)
编辑并运行代码