%Load data and make variables

load('../Data/KjonnOgHoyde.mat')

kjonn=KjonnOgHoyde(:,1);
hogde=KjonnOgHoyde(:,2);


%% Make vectors with hight of men and women

indexMann=find(kjonn==0);
indexKvinne=find(kjonn==1);
hogdeMann=hogde(indexMann);
hogdeKvinne=hogde(indexKvinne);

antkvinner=length(indexKvinne)


%% gjennomsnitt
gjsnittKvinner=mean(hogdeKvinne)


%% varians
varKvinner=var(hogdeKvinne)


stdKvinner=sqrt(varKvinner)



%% SME

%Data

nTrondelag=2;
nUte=8;
p=0.9;
pvaaredata=p^nUte*(1-p)^nTrondelag

%%

%%
%Likelihoodfunksjon
p=0:0.01:1;
L=p.^nUte.*(1-p).^nTrondelag;
plot(p,L,'LineWidth',2)
xlabel('p')
ylabel('L(p;x_1, ..., x_{10})')
title('Likelihoodfunksjon, studentjul')

%%
