การประมาณค่าพารามิเตอร์สำหรับ Multivariate Skew-Normal
การแจกแจงแบบ Multivariate Normal สามารถประมาณค่าพารามิเตอร์ได้จากค่าเฉลี่ยตัวอย่างและเมทริกซ์ความแปรปรวนร่วมของตัวอย่าง แต่สำหรับการแจกแจงแบบ Skew-Normal จำเป็นต้องใช้กระบวนการวนซ้ำในการประมาณค่า โดยสามารถใช้ฟังก์ชันที่มีอยู่ในแพ็กเกจ sn เพื่อดำเนินการนี้ได้
ในแบบฝึกหัดก่อนหน้า เราได้ตรวจสอบตัวแปร Ht และ Wt จากออบเจกต์ ais.female และพบว่าการแจกแจงแบบ Skew-Normal สามารถอธิบายการกระจายของข้อมูลได้ คราวนี้เราจะมาประมาณค่าพารามิเตอร์ของการแจกแจงนี้กัน
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การแจกแจงความน่าจะเป็นหลายตัวแปรใน R
คำแนะนำการฝึกหัด
- คำนวณพารามิเตอร์ของการแจกแจงแบบ Skew-Normal และกำหนดผลลัพธ์ให้กับออบเจกต์
fit.ais - แสดงเฉพาะพารามิเตอร์ความเบ้ (skewness parameters)
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Fit skew-normal parameters
fit.ais <- ___(y = cbind(ais.female$Ht, ais.female$Wt), opt.method = ___)
# Print the skewness parameters
fit.ais$dp$___