Analyse af solafskærmninger mht. termiske og ... - Viden om vinduer

title1 'GLM sp2_11';
*** Linear Models Analysis ***;
options pageno=1;
DATA udsyn.samlet_sort;
set udsyn.samlet;
if vejrtype=3 then delete;
if vejrtype=. then delete;
proc glm data=Udsyn.samlet_sort;
class vejrtype indstilling rum tid personnr alder;
model Sp2_11 = vejrtype indstilling(vejrtype) rum tid personnr(alder)
alder/
SS3 solution;
random personnr(alder) alder;
*Newman-Keuls test*;
means vejrtype rum tid personnr(alder) alder/SNK;
lsmeans indstilling(vejrtype)/pdiff;
* Create output data set for plots *;
output OUT=udsyn.resid_sp2_11 p=_pred r=_resid STDR=_STDR;
run; quit;
** Residual plots **;
Proc plot; plot _resid*sp2_11;
run;
quit;
proc gplot data=udsyn.resid_sp2_11;
goptions reset=symbol;
symbol1 color=black value=SQUARE height=1 cells;
label _resid = "Residualer";
label Sp2_11 = "Observeret";
plot _RESID * Sp2_11 /
vaxis=axis1 vminor=0 haxis=axis3 hminor=0 vref=0
cframe=white caxis=BLACK name="RESID"
description = "Residuals of sp2_11 by observed sp2_11";
run; quit;
filename grafout "C:\Documents and
settings\Helle\Skrivebord\Arbejdsmappe\SAS\ GLM";
goptions device=imggif gsfname=grafout gsfmode=replace ftext=SWISS
ctext=BLACK htext=1.8 cells; symbol1 c=BLACK h=1 cells;
run; quit;
PROC sort DATA=udsyn.resid_sp2_11 out=udsyn.plot_sp2_11;
by _Resid;
run;
**Kvantil plot**;
DATA udsyn.quantileplot_sp2_11;
set udsyn.plot_sp2_11;
p=(_n_ - 0.5)/822;
quantile=probit(p);
proc plot; plot quantile*_resid;
run; quit;
proc gplot data=udsyn.quantileplot_sp2_11;
goptions reset=symbol;
symbol1 color=black value=SQUARE height=1 cells;
label quantile = "Quantile";
label _resid = "Residualer";
plot quantile*_resid/
vaxis=axis1 vminor=0 haxis=axis3 hminor=0 vref=0 cframe=white
caxis=BLACK name="Quan"
description = "Quantiles som funktion af residualer sp2_11";
run; quit;
filename grafout "C:\Documents and Settings\Helle\
Skrivebord\Arbejdsmappe\SAS\GLM";
98