load('../Data/StudData2012.txt')
har=StudData2012(:,2);
auge=StudData2012(:,3);


subplot(2,1,1)
[nA,xoutA] = hist(auge,1:1:3)
bar(xoutA, nA)
title('Augefarge, blå (1), brun(2) og grøn(3)')
subplot(2,1,2)
[nH,xoutH]=hist(har,1:1:4)
bar(xoutH,nH)
title('Hårfarge, lys blond (1), brun (2), svart (3) og raud(4)')

flys=find(har==1);
fbrun = find(har==2);
fsvart=find(har==3);
fraud = find(har==4);

figure
subplot(2,2,1)
[nAlys,xoutA] = hist(auge(flys),1:1:3);
bar(xoutA,nAlys)
xlabel('Augefarge')
title('Hårfarge lys')
subplot(2,2,2)
[nAbrun,xoutAbrun] = hist(auge(fbrun),1:1:3)
bar(xoutA, nAbrun)
xlabel('Augefarge')
title('Hårfarge brun')
subplot(2,2,3)
[nAsvart,xoutA] = hist(auge(fsvart),1:1:3);
bar(xoutA,nAsvart)
xlabel('Augefarge')
title('Hårfarge svart')
subplot(2,2,4)
[nAraud,xoutA] = hist(auge(fraud),1:1:3);
bar(xoutA, nAraud)
xlabel('Augefarge')
title('Hårfarge raud')


figure()
bar3(xoutA,[nAlys;nAbrun;nAsvart;nAraud]')
colormap([0.9 0.9 0;0.2 0.1 0;0,0,0;1,0,0])
ylabel('Augefarge')
xlabel('Hårfarge')

figure()
bar3(xoutA,[nAlys;nAbrun;nAsvart;nAraud]'/length(har))
colormap([0.9 0.9 0;0.2 0.1 0;0,0,0;1,0,0])
ylabel('Augefarge')
xlabel('Hårfarge')
title('Simultan sannsynsfordeling for hårfarge (som 5-åring) og augefarge')

sannsford=[nAlys;nAbrun;nAsvart;nAraud]'/length(har)
%f(auge,hår)

antall=[nAlys;nAbrun;nAsvart;nAraud]'

%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')



