% hw2_11

% binomial model

clear
y_values = [20 30 40 50 60];
y_values_rep = repmat(y_values, 4,1);
f_values = [.032 .044 .027 .052]';
f_values_rep = repmat(f_values, 1,5);

    
pp = binopdf(y_values_rep,1000,f_values_rep);

p_post_1 = pp./repmat(sum(pp,1), 4,1);

disp(' ')
disp('y values:')
disp(y_values)
disp('posterior probability under uniform prior:')
disp(p_post_1)
disp(' ')


prior_values = [.3 .3 .2 .2]';
prior_values_rep = repmat(prior_values, 1,5);
pp_2 = pp.*prior_values_rep;
p_post_2 = pp_2./repmat(sum(pp_2,1), 4,1);

disp(y_values)
disp('posterior probability under non_uniform prior:')
disp(p_post_2)
disp(' ')


% Poisson model

clear
y_values = [20 30 40 50 60];
y_values_rep = repmat(y_values, 4,1);
f_values = [.032 .044 .027 .052]';
f_values_rep = repmat(f_values, 1,5);

    
pp = poisspdf(y_values_rep,1000*f_values_rep);

p_post_1 = pp./repmat(sum(pp,1), 4,1);

disp(' ')
disp('y values:')
disp(y_values)
disp('posterior probability under uniform prior:')
disp(p_post_1)
disp(' ')


prior_values = [.3 .3 .2 .2]';
prior_values_rep = repmat(prior_values, 1,5);
pp_2 = pp.*prior_values_rep;
p_post_2 = pp_2./repmat(sum(pp_2,1), 4,1);

disp(y_values)
disp('posterior probability under non_uniform prior:')
disp(p_post_2)
disp(' ')
