计算多元正态分布的密度
在许多统计任务中,例如假设检验、聚类和似然计算,您需要计算给定多元正态分布的密度。在本练习中,您将使用 dmvnorm() 函数,针对之前生成的样本 multnorm.sample 中的每个观测值,计算在指定均值与方差-协方差矩阵下的多元正态密度。
均值和方差-协方差矩阵已为您预加载为对象 mu.sim 和 sigma.sim。
本练习是课程的一部分
R 中的多元概率分布
练习说明
- 使用
dmvnorm()计算multnorm.sample中 100 个样本在二维正态分布下的密度高度。 - 使用
scatterplot3d()绘制 3D 散点图,展示每个生成样本点对应的密度高度。
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Calculate density
multnorm.dens <- dmvnorm(multnorm.sample, mean = ___, sigma = ___)
# Create scatter plot of density heights
___(cbind(___),
color="blue", pch="", type = "h",
xlab = "x", ylab = "y", zlab = "density")