建模粮食产量

R 函数编写入门

Richie Cotton

Data Evangelist at DataCamp

这是您在上个练习中绘制的玉米产量分面时间序列。共有九个面板，分别代表各普查地区。每个面板中，产量在约 1950 年前缓慢上升，之后快速上升。

线性模型 vs. 广义可加模型

线性模型


lm(
  response_var ~ explanatory_var1 + explanatory_var2, 
  data = dataset
)

广义可加模型

library(mgcv)
gam(
  response_var ~ s(explanatory_var1) + explanatory_var2, 
  data = dataset
)

predict_this <- data.frame(
  explanatory_var1 = c("some", "values"),
  explanatory_var2 = c("more", "values")
)

predicted_responses <- predict(model, predict_this, type = "response")

predict_this %>%
  mutate(predicted_responses = predicted_responses)

R 函数编写入门