Changeset 924 for applications
- Timestamp:
- 05/05/10 20:51:26 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
applications/bdmtoolbox/tutorial/userguide/frg_estim.m
r904 r924 11 11 A1.rv = y; 12 12 A1.rgr = RVtimes([y,u],[-3,-1]) ; 13 A1.log_level = 'logbounds ,logevidence';13 A1.log_level = 'logbounds'; 14 14 A1.frg = 0.98; 15 15 A1.name = 'A1'; … … 35 35 E.prior.class = 'eDirich'; % prior on non-linear part 36 36 E.prior.beta = [2 1]; % 37 E.log_level = 'logbounds ,logevidence';37 E.log_level = 'logbounds'; 38 38 E.name = 'MPF'; 39 39 … … 42 42 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 43 43 % plot results 44 ndat = size(M.DS_ u,1);44 ndat = size(M.DS_dt_u,1); 45 45 46 46 figure(1); 47 47 subplot(2,2,1); 48 plotestimates(true_theta, M.MPF_apost_mean_theta, M.MPF_apost_lb _theta, M.MPF_apost_ub_theta);48 plotestimates(true_theta, M.MPF_apost_mean_theta, M.MPF_apost_lbound_theta, M.MPF_apost_ubound_theta); 49 49 title(' Regression parameters \theta') 50 50 set(gca,'YLim',[-1.5,1]); 51 51 52 52 subplot(2,2,2); 53 plotestimates(true_R, M.MPF_apost_mean_r,M.MPF_apost_lb _r,M.MPF_apost_ub_r);53 plotestimates(true_R, M.MPF_apost_mean_r,M.MPF_apost_lbound_r,M.MPF_apost_ubound_r); 54 54 title('Variance parameters r') 55 55 56 56 subplot(2,2,3); 57 plotestimates(1, M.MPF_apost_mean_phi(:,1),M.MPF_apost_lb _phi(:,1),M.MPF_apost_ub_phi(:,1));57 plotestimates(1, M.MPF_apost_mean_phi(:,1),M.MPF_apost_lbound_phi(:,1),M.MPF_apost_ubound_phi(:,1)); 58 58 title('Forgetting factor') 59 59