diff --git a/hf.for b/hf.for
index 1771f3d9ee24ad2fa7db389f6103eb0072b1ec3f..6220c1a960f877e56fc503f5e53a16863eaafd1c 100644
--- a/hf.for
+++ b/hf.for
@@ -58,7 +58,7 @@ C Initial condition (stationary MS-VAR(1))
       CALL ergodic(nk,P,PE)
       inx = IYK(1,nx+1)
       DO J = 1,inx
-	 mu(J) = a(IYK(1,J),IS(4))+
+	 mu(J) = a(IYK(1,J),IS(4))
      #       + SUM(F(IYK(1,J),IYK(1,1:inx),IS(5))
      #       * yk(1,IYK(1,1:inx)))
      #       + SUM(c(IYK(1,J),1:nz,IS(1))*yk(1,nx+1:nx+nz))
@@ -92,7 +92,7 @@ C ----------------------------------------------------------
 c y(t) = c(t)z(t) + x(t)
 C x(t) = a(t)     + F(t)x(t-1) + R(t)u(t)
         DO J = 1,inx
-	   mu(J) = a(IYK(inobs,J),IS(4))+
+	   mu(J) = a(IYK(inobs,J),IS(4))
      #         + SUM(F(IYK(inobs,J),IYK(inobs,1:inx),IS(5))
      #         * yk(inobs-1,IYK(inobs-1,1:inx)))
      #         + SUM(c(IYK(inobs,J),1:nz,IS(1))*yk(inobs,nx+1:nx+nz))
@@ -149,7 +149,7 @@ C ------------------------------------------------------------------------------
 	  DO I=1,nk
          CALL int2seq(I,nv,INFOS,SEQ,IS)
          DO K = 1,inx
-	    mu(K) = a(IYK(inobs,K),IS(4))+
+	    mu(K) = a(IYK(inobs,K),IS(4))
      #          + SUM(F(IYK(inobs,K),IYK(inobs,1:inx),IS(5))
      #          * yk(inobs-1,IYK(inobs-1,1:inx)))
      #          + SUM(c(IYK(inobs,K),1:nz,IS(1))*yk(inobs,nx+1:nx+nz))
diff --git a/mengwong.for b/mengwong.for
index 2ff6c442a6e943a952fd737e5b1e0dce7f3a6030..a7ac08a03f511e98890ebf99ac13df15d6ae61f3 100644
--- a/mengwong.for
+++ b/mengwong.for
@@ -419,7 +419,7 @@ C PRIOR for PSI and Q(PSI)~Dirichlet(a1,a2,...,aN)
       MLMW(1,1)   = DLOG(MLMW(2,1)) + DET
       MLMW(1:2,2) = SS(1,1)*G/SUM(MAT(:,1))**2 +
      +              SS(2,2)*G/SUM(MAT(:,2))**2 +
-     +     	    - 2.D0*SS(1,2)*G/(SUM(MAT(:,1))*SUM(MAT(:,2)))
+     +     	    (- 2.D0)*SS(1,2)*G/(SUM(MAT(:,1))*SUM(MAT(:,2)))
 
       MLMW(2,1) = MLMW(1,1) ! log scale
       DO 500 I=1,10
diff --git a/mengwong2.for b/mengwong2.for
index 0c121c583d10ac32f89fbafcacbc21e4adf0b631..2778e1dd357bbc07505a597e5c6b6811b65e9ac9 100644
--- a/mengwong2.for
+++ b/mengwong2.for
@@ -418,7 +418,7 @@ C PRIOR for PSI and Q(PSI)~Dirichlet(a1,a2,...,aN)
       MLMW(1,1)   = DLOG(MLMW(2,1)) + DET
       MLMW(1:2,2) = SS(1,1)*G/SUM(MAT(:,1))**2 +
      +              SS(2,2)*G/SUM(MAT(:,2))**2 +
-     +     	    - 2.D0*SS(1,2)*G/(SUM(MAT(:,1))*SUM(MAT(:,2)))
+     +     	    (- 2.D0)*SS(1,2)*G/(SUM(MAT(:,1))*SUM(MAT(:,2)))
 
       MLMW(2,1) = MLMW(1,1) ! log scale
       DO 500 I=1,10