การรันโมเดล

เส้นแนวโน้มที่เรียบในกราฟผลผลิตตามช่วงเวลานั้นใช้ Generalized Additive Model (GAM) เพื่อกำหนดตำแหน่งของเส้น โมเดลประเภทนี้เหมาะสำหรับการ fit เส้นโค้งแบบ nonlinear เพื่อให้สามารถพยากรณ์ผลผลิตในอนาคตได้ มาลองรันโมเดลกันโดยตรง ไวยากรณ์สำหรับรัน GAM มีรูปแบบดังนี้

gam(response ~ s(explanatory_var1) + explanatory_var2, data = dataset)

โดยที่ s() หมายถึง "ทำให้ตัวแปรนั้น smooth" ซึ่ง smooth หมายถึงการเป็นแบบ nonlinear โดยประมาณ

โหลด mgcv และ dplyr ไว้แล้ว ชุดข้อมูล corn และ wheat พร้อมใช้งาน และได้เพิ่มคอลัมน์หน่วยเมตริกและข้อมูลภูมิภาคสำมะโนประชากรเรียบร้อยแล้ว