diff --git a/matlab/prior_posterior_statistics.m b/matlab/prior_posterior_statistics.m
index b388007d46363c19af6b3b53ce98cb1e61ceacc1..f06e0d433ea024755cd1d0069a03c935d2269685 100644
--- a/matlab/prior_posterior_statistics.m
+++ b/matlab/prior_posterior_statistics.m
@@ -106,7 +106,7 @@ if horizon
     MAX_nforc2 = min(B,ceil(MaxNumberOfBytes/((endo_nbr)*(horizon+maxlag))/ ...
                             8));
     IdObs    = bayestopt_.mfys;
-    
+
 end
 MAX_momentsno = min(B,ceil(MaxNumberOfBytes/(get_moments_size(options_)*8)));
 
@@ -232,7 +232,7 @@ else
                         'options_', options_, ...
                         'bayestopt_', bayestopt_, ...
                         'estim_params_', estim_params_, ...
-                        'oo_', oo_);    
+                        'oo_', oo_);
     % which files have to be copied to run remotely
     NamFileInput(1,:) = {'',[M_.fname '_static.m']};
     NamFileInput(2,:) = {'',[M_.fname '_dynamic.m']};
@@ -240,11 +240,10 @@ else
         NamFileInput(length(NamFileInput)+1,:)={'',[M_.fname '_steadystate.m']};
     end
     [fout] = masterParallel(options_.parallel, 1, B,NamFileInput,'prior_posterior_statistics_core', localVars,globalVars, options_.parallel_info);
-    
+
 end
 ifil = fout(end).ifil;
 
-% ??????????
 stock_gend=gend;
 stock_data=Y;
 save([DirectoryName '/' M_.fname '_data.mat'],'stock_gend','stock_data');
@@ -301,4 +300,4 @@ if ~isnumeric(options_.parallel),
     if leaveSlaveOpen == 0,
         closeSlave(options_.parallel,options_.parallel_info.RemoteTmpFolder),
     end
-end
+end
\ No newline at end of file
diff --git a/matlab/prior_posterior_statistics_core.m b/matlab/prior_posterior_statistics_core.m
index 41b5b6d33f579d4d8b342999e162083e371b2076..63443a6ae12bb054cd793df6e51bc64a110ac97e 100644
--- a/matlab/prior_posterior_statistics_core.m
+++ b/matlab/prior_posterior_statistics_core.m
@@ -133,14 +133,14 @@ end
 if run_smoother
   stock_smooth=NaN(endo_nbr,gend,MAX_nsmoo);
   stock_update=NaN(endo_nbr,gend,MAX_nsmoo);
-  stock_innov=NaN(M_.exo_nbr,gend,MAX_ninno);  
+  stock_innov=NaN(M_.exo_nbr,gend,MAX_ninno);
   if horizon
       stock_forcst_mean= NaN(endo_nbr,horizon+maxlag,MAX_nforc1);
       stock_forcst_point = NaN(endo_nbr,horizon+maxlag,MAX_nforc2);
   end
 end
 if nvn
-  stock_error = NaN(size(varobs,1),gend,MAX_nerro);
+  stock_error = NaN(length(varobs),gend,MAX_nerro);
 end
 if naK
     stock_filter_step_ahead =NaN(length(options_.filter_step_ahead),endo_nbr,gend+max(options_.filter_step_ahead),MAX_naK);
@@ -150,10 +150,6 @@ stock_logpo = NaN(MAX_nruns,1);
 stock_ys = NaN(MAX_nruns,endo_nbr);
 
 for b=fpar:B
-
-    %    [deep, logpo] = GetOneDraw(type);
-    %    set_all_parameters(deep);
-    %    dr = resol(oo_.steady_state,0);
     if strcmpi(type,'prior')
 
         [deep, logpo] = GetOneDraw(type);
@@ -304,22 +300,7 @@ for b=fpar:B
         end
         irun(7) = 1;
     end
-
-    % if moments_varendo && (irun(8) > MAX_momentsno || b == B)
-    %    stock = stock_moments(1:irun(8)-1);
-    %    ifil(8) = ifil(8) + 1;
-    %    save([DirectoryName '/' M_.fname '_moments' int2str(ifil(8)) '.mat'],'stock');
-    %    if RemoteFlag==1,
-    %    OutputFileName_moments = [OutputFileName_moments; {[DirectoryName filesep], [M_.fname '_moments' int2str(ifil(8)) '.mat']}];
-    %    end
-    %    irun(8) = 1;
-    % end
-
-    %   DirectoryName=TempPath;
-
-
     dyn_waitbar((b-fpar+1)/(B-fpar+1),h);
-
 end
 
 myoutput.ifil=ifil;
@@ -332,7 +313,6 @@ if RemoteFlag==1,
                         OutputFileName_param;
                         OutputFileName_forc_mean;
                         OutputFileName_forc_point];
-    % OutputFileName_moments];
 end
 
 dyn_waitbar_close(h);