ข้อมูลอุบัติเหตุรถจักรยานยนต์: แนวทางเชิงเส้น
ในแบบฝึกหัดแรกนี้ จะได้ fit โมเดลเชิงเส้นกับชุดข้อมูลและแสดงผลเป็นภาพ เพื่อดูว่าโมเดลจับความสัมพันธ์ในข้อมูลได้ดีเพียงใด ชุดข้อมูลนี้จัดเก็บอยู่ใน data frame ชื่อ mcycle ซึ่งบันทึกค่าความเร่งของหุ่นทดสอบการชนในอุบัติเหตุรถจักรยานยนต์ โดยมีค่าความเร่ง (g) อยู่ในคอลัมน์ accel และเวลา (มิลลิวินาที) อยู่ในคอลัมน์ times
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การสร้างโมเดลแบบไม่เชิงเส้นด้วย Generalized Additive Models (GAMs) ใน R
คำแนะนำการฝึกหัด
- ใช้ฟังก์ชัน
head()และplot()เพื่อดูข้อมูลใน data frame ชื่อmcycle - ใช้ฟังก์ชัน
lm()เพื่อ fit โมเดลกับข้อมูลmcycleโดยให้ตัวแปรaccelเป็นฟังก์ชัน เชิงเส้น ของtimes - แสดงผลการ fit โมเดลโดยใช้คำสั่ง
termplot()ที่เตรียมไว้ให้
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Examine the mcycle data frame
head(___)
plot(___)
# Fit a linear model
lm_mod <- lm(___, data = mcycle)
# Visualize the model
termplot(lm_mod, partial.resid = TRUE, se = TRUE)