diff --git a/matlab/ms-sbvar/cstz/fn_rnrprior_covres_dobs.m b/matlab/ms-sbvar/cstz/fn_rnrprior_covres_dobs.m
index 394d27b90ffe3faa779a82fd9eba06c86b7fa8f0..41a22cab13132fd375817ed06fe8d335b8b58834 100644
--- a/matlab/ms-sbvar/cstz/fn_rnrprior_covres_dobs.m
+++ b/matlab/ms-sbvar/cstz/fn_rnrprior_covres_dobs.m
@@ -141,6 +141,8 @@ for i = 1:lags
          sgpbid((i-1)*nvar+j) = lagdecay^2/sgsh(j);  % ith equation
       elseif (q_m==4)
          sgpbid((i-1)*nvar+j) = (1/i^mu(4))^2/sgsh(j);  % ith equation
+      elseif (q_m==1)
+         sgpbid((i-1)*nvar+j) = (1/(i*4)^mu(4))^2/sgsh(j);  % ith equation
       else
 			error('Incompatibility with lags, check the possible errors!!!')
          %warning('Incompatibility with lags, check the possible errors!!!')
diff --git a/matlab/ms-sbvar/cstz/fn_rnrprior_covres_dobs_tv2.m b/matlab/ms-sbvar/cstz/fn_rnrprior_covres_dobs_tv2.m
index e3d10d4094baedd3fb25aafb0e43d6b43fa4d718..bc4fec59764aeac7eb5b383859143a6c7cf2d6f7 100644
--- a/matlab/ms-sbvar/cstz/fn_rnrprior_covres_dobs_tv2.m
+++ b/matlab/ms-sbvar/cstz/fn_rnrprior_covres_dobs_tv2.m
@@ -163,6 +163,8 @@ for i = 1:lags
          sgpbid((i-1)*nvar+j) = lagdecay^2/sgsh(j);  % ith equation
       elseif (q_m==4)
          sgpbid((i-1)*nvar+j) = (1/i^mu(4))^2/sgsh(j);  % ith equation
+      elseif (q_m==1)
+         sgpbid((i-1)*nvar+j) = (1/(i*4)^mu(4))^2/sgsh(j);  % ith equation
       else
 			error('Incompatibility with lags, check the possible errors!!!')
          %warning('Incompatibility with lags, check the possible errors!!!')
diff --git a/matlab/ms-sbvar/msstart2.m b/matlab/ms-sbvar/msstart2.m
index 14bb315546647aa5c6f16c373e9836c8646c878b..e9dd08f6687c6ff655e5732e5b2907e5929a7de1 100644
--- a/matlab/ms-sbvar/msstart2.m
+++ b/matlab/ms-sbvar/msstart2.m
@@ -176,8 +176,11 @@ if (options_.ms.freq==12)
 elseif (options_.ms.freq==4)
    nStart=(yrStart-options_.ms.initial_year )*4+qmStart-options_.ms.initial_subperiod ;  % positive number of months at the start
    nEnd=(yrEnd-options_.ms.final_year )*4+qmEnd-options_.ms.final_subperiod ;     % negative number of months towards end
+elseif (options_.ms.freq==1)
+   nStart=(yrStart-options_.ms.initial_year )*1+qmStart-options_.ms.initial_subperiod ;  % positive number of months at the start
+   nEnd=(yrEnd-options_.ms.final_year )*1+qmEnd-options_.ms.final_subperiod ;     % negative number of months towards end
 else
-   error('Error: this code is only good for monthly/quarterly data!!!')
+   error('Error: this code is only good for monthly/quarterly/yearly data!!!')
    return
 end
 %