dm 'log;clear;output;clear'; 
options nodate nocenter nonumber ls=90 ps=512;
OPTIONS FORMCHAR="|----|+|---+=|-/\<>*";
ods html style=minimal body='Assignment 10 Example.html'; 
ods graphics on; 
ODS listing;   

************************************************;
*** James P Geaghan                          ***;
*** EXST 700X Analysis of Variance Examples  ***;
*** Assignment 10                            ***;
**********************************************************************;
*** Examples are from your text book                               ***;
*** Freund, Rudolph J., William J. Wilson and Donna L. Mohr. 2010. ***;
*** Statistical Methods, Academic Press (ELSIVIER), N.Y            ***;
**********************************************************************;
* Serious environmental problems arise from absorption into soil of  *;
* metals that escape into the air from different industrial          *;
* operations. To ascertain if absorption rates differ among soil     *;
* types, six soil samples were randomly selected from fields having  *;
* five different soil types (A, B, C, D, and E) in an area known to  *;
* have  relatively uniform exposure to the metals studied. The 30    *;
* soil samples were analyzed for cadmium (Cd) and lead (Pb) content. *;
* Data for the two metals were to be analyzed separately.            *;
**********************************************************************;

title1 "Assignment 10 - Cadmium content in soil example"; 

data MetalContent; 
   INFILE 'datatab_6_37.csv' dlm=',' dsd missover firstobs=2;
	  input soil $ Cadmium Lead;
datalines; run;
;
run;

*proc print data=MetalContent; run;

PROC BOXPLOT data=MetalContent; 
   plot Cadmium*soil; 
run;

PROC MIXED data=MetalContent cl covtest; 
   TITLE2 'Analysis of Variance using PROC MIXED'; 
   class soil;
   model Cadmium = soil / outp=outputstuff; 
   lsmeans soil / adjust=tukey cl;
 ods output diffs=ppp lsmeans=mmm;
 *ods exclude diffs lsmeans;
run;
%include 'pdmix800.sas';
%pdmix800(ppp,mmm,alpha=0.05,sort=yes);
RUN;

PROC UNIVARIATE data=outputstuff normal plots; 
   TITLE3 'Test of residuals from ANOVA'; 
   VAR resid;
   ods exclude BasicMeasures ExtremeObs ExtremeValues Modes
               MissingValues Quantiles TestsForLocation; 
RUN;

PROC MIXED data=MetalContent cl covtest; 
    TITLE2 'PROC MIXED with GROUP statement'; 
    class soil;
    model Cadmium = soil / outp=outputstuff ddfm=kr; 
    repeated / group=soil;
    lsmeans soil / adjust=tukey cl;
ods output diffs=ppp lsmeans=mmm;
 *ods exclude diffs lsmeans;
run;
%include 'pdmix800.sas';
%pdmix800(ppp,mmm,alpha=0.05,sort=yes);
RUN;




/*

*** the same analysis is run on the second variable below ***;
title1 "Assignment 10 - Lead content in soil example"; 

PROC BOXPLOT data=MetalContent; plot Lead*soil; run;

PROC MIXED data=MetalContent; 
   TITLE2 'Analysis of Variance using PROC MIXED'; 
   class soil;
   model Lead = soil / outp=outputstuff; 
   lsmeans soil / adjust=tukey cl;
 ods output diffs=ppp lsmeans=mmm;
 *ods exclude diffs lsmeans;
run;
%include 'pdmix800.sas';
%pdmix800(ppp,mmm,alpha=0.05,sort=yes);
RUN;

PROC UNIVARIATE data=outputstuff normal plots; 
   TITLE2 'Test of residuals from ANOVA'; 
   VAR resid;
   ods exclude BasicMeasures ExtremeObs ExtremeValues Modes
               MissingValues Quantiles TestsForLocation; 
RUN;

PROC MIXED data=MetalContent; 
   TITLE2 'Analysis of Variance using PROC MIXED'; 
   class soil;
   model Lead = soil / outp=outputstuff; 
   repeated / group=soil;
RUN;


*/

ods html close;