From 2400b6a0e799a164a9b99e392c950d77f16df81f Mon Sep 17 00:00:00 2001 From: draguar <33092622+draguar@users.noreply.github.com> Date: Wed, 31 Mar 2021 14:59:05 +0200 Subject: [PATCH] Corrected omega value update, as previously done for alpha See issue #2. --- ChrMod_main.m | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/ChrMod_main.m b/ChrMod_main.m index e0f06cb..6bd9ab1 100644 --- a/ChrMod_main.m +++ b/ChrMod_main.m @@ -99,7 +99,7 @@ function ChrMod_main(ChrNo,n) %Grid search for omega lnL_om=GS_omega(ChrIF,weight,Ensemble,EstLen,Dim); omega=find(lnL_om==max(lnL_om)); - omega=0.25+0.25*omega; + omega=0.25+0.25*mean(omega); OmegaRec(k)=omega; end