Changeset 700 for applications/bdmtoolbox/tutorial/userguide/frg_estim.m

Timestamp:

11/04/09 22:54:58 (15 years ago)

Author:

smidl

Message:

Making tutorial/userguide example work again (changes of mpdf and bayes)

Files:

• applications/bdmtoolbox/tutorial/userguide/frg_estim.m (modified) (3 diffs)

applications/bdmtoolbox/tutorial/userguide/frg_estim.m

@@ -17 +17 @@
 %%%%%% Random walk on frg - Dirichlet 
 phi_pdf.class = 'mDirich';         % random walk on coefficient phi
-phi_pdf.rv    = RV('phi',2);       % 2D random walk - frg is the first element
+phi_pdf.rv    = RV({'phi','1_phi'});       % 2D random walk - frg is the first element
 phi_pdf.k     = 0.01;              % width of the random walk
 phi_pdf.betac = [0.01 0.01];       % stabilizing elememnt of random walk
@@ -26 +26 @@
 E.parameter_pdf = phi_pdf;         % Random walk is the parameter evolution model
 E.res_threshold = 1.0;             % resampling parameter
-E.n = 20;                          % number of particles
+E.n = 10;                          % number of particles
 E.prior.class = 'eDirich';         % prior on non-linear part
 E.prior.beta  = [2 1]; % 
@@ -36 +36 @@
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % plot results
-ndat = size(M.u,1);
+ndat = size(M.DS_u,1);
 
 figure(1);
 subplot(2,2,1);
-plotestimates(true_theta, M.MPFmean_theta, M.MPFlb_theta, M.MPFub_theta);
+plotestimates(true_theta, M.MPF_apost_mean_theta, M.MPF_apost_lb_theta, M.MPF_apost_ub_theta);
 title(' Regression parameters \theta')
 set(gca,'YLim',[-1.5,1]);
 
 subplot(2,2,2);
-plotestimates(true_R, M.MPFmean_r,M.MPFlb_r,M.MPFub_r);
+plotestimates(true_R, M.MPF_apost_mean_r,M.MPF_apost_lb_r,M.MPF_apost_ub_r);
 title('Variance parameters r')
 
 subplot(2,2,3);
-plotestimates(1, M.MPFmean_phi(:,1),M.MPFlb_phi(:,1),M.MPFub_phi(:,1));
+plotestimates(1, M.MPF_apost_mean_phi(:,1),M.MPF_apost_lb_phi(:,1),M.MPF_apost_ub_phi(:,1));
 title('Forgetting factor')