%%
load('../Data/StudDataV16.txt')
%%
hogde=StudDataV16(:,1);
hogdeFar=StudDataV16(:,2);
hogdeMor=StudDataV16(:,3);
kjonn=StudDataV16(:,4);
har=StudDataV16(:,5);
auge=StudDataV16(:,6);

%% Make vectors with hight of men and women

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

antmenn=length(indexMann)
antkvinner=length(indexKvinne)



%% gjennomsnitt
gjsnittKvinner=mean(hogdeKvinne)


%% varians
varKvinner=var(hogdeKvinne)


stdKvinner=sqrt(varKvinner)



%% SME

%Data

nUtland=2;
nNorge=19;
p=0.9;
pvaaredata=p^nNorge*(1-p)^nUtland

%%

%%
%Likelihoodfunksjon
p=0:0.01:1;
L=p.^nNorge.*(1-p).^nUtland;
plot(p,L,'LineWidth',3);
set(gca,'FontSize',16);
t=xlabel('p');
set(t, 'FontSize', 20);
t=ylabel('L(p;x_1, ..., x_{n})');
set(t, 'FontSize', 20);
t=title('Likelihoodfunksjon, student p?ske');
set(t, 'FontSize', 20);
%%
