diff --git a/matlab/prior_draw.m b/matlab/prior_draw.m index 8cebcdfe479ab59dbdb08dc2806808820362237e..36bdc80d4fe159ae5136807e054174302a30718c 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