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)