%
% Written by:
% -- 
% John L. Weatherwax                2008-02-20
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

close all; 

tmp = []; 
for xi = 10:1000,
  tmp=[tmp;cmpt_ml_var_bias(xi)];
end
figure; hist(tmp); 
xlabel('bias=E[T]-1'); ylabel('number of samples'); 

%title( 'unbiased estimate of variance' );
%saveas( gcf, 'unbiased_est_var', 'epsc' ); 

title( 'biased estimate of variance' );
saveas( gcf, 'biased_est_var', 'epsc' );