dnorm和dmvnorm結果不匹配

Question

我注意到從dmvnorm（mvtnorm包）返回到dnorm（baseR）的值之間的實質區別。哪一個是正確的結果，爲什麼？在獨立正態隨機變量的情況下，我應該使用哪一個？

> mvtnorm::dmvnorm(rep(1,2),rep(0,2), diag(2,2)) 
[1] 0.04826618 
> prod(dnorm(rep(1,2),0,2^2)) 
[1] 0.009344515 

Answer 1

對於dmvnorm，指定協方差;對於dnorm，您設置了標準偏差。因此，如果您向前者輸入diag(2,2)，則後者需要sqrt(2)。

prod(dnorm(rep(1,2),0,sqrt(2))) 
# [1] 0.04826618