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)
的向量,为什么输出不是概率向量?
答案 0 :(得分:1)
以双变量正常(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)
在这种情况下,矩阵的每一行都是并行处理的。