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)