MATLAB kommandoer ved regresjon om det skjeve t?rnet i Pisa (II)

 

% Vi leser inn ?rstallene og skjevheten i s?ylevektorer

 

aar=(75:87)'

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

 

% Vi bestemmer minste kvadraters estimater, residualene, m.m.

 

[b,bint,res,resint,stats]=regress(skjev,[ones(size(aar)),aar]);

 

% Vi plotter residualene mot ?r og lager et normalfordelingplott av dem

 

scatter(aar,res)

 

probplot(res)

 

% Variabelen 'stats' inneholder fire tall.

% Det f?rste av dem er R^2 (kvadratet av korrelasjonen)

% og det siste er S^2 (estimatet for residual variansen)

 

stats