% Foerst importerer du doedssannsynlighetene.
% Du gjoer det ved aa (instruksjonene er testet paa pc):
%  1) Laste ned filen "doedelighet.txt" fra kurssiden
%  2) Klikke paa "Import Data" i MATLAB
%  3) Velge filen "doedelighet.txt" 
%  4) Klikke paa "Import Selection" 
%     (importer dataene som "Column vectors")
% Merk at paa filen er doedssannsynlighetene gitt i promille.

%%%% a)
% Plott doedsssannsynligheter per 1000 menn paa logaritmisk skala:
semilogy(alder, kvinner)
hold on;
semilogy(alder, menn)
hold off;

% Setter navn paa aksene:
xlabel('Alder')
ylabel('Doedsssannsynligheter per 1000') 

%%%% b)
% Bestem de ettaarige doedssannsynlighetene (qx)
q_kvinner = kvinner / 1000;
q_menn = menn / 1000;

% Bestem overlevelsessannsynlighetene 
% S(x) = P(X > x) for x = 0, 1, 2, ..., 106:
S_kvinner = cumprod(1 - q_kvinner);
S_menn = cumprod(1 - q_menn);

% Bestem den kumulative fordelingsfunksjonen F(x) = P(X <= x) = 1 - S(x)
% for x = 0, 1, 2, ..., 106:
F_kvinner = 1 - S_kvinner;
F_menn = 1 - S_menn;

% Plott av den kumulative fordelingsfunksjonen:
hold on;
stairs(alder, F_kvinner)
stairs(alder, F_menn)
hold off;

% Setter navn paa aksene 
xlabel('Alder')
ylabel('F(x)') 

%%%% c)
% P(X >= 60) = 1 - P(X <= 59)
% Siden X begynner fra 0, så finnes P(X <= 59) på plass 60 i vektorene.
P_minst60_kvinner = 1 - F_kvinner(60)
P_minst60_menn = 1 - F_menn(60)

% P(X >= 70)
P_minst70_kvinner = 1 - F_kvinner(70)
P_minst70_menn = 1 - F_menn(70)

% P(X >= 80)
P_minst80_kvinner = 1 - F_kvinner(80)
P_minst80_menn = 1 - F_menn(80)

% P(X >= 90)
P_minst90_kvinner = 1 - F_kvinner(90)
P_minst90_menn = 1 - F_menn(90)







% Bestem punktsannsynligheten p(x)=P(X=x)=F(x)-F(x-1) 
% for x=0,1,2,...,106:
px = Fx - [0; Fx(1:106)];

%Plott av punktsannsynligheten:
bar(alder, px)

% Setter navn paa aksene 
xlabel('Alder')
ylabel('p(x)') 