From 08375e7ba6624ec2f8aeb8a49e869ddc61e56b1b Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?St=C3=A9phane=20Adjemian=20=28Lupi=29?=
 <stephane.adjemian@univ-lemans.fr>
Date: Fri, 4 Mar 2016 16:45:06 +0100
Subject: [PATCH] Fixed wrong ordering of the hyperparameters when walling the
 routine for drawing random deviates from the Weibull distribution.

---
 matlab/prior_draw.m | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/matlab/prior_draw.m b/matlab/prior_draw.m
index 8cebcdfe4..36bdc80d4 100644
--- a/matlab/prior_draw.m
+++ b/matlab/prior_draw.m
@@ -166,10 +166,10 @@ if inverse_gamma_2_draws
 end
 
 if weibull_draws
-    pdraw(weibull_index) = wblrnd(p6(weibull_index), p7(weibull_index)) + p3(weibull_index);
+    pdraw(weibull_index) = wblrnd(p7(weibull_index), p6(weibull_index)) + p3(weibull_index);
     out_of_bound = find( (pdraw(weibull_index)'>ub(weibull_index)) | (pdraw(weibull_index)'<lb(weibull_index)));
     while ~isempty(out_of_bound),
-        pdraw(weibull_index(out_of_bound)) = wblrnd(p6(weibull_index(out_of_bound)),p7(weibull_index(out_of_bound)))+p3(weibull_index(out_of_bound));
+        pdraw(weibull_index(out_of_bound)) = wblrnd(p7(weibull_index(out_of_bound)),p6(weibull_index(out_of_bound)))+p3(weibull_index(out_of_bound));
         out_of_bound = find( (pdraw(weibull_index)'>ub(weibull_index)) | (pdraw(weibull_index)'<lb(weibull_index)));
     end
 end