%% 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
%% Plot hight
hist(hogde)
set(gca,'FontSize',20)

%%
boxplot(hogde)
set(gca,'FontSize',20)
%%
boxplot(hogde,kjonn)
set(gca,'FontSize',20)
%% Make vectors with hight of men and women

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

%% histograms
subplot(1,2,1)
hist(hogdeMann)
title('Hogde menn')
set(gca,'FontSize',20)
subplot(1,2,2)
hist(hogdeKvinne)
title('Hogde kvinner')
set(gca,'FontSize',20)
%% mean
mean(hogdeMann)
mean(hogdeKvinne)

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


%%

% Antar at hogde er normalfordelt med
mu=179.8
sigma=6.5;

% Lager vektor vi skal plotte for
del=8*sigma/1000;
xvec=[mu-3*sigma:del:mu+3*sigma];

pdfpopulasjon=normpdf(xvec,mu,sigma);

figure
plot(xvec, pdfpopulasjon, 'k','LineWidth',3)
set(gca,'FontSize',20)
hold
%% Utvalgsfordeling gj.snitt av 5
pdfgj5=normpdf(xvec,mu,sigma/sqrt(5));
plot(xvec, pdfgj5, 'b','LineWidth',3 )
set(gca,'FontSize',20)
%% Utvalgsfordeling gj.snitt av 20
pdfgj20=normpdf(xvec,mu,sigma/sqrt(20));
plot(xvec, pdfgj20, 'g','LineWidth',3 )
set(gca,'FontSize',20)
%% Utvalgsfordeling gj.snitt av antmenn
pdfgjantmenn=normpdf(xvec,mu,sigma/sqrt(antmenn));
plot(xvec, pdfgjantmenn, 'r','LineWidth',3 )
set(gca,'FontSize',20)

h_legend=legend('Populasjon','n=5','n=20','n=antmenn (216)');
set(h_legend,'FontSize',18);

t=title('Fordelinga til populasjon og gjennomsnittet');
set(t,'FontSize',18);
set(gca,'FontSize',20)
%% Sample datasets with antmenn persons from distribution of hights for young Norwegian men, and find
%means

figure

nsample=200; %Number of means
samplemeans=zeros(nsample,1);
for(i=1:nsample)
    ASample=random('Normal',mu,sigma,1,antmenn);
    subplot(1,2,1)
    hist(ASample);
    t=title('Eit utvalg');
    set(t, 'FontSize', 18);
    set(gca,'FontSize',20)
    samplemeans(i)=mean(ASample);
    subplot(1,2,2) 
    hist(samplemeans(1:i));
    t=title('Utvalsfordeling for gjennomsnitt');
    set(get(gca,'child'),'FaceColor','r','EdgeColor','k');
    set(t, 'FontSize', 18);
    set(gca,'FontSize',20)
    pause(0.2)
end

%% Rask
nsample=10000;
samplemeans=zeros(nsample,1);
Samples=random('Normal',mu,sigma,nsample,antmenn);
samplemeans=mean(Samples');

figure
hist(samplemeans);
 title('Utvalsfordeling for gjennomsnitt');
set(gca,'FontSize',20)

f=find(samplemeans>mean(hogdeMann));
length(f)/nsample





%% Probability hypothesis test
1-normcdf(mean(hogdeMann),mu,sqrt(sigma^2/antmenn))
1-normcdf(mean(hogdeKvinne),mu,sqrt(sigma^2/antkvinner))



%% Demonstrasjon Sentralgrenseteoremet Simulerer med Poisson

nsample=10000; %Number of datasets
ineachsample=20;
samplemeans=zeros(nsample,1);
Samples=random('Poisson',1,nsample,ineachsample);
samplemeans=mean(Samples');

figure
subplot(1,2,1)
hist(Samples(1,:))
t=title('Eit utvalg, Poisson')
set(t, 'FontSize', 18);
subplot(1,2,2)
hist(samplemeans)
t=title('Utvalsfordeling for gjennomsnitt, Poisson')
set(t, 'FontSize', 18);
set(get(gca,'child'),'FaceColor','r','EdgeColor','k');

%% Simulerer med eksponensial
nsample=10000; %Number of datasets
ineachsample=100;
samplemeans=zeros(nsample,1);
Samples=random('exp',1,nsample,ineachsample);
samplemeans=mean(Samples');
figure
subplot(1,2,1)
hist(Samples(1,:))
t=title('Eit utvalg, eksponensial')
set(t, 'FontSize', 18);
subplot(1,2,2)
hist(samplemeans)
t=title('Utvalsfordeling for gjennomsnitt, eksponensial')
set(get(gca,'child'),'FaceColor','r','EdgeColor','k');
    set(t, 'FontSize', 18);
