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

การคำนวณ dmvnorm บนกริด

ในการแสดงภาพพื้นผิวความหนาแน่นของการแจกแจงปกติแบบสองตัวแปร จำเป็นต้องคำนวณความหนาแน่นบนกริดที่มีพิกัด x และ y ละเอียด จากนั้นใช้ฟังก์ชันพล็อตพื้นผิว 3 มิติ เช่น persp() เพื่อแสดงผล

กริดขนาด 40 x 40 ที่สร้างจากคำสั่ง mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) ได้ถูกโหลดไว้ให้แล้ว

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

การแจกแจงความน่าจะเป็นหลายตัวแปรใน R

ดูคอร์ส

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

  • คำนวณความหนาแน่นของการแจกแจงปกติที่มีค่าเฉลี่ย mu.sim และเมทริกซ์ความแปรปรวนร่วม sigma.sim บนค่ากริดในออบเจกต์ mvals
  • ใช้ฟังก์ชันพล็อต persp() เพื่อแสดงพื้นผิวความหนาแน่นบนกริด

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

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

# Calculate density over the specified grid
mvds <- ___
matrix_mvds <-  matrix(___, nrow = 40)

# Create a perspective plot
persp(___, theta = 80, phi = 30, expand = 0.6, shade = 0.2, col = "lightblue", xlab = "x", ylab = "y", zlab = "dens")
แก้ไขและรันโค้ด