Bug fix

Fixed the beta distribution in updatehyperpars.
WTCN-computational-anatomy-group · Sep 21, 2020 · 6293a9f · 6293a9f
1 parent 0c81787
commit 6293a9f
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/spm_gmm_lib.m b/spm_gmm_lib.m
@@ -1502,12 +1502,12 @@
     alpha0 = b0_priors{1};
     beta0  = b0_priors{2};
     b0(k)  = 0;
-    alph   = N*S;
-    bet    = 1;
+    alph   = 0.5*N*S+1;
+    bet    = 0;
     for s=1:S
         [m,b,V,n] = get_posteriors(cluster,s);
         m1  = m(:,k) - m0(:,k);
-        bet = bet + m1.' * (n(k)*V(:,:,k)) * m1 + N/b(k);
+        bet = bet + 0.5*(m1.' * (n(k)*V(:,:,k)) * m1 + N/b(k));
     end
     b0(k) = (alph+alpha0-1)/(bet+beta0);
     % ---------------------------------------------------------------------