ODS RTF FILE = 'C:\Documents and Settings\richc\My Documents\CPH931F08\WA2931F08SAS.rtf';

PROC IMPORT DATAFILE = 'C:\Documents and Settings\richc\My Documents\CPH931F08\DSST.xls'
     OUT = DSST  DBMS = EXCEL REPLACE;
     SHEET = DSST;
     GETNAMES = YES;
RUN;

PROC MIXED DATA = DSST method = ML covtest;
CLASS SUBJECT DRUG ACTIVITY SSSTAT TIME;
MODEL CORRECTTRIALS = DRUG ACTIVITY TIME SSSTAT
		      SSSTAT*DRUG SSSTAT*ACTIVITY SSSTAT*TIME CORRECTTRIALS_B
			 DRUG*ACTIVITY SSSTAT*DRUG*ACTIVITY  / SOLUTION;
RANDOM SUBJECT; 
ESTIMATE 'a: Low SS Placebo Low Activity' 	Intercept 1  TIME 1 0 0   CORRECTTRIALS_B 60
					DRUG 1 0 ACTIVITY 1 0  SSSTAT 1 0
					DRUG*ACTIVITY 1 0 0 0 
					DRUG*ACTIVITY*SSSTAT  1 0 0 0 0 0 0 0
					DRUG*SSSTAT  1 0 0 0	
					ACTIVITY*SSSTAT  1 0 0 0 
					SSSTAT*TIME   1 0 0 0 0 0 ;
ESTIMATE 'b: Low SS Placebo High Activity' 	Intercept 1  TIME 1 0 0   CORRECTTRIALS_B 60
					DRUG 1 0 ACTIVITY 0 1 SSSTAT 1 0
					DRUG*ACTIVITY 0 1 0 0 
					DRUG*ACTIVITY*SSSTAT  0 0 1 0 0 0 0 0
					DRUG*SSSTAT  1 0 0 0	
					ACTIVITY*SSSTAT  0 0 1 0 
					SSSTAT*TIME   1 0 0 0 0 0 ;
ESTIMATE 'c: Low SS Drug Low Activity' 	Intercept 1  TIME 1 0 0   CORRECTTRIALS_B 60
					DRUG 0 1 ACTIVITY 1 0  SSSTAT 1 0
					DRUG*ACTIVITY 0 0 1 0 
					DRUG*ACTIVITY*SSSTAT   0 0 0 0 1 0 0 0
					DRUG*SSSTAT  0 0 1 0
					ACTIVITY*SSSTAT  1 0 0 0 
					SSSTAT*TIME   1 0 0 0 0 0 ;	
ESTIMATE 'd: Low SS Drug High Activity'		Intercept 1  TIME 1 0 0   CORRECTTRIALS_B 60
					DRUG 0 1 ACTIVITY 0 1  SSSTAT 1 0
					DRUG*ACTIVITY 0 0 0 1 
					DRUG*ACTIVITY*SSSTAT   0 0 0 0 0 0 1 0
					DRUG*SSSTAT  0 0 1 0  
					ACTIVITY*SSSTAT  0 0 1 0 
					SSSTAT*TIME   1 0 0 0 0 0 ;		
ESTIMATE 'e: High SS Placebo Low Activity' 	Intercept 1  TIME 1 0 0   CORRECTTRIALS_B 60
					DRUG 1 0 ACTIVITY 1 0  SSSTAT 0 1
					DRUG*ACTIVITY 1 0 0 0 
					DRUG*ACTIVITY*SSSTAT  0 1 0 0 0 0 0 0
					DRUG*SSSTAT  0 1 0 0	
					ACTIVITY*SSSTAT  0 1 0 0 
					SSSTAT*TIME   0 0 0 1 0 0 ;
ESTIMATE 'f: High SS Placebo High Activity' 	Intercept 1  TIME 1 0 0   CORRECTTRIALS_B 60
					DRUG 1 0 ACTIVITY 0 1 SSSTAT 0 1
					DRUG*ACTIVITY 0 1 0 0 
					DRUG*ACTIVITY*SSSTAT  0 0 0 1 0 0 0 0
					DRUG*SSSTAT  0 1 0 0	
					ACTIVITY*SSSTAT  0 0 0 1 
					SSSTAT*TIME   0 0 0 1 0 0 ;
ESTIMATE 'g: High SS Drug Low Activity' 	Intercept 1  TIME 1 0 0   CORRECTTRIALS_B 60
					DRUG 0 1 ACTIVITY 1 0  SSSTAT 0 1
					DRUG*ACTIVITY 0 0 1 0 
					DRUG*ACTIVITY*SSSTAT   0 0 0 0 0 1 0 0
					DRUG*SSSTAT  0 0 0 1
					ACTIVITY*SSSTAT  0 1 0 0  
					SSSTAT*TIME   0 0 0 1 0 0 ;	
ESTIMATE 'h: High SS Drug High Activity'	Intercept 1  TIME 1 0 0   CORRECTTRIALS_B 60
					DRUG 0 1 ACTIVITY 0 1  SSSTAT 0 1
					DRUG*ACTIVITY 0 0 0 1 
					DRUG*ACTIVITY*SSSTAT   0 0 0 0 0 0 0 1
					DRUG*SSSTAT  0 0 0 1 
					ACTIVITY*SSSTAT  0 0 0 1  
					SSSTAT*TIME   0 0 0 1 0 0 ;	
ESTIMATE 'DRUG/ACTIVITY interaction for Low SS: a+d-b-c'
					DRUG*ACTIVITY 1 -1 -1 1 
					DRUG*ACTIVITY*SSSTAT  1 0 -1 0 -1 0 1 0;
ESTIMATE 'DRUG/ACTIVITY interaction for High SS: e+h-f-g'
					DRUG*ACTIVITY 1 -1 -1 1 
					DRUG*ACTIVITY*SSSTAT  0 1 0 -1 0 -1 0 1;
CONTRAST 'Test No DRUG/ACTIVITY interaction in either SS stratum: a+d-b-c=0, e+h-f-g=0'
	DRUG*ACTIVITY 1 -1 -1 1 	DRUG*ACTIVITY*SSSTAT  1 0 -1 0 -1 0 1 0,
	DRUG*ACTIVITY 1 -1 -1 1 	DRUG*ACTIVITY*SSSTAT  0 1 0 -1 0 -1 0 1;
CONTRAST 'Test for relevance of ACTIVITY when DRUG = 0 and SSSTAT = 0: a-b=0'
	ACTIVITY  1 -1  DRUG*ACTIVITY  1 -1 0 0   
	ACTIVITY*SSSTAT  1 0 -1 0  DRUG*ACTIVITY*SSSTAT  1 0 -1 0 0 0 0 0 ;
CONTRAST 'Test for relevance of ACTIVITY when DRUG = 1 and SSSTAT = 0: c-d=0'
	ACTIVITY  1 -1  DRUG*ACTIVITY  0 0 1 -1   
	ACTIVITY*SSSTAT  1 0 -1 0  DRUG*ACTIVITY*SSSTAT  0 0 0 0 1 0 -1 0 ;
CONTRAST 'Test for relevance of ACTIVITY when DRUG = 0 and SSSTAT = 1: e-f=0'
	ACTIVITY  1 -1  DRUG*ACTIVITY  1 -1 0 0   
	ACTIVITY*SSSTAT  0 1 0 -1  DRUG*ACTIVITY*SSSTAT  0 1 0 -1 0 0 0 0 ;
CONTRAST 'Test for relevance of ACTIVITY when DRUG = 1 and SSSTAT = 1: g-h=0'
	ACTIVITY  1 -1  DRUG*ACTIVITY  0 0 1 -1  
	ACTIVITY*SSSTAT  0 1 0 -1  DRUG*ACTIVITY*SSSTAT  0 0 0 0 0 1 0 -1  ;
CONTRAST 'Test for overall relevance of ACTIVITY: a-b=0, c-d=0, e-f=0, g-h=0'
	ACTIVITY  1 -1  DRUG*ACTIVITY  1 -1 0 0   
	ACTIVITY*SSSTAT  1 0 -1 0  DRUG*ACTIVITY*SSSTAT  1 0 -1 0 0 0 0 0,
	ACTIVITY  1 -1  DRUG*ACTIVITY  0 0 1 -1   
	ACTIVITY*SSSTAT  1 0 -1 0  DRUG*ACTIVITY*SSSTAT  0 0 0 0 1 0 -1 0,
	ACTIVITY  1 -1  DRUG*ACTIVITY  1 -1 0 0   
	ACTIVITY*SSSTAT  0 1 0 -1  DRUG*ACTIVITY*SSSTAT  0 1 0 -1 0 0 0 0,
	ACTIVITY  1 -1  DRUG*ACTIVITY  0 0 1 -1  
	ACTIVITY*SSSTAT  0 1 0 -1  DRUG*ACTIVITY*SSSTAT  0 0 0 0 0 1 0 -1  ;
ESTIMATE  'Item 1d: b-c from above and  TIME  1 0 -1  SSSTAT*TIME  1 0 -1 0 0 0'
	TIME 1 0 -1  SSSTAT*TIME  1 0 -1 0 0 0  DRUG 1 -1   ACTIVITY -1 1  DRUG*ACTIVITY  0 1 -1 0
	DRUG*ACTIVITY*SSSTAT   0 0 1 0 -1 0 0 0   DRUG*SSSTAT  1 0 -1 0
	ACTIVITY*SSSTAT  -1 0 1 0;
ESTIMATE  'Item 1e: h-e from above and  TIME  0 1 -1  SSSTAT*TIME  0 0 0 0 1 -1'
	TIME 0 1 -1  SSSTAT*TIME  0 0 0 0 1 -1  DRUG -1 1   ACTIVITY -1 1  DRUG*ACTIVITY  -1 0 0 1
	DRUG*ACTIVITY*SSSTAT   0 -1 0 0 0 0 0 1   DRUG*SSSTAT  0 -1 0 1  
	ACTIVITY*SSSTAT 0 -1 0 1;
RUN;
  
PROC MIXED DATA = DSST method = ML covtest;
CLASS SUBJECT DRUG ACTIVITY SSSTAT TIME;
MODEL CORRECTTRIALS = DRUG   TIME SSSTAT
		      SSSTAT*DRUG   SSSTAT*TIME CORRECTTRIALS_B
			   / SOLUTION;
RANDOM SUBJECT; 
RUN;

PROC IMPORT DATAFILE = 'C:\Documents and Settings\richc\My Documents\CPH931F08\Mercury.xls'
     OUT = Mercury  DBMS = EXCEL REPLACE;
     SHEET = Mercury;
     GETNAMES = YES;
RUN;

DATA MERCURY;
SET  MERCURY;
Y = LOG( __YR_STANDARD_MERCURY);
X1 = PH;
X2 = LOG(CHLOROPHYLL);
Z1 = LOG(AVG_MERCURY);
Z2 = LOG(MIN);
Z3 = LOG(MAX);
RUN;

PROC REG DATA = Mercury;
MODEL Y = X1 X2;
RUN;

PROC MI DATA=MERCURY SEED=12345 nimpute=5 out=imputed;
	mcmc chain=multiple displayinit initial=em(itprint);
var Y X1 X2;
RUN;
 
PROC PRINT DATA=imputed;
var Y X1 X2;
RUN;
 
PROC REG DATA=imputed outest=outreg covout;
MODEL Y = X1 X2;
BY _IMPUTATION_;
RUN;

PROC MIANALYZE DATA=OUTREG;
VAR INTERCEPT X1 X2;
RUN;
 

PROC MI DATA=MERCURY SEED=12345 nimpute=5 out=imputed;
	mcmc chain=multiple displayinit initial=em(itprint);
var Y X1 X2 Z1 Z2 Z3;
RUN;
 
PROC PRINT DATA=imputed;
var Y X1 X2 Z1 Z2 Z3;
RUN;
 
PROC REG DATA=imputed outest=outreg covout;
MODEL Y = X1 X2;
BY _IMPUTATION_;
RUN;

PROC MIANALYZE DATA=OUTREG;
VAR INTERCEPT X1 X2;
RUN;
 

ODS RTF CLOSE;
RUN;