fix fixed effect sign for longitudinal likelihood

premium-profile-regression · Mar 22, 2019 · 919e744 · 919e744
1 parent 6d23d2a
commit 919e744
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/PReMiuM/src/include/PReMiuMModel.h b/PReMiuM/src/include/PReMiuMModel.h
@@ -2535,7 +2535,7 @@ double logPYiGivenZiWiLongitudinal(const pReMiuMParams& params, const pReMiuMDat
 				timesk[counter+j] = times[tStart[i]-1+j];
 				meanVec(counter+j) = 0.0;
 				for(unsigned int b=0;b<nFixedEffects;b++){
-					yk(counter+j)+=params.beta(b,0)*dataset.W(i,b);
+					yk(counter+j) -= params.beta(b,0)*dataset.W(i,b);
 				}
 			}
 			//RJ tidy up vector copying