that's the bvnpdf() code:
bvnpdf <- function(xy, mu, Siginv, detSig){
f<-numeric()
x <- xy[[1]]
if(xy[[1]]>xy[[2]] & (1==2)) {
f <- 0} else {
v <- rbind(xy[1], xy[2])
e <- t(v-mu) %*% Siginv %*% (v-mu)
f <- as.numeric(exp(-e/2)/(2*pi)/sqrt(detSig))
}
f
}
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
