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)