MATLAB kommandoer for regresjon om f?dselsvekter (II)

 

% F?rst leser vi inn dataene.

% (Kommandoen forutsetter at filen 'fvekt.txt' er

% lagret i en mappe der MATLAB kan finne den.)

 

fodsel=load('fvekt.txt') ;

 

 

% S? definerer vi en vektor for f?dselsvektene og en vektor for varigheten av svangerskapene

 

vekt=fodsel(:,5);

varighet=fodsel(:,2);

 

 

% Vi bestemmer minste kvadraters estimater, residualene, m.m.

 

[b,bint,res,resint,stats]=regress(vekt,[ones(size(varighet)),varighet]);

 

 

% Vi plotter residualene mot varighet og lager et normalfordelingplott av dem

 

scatter(varighet,res)

 

probplot(res)

 

% Variabelen 'stats' inneholder fire tall.

% Det f?rste av dem er R^2 (kvadratet av korrelasjonen)

% og det siste er S^2 (estimatet for residual variansen)

 

stats