การพล็อตโมเดล Bayesian
ในแบบฝึกหัดก่อนหน้า เราได้ประมาณค่าโมเดล Bayesian เพื่อพยากรณ์ความนิยมของเพลง (popularity) จากอายุของเพลง (song_age) มาลองสร้างภาพแสดงผลโมเดลกัน โดยใช้ชุดข้อมูล songs และออบเจ็กต์ stan_model ที่โหลดไว้แล้ว สร้างกราฟแสดงข้อมูลพร้อมเส้นถดถอยที่ประมาณได้ด้วย ggplot2
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การสร้างโมเดลการถดถอยแบบ Bayesian ด้วย rstanarm
คำแนะนำการฝึกหัด
- บันทึกสรุปพารามิเตอร์ของโมเดลในรูปแบบ tidy ลงใน
tidy_coef - ดึงค่าประมาณ intercept และ slope จาก
tidy_coef - สร้างกราฟแสดงข้อมูลและเส้นถดถอยที่ประมาณได้ โดยให้
song_ageอยู่บนแกน x และpopularityอยู่บนแกน y
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Save the model parameters
tidy_coef <- ___(stan_model)
# Extract intercept and slope
model_intercept <- tidy_coef$___[1]
model_slope <- tidy_coef$___[2]
# Create the plot
ggplot(songs, aes(x = ___, y = ___)) +
geom_point() +
geom_abline(intercept = ___, slope = ___)