Changeset 664d021 in tol

Timestamp:

Nov 20, 2008, 5:04:40 PM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

f9ee626

Parents:

af055cb

Message:

At first iteration of Logit is B=0 and so p=X*B=0

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

File:

• tol/bmath/bstat/logit.cpp (modified) (1 diff)

tol/bmath/bstat/logit.cpp

@@ -74 +74 @@
   clock_t tm_1, tm_0 = BTimer::Clocks();
   BReal tm_dif;
+
   for(iter=0; iter<BDat::MaxIter(); iter++)
   {
     H_.SetAllValuesTo(0);
     G.SetAllValuesTo(0);
-    p = X*B;
+    if(iter==0)
+    {
+      p.SetAllValuesTo(0);
+    }
+    else
+    {
+    //Std(BText("\nTRACE Logit B=\n")+B.Name());
+    //Std(BText("\nTRACE Logit X=\n")+X.Name());
+      p = X*B;
+    //Std(BText("\nTRACE Logit p=\n")+p.Name());
+    }
     oldLlh = lnLikelyhood;
     lnLikelyhood=0;