Auge=StudDataV2011(:,5);
Har=StudDataV2011(:,6);

subplot(2,1,1)
[nA,xoutA] = hist(Auge,1:1:4)
bar(xoutA, nA)
title('Augefarge, blå (1), brun(2), grøn(3) og grå(4)')
subplot(2,1,2)
[nH,xoutH]=hist(Har,1:1:2)
bar(xoutH,nH)
title('Hårfarge, lys blond (1) eller brun/mørk blond/svart (2)')

fblond=find(Har==1);
fmork = find(Har==2);

figure
subplot(2,1,1)
[nAblond,xoutA] = hist(Auge(fblond),1:1:4)
bar(xoutA,nAblond)
subplot(2,1,2)
[nAmork,xoutAmork] = hist(Auge(fmork),1:1:4)
bar(xoutA, nAmork)

figure()
bar3(xoutA,[nAblond;nAmork]')
colormap([0.9 0.9 0;0.2 0.1 0])
ylabel('Augefarge')
xlabel('Hårfarge')

figure()
bar3(xoutA,[nAblond;nAmork]'/length(Har))
colormap([0.9 0.9 0;0.2 0.1 0])
ylabel('Augefarge')
xlabel('Hårfarge')
title('Simultan sannsynsfordeling for hårfarge (som 5-åring) og augefarge')

sannsford=[nAblond;nAmork]'/length(Har)
%f(auge,hår)

antall=[nAblond;nAmork]'

%Finn betingafordeling, hårfarge|auge=blå
sannsford(1,:) %tar ut den rada

betford=sannsford(1,:)/sum(sannsford(1,:))


%Finn marinalfordelingar

MarginalAuge=sum(sannsford')
nA/length(Auge)
subplot(2,1,1)
bar(xoutA,MarginalAuge)
title('(Marginal)fordeling til augefarge')


MarginalHar=sum(sannsford)
nH/length(Har)
subplot(2,1,2)
bar(xoutH,MarginalHar)
title('(Marginal)fordeling til hårfarge')

%Finn betingafordeling, hårfarge|auge=blå
sannsford(1,:) %tar ut den rada

betford=sannsford(1,:)/sum(sannsford(1,:))


