Skip to content
Snippets Groups Projects
Commit 7cb9854b authored by Houtan Bastani's avatar Houtan Bastani
Browse files

dmm: replace single quotes with double quotes in mL file

parent e6fb6c73
Branches
Tags
No related merge requests found
......@@ -54,7 +54,7 @@ for i=1:nt
error(['Only inv_gamma, normal, and beta are supported for the ' ...
'prior distributions of dmm.']);
end
fprintf(fid, 'pdftheta(%d) = ''%s'' hyptheta(1, %d) = %d %d %d %d\n', i, ...
fprintf(fid, 'pdftheta(%d) = "%s" hyptheta(1, %d) = %d %d %d %d\n', i, ...
shape, i, estimation_info.parameter(i).prior(1).mean, ...
estimation_info.parameter(i).prior(1).stdev, ...
estimation_info.parameter(i).prior(1).interval);
......@@ -63,7 +63,7 @@ fprintf(fid, '&end\n');
%% MCMC
fprintf(fid, '\n&mcmc\n');
fprintf(fid, 'seed=%d thin=%d burnin=%d simulrec=%d hbl=%d MargLik=''%s''', ...
fprintf(fid, 'seed=%d thin=%d burnin=%d simulrec=%d hbl=%d MargLik="%s"', ...
options_.dmm.seed, options_.dmm.thinning_factor, options_.mcmc.drop, ...
options_.mcmc.replic, options_.dmm.block_length, ...
options_.dmm.calc_marg_lik);
......@@ -101,7 +101,7 @@ for i=1:size(options_.dmm.S, 2)
params = estimation_info.transition_probability(tpidx).prior.params;
fprintf(fid, '\n&S%d\n', i);
fprintf(fid, 'dynS%d=''%s'' nS%d=%d hypS%d(1,1)=%s matS%d=''%s''',i,'I',i, ...
fprintf(fid, 'dynS%d="%s" nS%d=%d hypS%d(1,1)=%s matS%d="%s"',i,'I',i, ...
numStates, i, num2str(params), i, options_.dmm.S(i).mat);
fprintf(fid, '\n&end\n');
wroteLatVarFor = [wroteLatVarFor options_.dmm.S(i).mat];
......@@ -110,7 +110,7 @@ end
%% SSM
fprintf(fid, '\n&ssm\n');
fprintf(fid, 'nu=%d nv=%d nx=%d d=%d %d dllname=''%s'' check=''%s''', M_.exo_nbr, ...
fprintf(fid, 'nu=%d nv=%d nx=%d d=%d %d dllname="%s" check="%s"', M_.exo_nbr, ...
length(wroteLatVarFor), options_.dmm.nx, ...
options_.dmm.max_order_of_integration, options_.dmm.num_nonstationary, ...
[cd filesep M_.fname '_dmm.m'], options_.dmm.check_mats);
......@@ -131,7 +131,7 @@ else
Tobs = Tobs - options_.dmm.num_forecasts;
end
fprintf(fid, 'T=%d ny=%d nz=%d nf=%d datasim=''%s'' obs=\n', ...
fprintf(fid, 'T=%d ny=%d nz=%d nf=%d datasim="%s" obs=\n', ...
Tobs, options_.varobs_nbr, M_.exo_det_nbr, options_.dmm.num_forecasts, options_.dmm.simulate_data);
for i=1:size(data,1)
fprintf(fid, '%f ', data.data(i,:));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment