load('/home/ingelins/Documents/Fag/Statistikk/TMA4245V11/Data/StudDataV2011.txt')
%Finn indeks til med oppgitt høgde på både mor og far
f=find(1-isnan(StudDataV2011(:,3)) & (1-isnan(StudDataV2011(:,4))));
kjonn=StudDataV2011(f,1)-1;
%Hogde barn, jenter justert til gutar
hogdeBarn=StudDataV2011(f,2).*(1+0.08*kjonn);
%Gj.snitt hogda til foreldra. Mor*1.08
hogdeForeldre=0.5*StudDataV2011(f,3)+0.5*1.08*StudDataV2011(f,4);

plot(hogdeForeldre, hogdeBarn,'.')
xlabel('Høgde foreldre','FontSize',16)
ylabel('Høgde barn','FontSize',16)
%estimerer regresjonsparametra
[b,bint,r,rint,stats]=regress(hogdeBarn,[ones(length(f),1),hogdeForeldre]);
hold
plot([160,200],[b(1)+b(2)*160, b(1)+b(2)*200],'r','LineWidth',2)

gtext('y=a+bx','Color','r','FontSize',12)