% Script de Matlab/Octave -- Ricard Torres, 2015
%
% Calculos lista de ejercicios 2, MicroApl II (sin acentos)

% Matriz autentica de probabilidades
% Dimension 1: Comida (normal, organica)
% Dimension 2: Mortandad (muerte antes 60, sobrevivencia)
% Dimension 3: Factor genetico (presente, ausente)
a = zeros(2,2,2);
a(:,:,1) = [ .0008 .0092 
             .0198 .1602];
a(:,:,2) = [ .7120 .0880 
             .0099 .0001];

% Matriz agregada si no distinguimos el factor genetico
b = zeros(2,2);
b = sum(a,3) ;  % ie, b = a(:,:,1) + a(:,:,2)

% Probabilidades condicionadas de sobrevivencia para la matriz agregada,
% y para cada factor genetico
for i = 1 : 2  % 1: normal, 2: organica
  cond_agg(i) = b(i,2)/sum(b(i,:));
  cond_pres(i) = a(i,2,1)/sum(a(i,:,1));
  cond_aus(i) = a(i,2,2)/sum(a(i,:,2));
end

disp('Prob cond superv dada org si factor genetico no se toma en cuenta');
disp(cond_agg);
disp('Prob cond superv dada org si factor genetico esta presente');
disp(cond_pres);
disp('Prob cond superv dada org si factor genetico esta ausente');
disp(cond_aus);

% Para entender mejor los resultados anteriores, veamos las probabilidades 
% de comer organico cuando condicionamos en la supervivencia
disp('Prob cond comer organico dada supervivencia sin considerar factor');
disp(b(2,2)/sum(b(:,2)));
disp('Prob cond comer organico dada supervivencia con presencia factor');
disp(a(2,2,1)/sum(a(:,2,1)));
disp('Prob cond comer organico dada supervivencia con ausencia factor');
disp(a(2,2,2)/sum(a(:,2,2)));

%% Output
%Prob cond superv dada org si factor genetico no se toma en cuenta
%   0.12000   0.84368
%Prob cond superv dada org si factor genetico esta presente
%   0.92000   0.89000
%Prob cond superv dada org si factor genetico esta ausente
%   0.110000   0.010000
%Prob cond comer organico dada supervivencia sin considerar factor
% 0.62252
%Prob cond comer organico dada supervivencia con presencia factor
% 0.94569
%Prob cond comer organico dada supervivencia con ausencia factor
% 0.0011351

