MATLAB kommandoer for multippel
regresjon om f�dselsvekter
% Tidligere har vi brukt
kommandoen "regress" i forbindelse med line�r regresjon.
% Vi vil n� bruke
kommandoen "regstats" som gir oss flere
muligheter enn "regress".
% Kommandoen "regstats" tar automatisk med et konstantledd i
regresjonen,
% s� vi� trenger ikke � gi dette eksplisitt
slik tilfellet er for "regress".
% Vi
leser inn dataene.
% (Kommandoen forutsetter at filen 'fvekt.txt' er
% lagret
i en mappe der MATLAB kan finne den.)
fodsel=load('fvekt.txt') ;
% Vi lager et matriseplott for � f�
oversikt over alle variablene
plotmatrix(fodsel)
% Vi gj�r s� en multippel line�r
regresjon av f�dselsvekt mot de fire andre variablene
% (mors alder, svangerskapets
varighet, barnets kj�nn og barnets paritet)
regstats(fodsel(:,5), fodsel(:,1:4))
% N�r du gir denne kommandoen,
kommer det opp et vindu det du kan angi hvilke resultater du vil ta vare p�.
% Her er det flere muligheter enn
vi vil komme inn p� i STK1110.
% Vi haker av for "Coefficients",� "Fitted
Values", "Residuals",
"R-square Statistic"
og "t Statistics" og klikker OK
% Vi f�r da beregnet og lagret
disse variablene med de navnene som er gitt i vinduet.
% Vi ser p� de estimerte regresjonskoeffisientene:
beta
% Vi f�r flere detaljer ved � se p�
"tstat".
% Denne innholder flere variable:
tstat���������������� % Gir oversikt over variablene i "tstat"
tstat.beta�������� %
Estimatene
tstat.se������������ %
Standardfeilen til estimatene
tstat.t�������������� %
t-test for nullhypotesene at regresjonskoeffisientene er lik null
tstat.pval�������� %
P-verdier for t-testene
tstat.dfe���������� %
Antall frihetsgrader (n-p)
% Vi plotter residualene
mot hver av x-variablene
scatter(fodsel(:,1), r)��� ���������� %
Plott mot mors alder
scatter(fodsel(:,2), r)��� ���������� %
Plott mot svangerskapets varighet
scatter(fodsel(:,3), r)��� ���������� %
Plott mot barnets kj�nn
scatter(fodsel(:,4), r)��� ���������� %
Plott mot barnets pariet
% Vi plotter residualene
mot det tilpassede verdiene
scatter(yhat, r)���
% Vi lager normalfordelingsplott av �residualene
probplot(r)
% Vi skriver ut R^2 (dvs "forklart" variasjon)
rsquare