library(mvtnorm)
dmvnorm(x, mean = rep(0, p), sigma = diag(p), log = FALSE)dmvnorm为多元正态分布提供了密度函数。第一个参数x到底代表什么?文档中写着“矢量或分位数矩阵,如果x是矩阵,则每一行都被视为分位数。”
> dmvnorm(x=c(0,0), mean=c(1,1))
[1] 0.0585下面是帮助页面上的示例代码。在这种情况下,您正在生成以平均值1和sd 1为正态分布的分位数0的概率(假设这是默认的)。既然这是一个多元正态密度函数,并且传递了分位数(0, 0)的向量,为什么输出不是概率的向量呢?
发布于 2017-06-05 13:28:12
仅以二元正常(X1, X2)为例,通过传入x = (0, 0),您可以得到P(X1 = 0, X2 = 0),它是一个值。为什么你希望得到一个向量?
如果你想要一个向量,你需要传递一个矩阵。例如,x = cbind(c(0,1), c(0,1))给出了
P(X1 = 0, X2 = 0)
P(X1 = 1, X2 = 1)在这种情况下,矩阵的每一行都被并行处理。
https://stackoverflow.com/questions/44369959
复制相似问题