%Load data and make variables
%%
load('../Data/StudData2014.txt')

hogde=StudData2014(:,1);
hogdeMor=StudData2014(:,2);
hogdeFar=StudData2014(:,3);
kjonn=StudData2014(:,4);
har=StudData2014(:,5);
auge=StudData2014(:,6);
Matlab=StudData2014(:,7);


%Make vectors with hight of men and women

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

%% Finn gjennomsnitt og varians kvinner
gjKvinner=mean(hogdeKvinne)
%varians
varKvinner=var(hogdeKvinne)
stdKvinner=sqrt(var(hogdeKvinne))


antkvinner=length(indexKvinne)


%% SME

%Data

nBort=9;
nTron=3;
p=0.8;
pvaaredata=p^nBort*(1-p)^nTron

%%

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

%%
