From: N. Brouard <brouard@ined.fr>
Date: Wed, 20 Nov 2002 07:56:21 +0000 (+0000)
Subject: just a sign error in mlikeli
X-Git-Tag: imach-099s7~653
X-Git-Url: https://henry.ined.fr/git/?a=commitdiff_plain;h=964a040c8244e786d0264233ef164f87e0b543c8;p=.git

just a sign error in mlikeli
---

diff --git a/src/imach.c b/src/imach.c
index a18c3c7..2597f59 100644
--- a/src/imach.c
+++ b/src/imach.c
@@ -969,7 +969,7 @@ double func( double *x)
 	s1=s[mw[mi][i]][i];
 	s2=s[mw[mi+1][i]][i];
 	bbh=(double)bh[mi][i]/(double)stepm;
-	lli= (savm[s1][s2]>(double)1.e-8 ?(1.+bbh)*log(out[s1][s2])- bbh*log(savm[s1][s2]):log((1.-bbh)*out[s1][s2]));
+	lli= (savm[s1][s2]>(double)1.e-8 ?(1.+bbh)*log(out[s1][s2])- bbh*log(savm[s1][s2]):log((1.+bbh)*out[s1][s2]));
 	/*lli= (savm[s1][s2]>1.e-8 ?(1.+bbh)*log(out[s1][s2])- bbh*log(savm[s1][s2]):log((1.-bbh)*out[s1][s2]));*/
 	/*lli=(1.+bbh)*log(out[s1][s2])- bbh*log(savm[s1][s2]);*/
 	/*if(lli ==000.0)*/