เริ่มต้นใช้งานเริ่มต้นใช้งานได้ฟรี

วิธี Percentile - Bootstrap CI สำหรับค่าความชัน

อีกวิธีหนึ่งในการสร้าง CI สำหรับค่าความชัน คือการใช้ค่า percentile จากการกระจายตัวของสถิติความชันที่ได้จาก bootstrap โดยทั่วไป CI จะถูกสร้างขึ้นในลักษณะที่อัตราครอบคลุม (coverage rate) ตลอดช่วงการวิเคราะห์จะเท่ากับ (1-alpha)*100% หาก alpha = 0.05 ช่วงความเชื่อมั่น 95% จะครอบคลุมพารามิเตอร์ที่สนใจ 95% ของเวลาทั้งหมด ส่วนอีก 5% ที่ช่วงนั้นพลาด บางครั้งช่วงจะสูงเกินไป (2.5%) และบางครั้งต่ำเกินไป (2.5%)

ค่าประมาณ bootstrap ของ slope ที่ชื่อว่า boot_slope ได้ถูกโหลดไว้ใน workspace แล้ว

แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร

การอนุมานสำหรับการถดถอยเชิงเส้นใน R

ดูคอร์ส

คำแนะนำการฝึกหัด

  • กำหนดให้ alpha เท่ากับ 0.05 (สำหรับงานของคุณเอง สามารถเลือกระดับความเชื่อมั่นที่ต่างออกไปได้)
  • คำนวณค่า percentile ที่จำเป็นสำหรับการสร้างช่วงความเชื่อมั่น
    • ค่าตัดต่ำสุดอยู่ที่ครึ่งหนึ่งของ alpha
    • ค่าตัดสูงสุดอยู่ที่หนึ่งลบด้วยครึ่งหนึ่งของ alpha
  • สร้างช่วงความเชื่อมั่นของ stat โดยใช้ quantile() และค่า percentile ที่คำนวณไว้ โดยปลายทั้งสองด้านของช่วงต้องตั้งชื่อว่า lower และ upper

แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ

ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์

# Set alpha
alpha <- ___

# Set the lower percentile cutoff
p_lower <- ___

# Set the upper percentile cutoff
p_upper <- ___

# Create a confidence interval of stat using quantiles
boot_slope %>% 
  ___
แก้ไขและรันโค้ด