运行一个模型
您在按年份绘制产量时看到的平滑趋势线,是通过广义加性模型(GAM)来确定曲线位置的。这类模型非常适合拟合非线性曲线。为了能对未来产量做出预测,我们来显式运行这个模型。运行该 GAM 的语法如下所示。
gam(response ~ s(explanatory_var1) + explanatory_var2, data = dataset)
这里,s() 表示"对变量进行平滑",平滑非常粗略地说就是非线性。
已加载 mgcv 和 dplyr;corn 与 wheat 数据集可用,且已扩展包含公制单位与人口普查地区信息。
本练习是课程的一部分
R 函数编写入门
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Run a generalized additive model of yield vs. smoothed year and census region
___(___ ~ s(___) + ___, data = ___)