From f1254ca2d1fbc7d0ae024da5efcf98007c11de3a Mon Sep 17 00:00:00 2001
From: Christian Roever <christian.roever@med.uni-goettingen.de>
Date: Fri, 6 Aug 2021 08:50:21 +0200
Subject: [PATCH] small bug fix
---
R/bmr.R | 12 ++++++++----
1 file changed, 8 insertions(+), 4 deletions(-)
diff --git a/R/bmr.R b/R/bmr.R
index 3d84496..5b8d819 100644
--- a/R/bmr.R
+++ b/R/bmr.R
@@ -386,10 +386,14 @@ bmr.default <- function(y, sigma, labels = names(y),
}
logdens <- log(sum(exp(log(support$weight) + logdens)))
} else { # joint density of (tau,beta)
- cm <- conditionalmoments(tau=tau)
- logdens <- (taulogpostdens(tau)
- + mvtnorm::dmvnorm(x=beta[idx], mean=cm$mean[idx],
- sigma=cm$covariance[idx,idx,drop=FALSE], log=TRUE))
+ if (tau>=0) {
+ cm <- conditionalmoments(tau=tau)
+ logdens <- (taulogpostdens(tau)
+ + mvtnorm::dmvnorm(x=beta[idx], mean=cm$mean[idx],
+ sigma=cm$covariance[idx,idx,drop=FALSE], log=TRUE))
+ } else {
+ logdens <- -Inf
+ }
}
return(logdens)
}
--
GitLab