Density and log-density for the multivariate normal distribution
with mean equal to mu
and variance matrix equal to sigma
.
dmvnorm(y, mu=NULL, sigma=NULL)
logdmvnorm(y, mu=NULL, sigma=NULL)
dmvnorm
gives the densities, while
logdmvnorm
gives the logarithm of the densities.
Either a
This
This code is written to be efficient, using the qr-decomposition of the
covariance matrix (and using it only once, rather than recalculating it
for both the determinant and the inverse of sigma
).