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)