穀物產量建模

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 函式撰寫入門