การใช้ฟังก์ชัน pmvnorm
นอกจากการคำนวณความหนาแน่นของการแจกแจงปกติหลายตัวแปรแล้ว บางครั้งยังต้องคำนวณการแจกแจงสะสม (cumulative distribution) เพื่อหาปริมาตรใต้เส้นโค้งความหนาแน่นระหว่างค่าสองค่าที่กำหนด ในแบบฝึกหัดนี้จะใช้ฟังก์ชัน pmvnorm() เพื่อคำนวณการแจกแจงสะสมของการแจกแจงปกติสองตัวแปร (bivariate normal) ที่กำหนดให้
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การแจกแจงความน่าจะเป็นหลายตัวแปรใน R
คำแนะนำการฝึกหัด
- คำนวณปริมาตรใต้การแจกแจงปกติสองตัวแปรมาตรฐาน ระหว่าง \(\begin{pmatrix} -1\\ -1 \end{pmatrix}\) และ \(\begin{pmatrix} 1 \\ 1 \end{pmatrix}\)
- คำนวณปริมาตรระหว่าง \(\begin{pmatrix} -5 \\ -5\end{pmatrix}\) และ \(\begin{pmatrix} 5 \\ 5 \end{pmatrix}\) สำหรับการแจกแจงปกติสองตัวแปรที่มีค่าเฉลี่ย
mu.simและความแปรปรวนsigma.sim
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Volume under a bivariate standard normal
pmvnorm(lower = ___, upper = ___)
# Volume under specified mean and variance-covariance matrix
pmvnorm(___, mean = ___, sigma = ___)