%Load data and make variables
%%
load('../Data/StudDataV16.txt')

hogde=StudDataV16(:,1);
hogdeFar=StudDataV16(:,2);
hogdeMor=StudDataV16(:,3);
kjonn=StudDataV16(:,4)-1;


%% Finn dei som har begge foreldre og seg sjolv, og justerer kvinnehogde til mannshogde

f=find(1-isnan(hogdeMor) & (1-isnan(hogdeFar)) & (1-isnan(hogde)));
% Hogde barn, jenter justert til gutar
hogdeBarn=hogde(f).*(1+0.08*(kjonn(f)));

%Gj.snitt hogda til foreldra. Mor*1.08
hogdeForeldre=0.5*hogdeFar(f)+0.5*1.08*hogdeMor(f);

h=plot(hogdeForeldre, hogdeBarn,'.','MarkerSize',12);
set(gca,'FontSize',18);
xlabel('Hogde foreldre','FontSize',20)
ylabel('Hogde barn','FontSize',20)

%%

%% estimerer regresjonsparametra, og plotter
[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)

a=b(1)
bb=b(2)
sigmae=stats(4)
R2=stats(1)


%% Finn PI for mangen x0 og plottar
x0=160:1:200;
ypred=a+bb*x0;
predvar=sigmae*(1+1/length(hogdeForeldre)+(x0-mean(hogdeBarn)).^2/Sxx);
PInedre=ypred-1.96*sqrt(predvar);
PIovre=ypred+1.96*sqrt(predvar);
plot(x0,PIovre,'g')
plot(x0,PInedre,'g')

%% Finn KI for mangen x0 og plottar
x0=160:1:200;
mu=a+bb*x0;
muvar=sigmae*(1/length(hogdeForeldre)+(x0-mean(hogdeBarn)).^2/Sxx);
KInedre=mu-1.96*sqrt(muvar);
KIovre=mu+1.96*sqrt(muvar);
plot(x0,KIovre,'r--')
plot(x0,KInedre,'r--')

%%

figure

[b,bint,r,rint,stats]=regress(hogdeBarn,[ones(length(hogdeForeldre),1),hogdeForeldre]);

h=plot(r,'x','MarkerSize',16);
set(gca,'FontSize',18);
xlabel('Residual nr.(i)','FontSize',20)
ylabel('Residual = y_i-(a+bx_i)','FontSize',20)


%%
figure

hist(r);
set(gca,'FontSize',18);

%%


figure

normplot(r)

set(gca,'FontSize',18);
%% Hogde mor og hogde far

mor=hogdeMor(f);
far=hogdeFar(f);

gjmor=mean(mor)
gjfar=mean(far)
stdmor=std(mor)
stdfar=std(far)
corr(mor,far)

figure
plot(far,mor,'.','MarkerSize',16);
set(gca,'FontSize',18);

xlabel('Hogde far','FontSize',20)
ylabel('Hogde mor','FontSize',20)

%% Model for mors hogde, bruker fars hogde som forklaringsvariabel

[b,bint,r,rint,stats]=regress(mor,[ones(length(far),1),far]);
hold
plot([165,205],[b(1)+b(2)*165, b(1)+b(2)*205],'r','LineWidth',2)
a=b(1)
bb=b(2)
bint

%%
figure

plot(r,'.')
%% Model for mors hogde, bruker fars hogde som forklaringsvariabel

[b,bint,r,rint,stats]=regress(far,[ones(length(far),1),mor]);


normplot(r)


