quick hotfix gaussian deviance

stan-dev · Mar 8, 2021 · 0cabd73 · 0cabd73
1 parent c70dacb
commit 0cabd73
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/extend_family.R b/R/extend_family.R
@@ -136,7 +136,7 @@ extend_family_gaussian <- function(family) {
     if (NCOL(y) < NCOL(mu)) {
       y <- matrix(y, nrow = length(y), ncol = NCOL(mu))
     }
-    -2 * weights * (-0.5 / dis * (y - mu)^2 - log(dis))
+    -2 * weights * (-0.5 / dis^2 * (y - mu)^2 - log(dis))
   }
   ppd_gauss <- function(mu, dis, weights = 1) rnorm(length(mu), mu, dis)