Changeset 7a0e5ad8 in tol

Timestamp:

Sep 3, 2009, 10:44:56 PM (17 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

1941-docker-image-to-test, devel, master

Children:

6b49d63

Parents:

e7e91d9

Message:

Fixed bug in draw.mh

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@1324 be4846ae-f051-483a-831f-b64e7d08527f

File:

• tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_stdLinearBlock.tol (modified) (1 diff)

tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_stdLinearBlock.tol

@@ -471 +471 @@
       Real a = Rand(0,1);
       Real c = Exp((p1+q0)-(p0+q1));
-      Real s = Sqrt(gsl_ran_chisq(m)/m);
+      Real s = s0*Sqrt(gsl_ran_chisq(m)/m);
       Real accept := Case(
-        s1 <= 1, True,
+      //s1 <= 1, True,
         s1 <= s, True,
         1  == 1, a<c);