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