
%% last inn 'spund1ed.csv' vha import data
%%
hogde=nhoyde;
kjonn = zeros(length(nhoyde),1);

for(i=1:length(nhoyde))
    if(strcmp(char(VarName2(i)),'Mann')==1)
        kjonn(i)=1;
    end
    if(strcmp(char(VarName2(i)),'Kvinne')==1)
        kjonn(i)=2;
    end
end

%% Make vectors with hight of men and women

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

%% gjennomsnitt
gjsnittKvinner=mean(hogdeKvinne)


%% varians
varKvinner=var(hogdeKvinne)


stdKvinner=sqrt(varKvinner)



%% SME

%Data

nTrondelag=1;
nUtanfor=9;
p=0.9;
pvaaredata=(1-p)^nTrondelag*p^nUtanfor

%%

%%
%Likelihoodfunksjon / Sanns for vaare data
p=0:0.01:1;
L=(1-p).^nTrondelag.*(p).^nUtanfor;
plot(p,L,'LineWidth',3);
xlabel('Andel som reiser bort (p)');
ylabel('Sannsyn for vaare data, L(p;x_1, ..., x_{n})');
t=title('Studentjul (likelihoodfunksjon)');
set(gca, 'FontSize', 20);
%%
