*** SAS file generated by sastangle ***; options nocenter nodate nonumber ls=80 formdlim=''; *** Code chunk number 1 ***; proc import datafile='c:/book/help.dta' out=help_a dbms=dta; run; data help; set help_a; if female; run; *** Code chunk number 2 ***; axis1 order = (0 to 40 by 10) minor=none; axis2 minor=none; legend1 label=none value=(h=1.5) shape=symbol(10,1.2) down=3 position=(top right inside) frame mode=protect; symbol1 v=circle i=sm70s c=black l=1 h=1.1 w=5; symbol2 v=diamond i=sm70s c=black l=33 h=1.1 w=5; symbol3 v=square i=sm70s c=black l=8 h=1.1 w=5; proc gplot data=help; plot i1*age = substance / vaxis=axis1 haxis=axis2 legend=legend1; run; quit; *** Code chunk number 4 ***; options ls=74; /* reduces width of output to make it fit in gray area */ ods select overallanova; proc glm data=help; class substance; model i1 = age substance age * substance / solution; output out=helpout cookd=cookd_ch4 dffits=dffits_ch4 student=sresids_ch4 residual=resid_ch4 predicted=pred_ch4 h=lev_ch4; run; quit; ods select all; *** Code chunk number 5 ***; options ls=74; /* reduces width of output to make it fit in gray area */ ods select modelanova; proc glm data=help; class substance; model i1 = age substance age * substance / solution; output out=helpout cookd=cookd_ch4 dffits=dffits_ch4 student=sresids_ch4 residual=resid_ch4 predicted=pred_ch4 h=lev_ch4; run; quit; ods select all; *** Code chunk number 6 ***; options ls=74; /* reduces width of output to make it fit in gray area */ ods select parameterestimates; proc glm data=help; class substance; model i1 = age substance age * substance / solution; output out=helpout cookd=cookd_ch4 dffits=dffits_ch4 student=sresids_ch4 residual=resid_ch4 predicted=pred_ch4 h=lev_ch4; run; quit; ods select all; *** Code chunk number 7 ***; ods select none; ods output parameterestimates=helpmodelanova; proc glm data=help; class substance; model i1 = age|substance / solution; output out=helpout cookd=cookd_ch4 dffits=dffits_ch4 student=sresids_ch4 residual=resid_ch4 predicted=pred_ch4 h=lev_ch4; run; quit; ods select all; *** Code chunk number 8 ***; proc print data=helpmodelanova; var parameter estimate stderr tvalue probt; format _numeric_ 6.3; run; *** Code chunk number 9 ***; options ls=74; proc means data=helpout min q1 median q3 max maxdec=2; var resid_ch4; run; *** Code chunk number 10 ***; proc sort data=helpout; by resid_ch4; run; proc print data=helpout (obs=1); var id age i1 substance pred_ch4 resid_ch4; run; *** Code chunk number 11 ***; proc sort data=helpout; by descending resid_ch4; run; proc print data=helpout (obs=1); var id age i1 substance pred_ch4 resid_ch4; run; *** Code chunk number 12 ***; goptions reset=all; *** Code chunk number 13 ***; ods select univar; proc univariate data=helpout; qqplot resid_ch4 / normal(mu=est sigma=est color=black); run; ods select all; *** Code chunk number 14 ***; axis1 label=("Standardized residuals"); ods select "Histogram 1"; proc univariate data=helpout; var sresids_ch4; histogram sresids_ch4 / normal(mu=est sigma=est color=black) kernel(color=black) haxis=axis1; run; ods select all; *** Code chunk number 15 ***; ods select none; proc sort data=help; by substance; run; ods output parameterestimates=helpsubstparams; proc glm data=help; by substance; model i1 = age / solution; run; ods select all; *** Code chunk number 16 ***; options ls=74; proc print data=helpsubstparams; run; *** Code chunk number 17 ***; libname k 'c:/book'; proc sort data=k.help; by substance female; run; ods select none; proc means data=k.help; by substance female; var cesd; output out=helpmean mean=; run; ods select all; *** Code chunk number 18 ***; axis1 minor=none; symbol1 i=j v=none l=1 c=black w=5; symbol2 i=j v=none l=2 c=black w=5; proc gplot data=helpmean; plot cesd*substance = female / haxis=axis1 vaxis=axis1; run; quit; *** Code chunk number 19 ***; data h2; set k.help; if female eq 1 then sex='F'; else sex='M'; run; proc sort data=h2; by sex; run; symbol1 v='x' c=black; proc boxplot data=h2; plot cesd * substance(sex) / notches boxwidthscale=1; run; *** Code chunk number 20 ***; options ls=74; ods select modelanova; proc glm data=k.help; class female substance; model cesd = female substance female*substance / ss3; run; *** Code chunk number 22 ***; options ls=74; /* stay in gray box */ ods select overallanova; proc glm data=k.help; class female substance; model cesd = female substance / ss3 solution; run; *** Code chunk number 23 ***; options ls=74; /* stay in gray box */ ods select parameterestimates; proc glm data=k.help; class female substance; model cesd = female substance / ss3 solution; run; *** Code chunk number 24 ***; ods select fitstatistics; proc mixed data=k.help method=ml; class female substance; model cesd = female|substance; run; quit; *** Code chunk number 25 ***; ods select fitstatistics; proc mixed data=k.help method=ml; class female substance; model cesd = female substance; run; quit; ods select all; *** Code chunk number 27 ***; ods select diff; proc glm data=k.help; class substance; model cesd = substance; lsmeans substance / pdiff adjust=tukey cl lines; run; quit; ods select all; *** Code chunk number 28 ***; ods select lsmeandiffcl; proc glm data=k.help; class substance; model cesd = substance; lsmeans substance / pdiff adjust=tukey cl lines; run; quit; ods select all; *** Code chunk number 29 ***; ods select lsmlines; proc glm data=k.help; class substance; model cesd = substance; lsmeans substance / pdiff adjust=tukey cl lines; run; quit; ods select all; *** Code chunk number 30 ***; ods select diffplot; proc glm data=k.help; class substance; model cesd = substance; lsmeans substance / pdiff adjust=tukey cl lines; run; quit; ods pdf close; ods graphics off; *** Code chunk number 32 ***; ods select contrastcoef; proc glm data=k.help; class female substance; model cesd = female substance; output out=outanova residual=resid_ch4anova; estimate 'A+H = C?' substance 1 -2 1 / e; run; quit; ods select all; *** Code chunk number 33 ***; options ls=74; ods select estimates; proc glm data=k.help; class female substance; model cesd = female substance; output out=outanova residual=resid_ch4anova; estimate 'A+H = C?' substance 1 -2 1 / e; run; quit; ods select all;