MATLAB kommandoer for regresjon om det skjeve t�rnet i Pisa

 

% Vi leser f�rst inn �rstallene og skjevheten i s�ylevektorer

% (Det er viktig at vi bruker s�ylevektorer her.)

 

aar=(75:87)'

skjev=[642,644,656,667,673,688,696,698,713,717,725,742,757]'

 

 

% Vi plotter skjevhet mot �r:

 

scatter(aar,skjev)

 

 

% 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(aar))".)

 

[b,bint]=regress(skjev,[ones(size(aar)),aar])

 

 

% Endelig tegner vi regresjonslinja i samme figur som observasjonene:

 

hold on

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