R 多元正态分布

我试图通过使用库

MASS

和函数

dmvnorm

在R中使用多元正态分布。 我有向量：

Y = c(26.385112, 17.108580, 11.907650,  4.737202)
Mu = c(31.19789, 30.33983, 29.48177, 28.62371)

和协方差矩阵：

> Sigma
           [,1]      [,2]      [,3]      [,4]
[1,]  10.861298 -3.910141  -9.40906 -14.90798
[2,]  -3.910141 24.890056  35.14521  54.67289
[3,]  -9.409060 35.145213  88.97201 124.25376
[4,] -14.907979 54.672890 124.25376 203.10715

> dput(Sigma)
structure(c(10.8612978825767, -3.91014077427266, -9.4090596651869, 
-14.9079785561011, -3.91014077427266, 24.8900558453267, 35.1452129330559, 
54.6728897867201, -9.4090596651869, 35.1452129330559, 88.9720052972337, 
124.253758129541, -14.9079785561011, 54.6728897867201, 124.253758129541, 
203.107146238298), .Dim = c(4L, 4L))

所以，当我应用

dmvnorm

时，我只得到一个密度，它假设每个值有4个密度

dmvnorm(Y,mean=Mu,Sigma)
[1] 2.756653e-07

---

还是我错了？如果没有，你能帮我吗？

我在包

MASS

中没有看到函数

dmvnorm

。哦，对不起！！你是对的。这是一个库（mvtnorm），您能

dput（Sigma）

并将结果发布到您的问题中吗？我在dmvnorm（Y，mean=Mu，Sigma）中得到了

错误：Sigma必须是一个对称矩阵，与您发布的矩阵一致。四元向量表示来自四维MVN分布的单个观测值，具有单个概率密度值。