- Timestamp:
- 05/11/10 16:25:24 (14 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
applications/bdmtoolbox/tutorial/userguide/dist_estim_example.m
r871 r933 27 27 A1.rv = y; 28 28 A1.rgr = RVtimes([y,u1],[-3,0]) ; % correct structure is {y,y} 29 A1. options='logbounds,logevidence';29 A1.log_level ='logbounds,logevidence'; 30 30 A1.frg = 0.95; 31 31 A1.constant = 0; … … 38 38 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 39 39 % plot results 40 ndat = size(M.DS_ u1,1);40 ndat = size(M.DS_dt_u1,1); 41 41 42 42 true_theta1 = fy.A([1,2]); … … 46 46 figure(1); 47 47 subplot(2,1,1); 48 plot(M.DS_ y);48 plot(M.DS_dt_y); 49 49 title('Output'); 50 50 51 51 subplot(2,1,2); 52 plot(M.DS_ u1);52 plot(M.DS_dt_u1); 53 53 hold on 54 plot(M.DS_ u2);54 plot(M.DS_dt_u2); 55 55 title('Input'); 56 56 … … 62 62 plotestimates(true_theta1, ... 63 63 M.Est0_apost_mean_theta, ... 64 M.Est0_apost_lb _theta, ...65 M.Est0_apost_ub _theta);64 M.Est0_apost_lbound_theta, ... 65 M.Est0_apost_ubound_theta); 66 66 set(gca,'YLim',[-1.5,1]); 67 67 … … 70 70 plotestimates(true_R, ... 71 71 M.Est0_apost_mean_r, ... 72 M.Est0_apost_lb _r, ...73 M.Est0_apost_ub _r);72 M.Est0_apost_lbound_r, ... 73 M.Est0_apost_ubound_r); 74 74 75 75 title('Variance parameters r') … … 79 79 plotestimates(true_theta2, ... 80 80 M.Est1_apost_mean_theta, ... 81 M.Est1_apost_lb _theta, ...82 M.Est1_apost_ub _theta);81 M.Est1_apost_lbound_theta, ... 82 M.Est1_apost_ubound_theta); 83 83 set(gca,'YLim',[-1.5,1]); 84 84 … … 87 87 plotestimates(true_R, ... 88 88 M.Est1_apost_mean_r, ... 89 M.Est1_apost_lb _r, ...90 M.Est1_apost_ub _r);89 M.Est1_apost_lbound_r, ... 90 M.Est1_apost_ubound_r); 91 91 92 92 title('Variance parameters r')