MATLAB kommandoer for regresjon om f�dselsvekter

 

% 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 plotter f�dselsvektene mot varighet:

 

scatter(varighet,vekt)

 

 

% Vi bestemmer minste kvadraters estimater med konfidensintervall

% (Merk at vi m� si eksplisitt at vi skal ha med et konstantledd

% i regresjonen ved � skrive "ones(size(varighet))".)

 

[b,bint]=regress(vekt,[ones(size(varighet)),varighet])

 

 

% Endelig tegner vi regresjonslinja i samme figur som f�dselsvektene:

 

hold on

line(varighet,b(1)+b(2)*varighet)