var rchl; 
run; 
*( step to get rid of heterogeneous variances; 
Proc Sort data=reschl; 
by SegGrp year; 
run; 
Proc Means data=reschl noprint; 
by SegGrp year; 
var rchl; 
output out=sdchl StdDev = sdchl n=n; 
run; 
*(proc print data=sdchl; 
*( title "standard deviation results"; 
*( var SegGrp year sdchl n; 
data ResChl; 
merge ResChl sdchl; 
by SegGrp year; 
StdResChl = rchl/sdchl; 
run; 
proc Univariate normal plot data=work.reschl; 
title2 "Normality test on standardized residuals"; 
var StdResChl; 
run; 
* (proc contents data=work.reschl; 
*(run; 
data four; 
set work.reschl; 
aRChl= abs(RChl); 
aStdResChl= abs(StdResChl); 
run; 
proc glm data=work.four; 
title2 "Levene's test for both types of residuals"; 
class year cbseg_2003; 
model aRChl aStdResChl=year cbseg_2003; 
*(means cbseg_2003/snk; 
*(lsmeans cbseg_2003; 
run; 
★ ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 
* PROGRAM: JAMES_RIVER.sas 
* This program will TEST CHLOROPHYLL DATA FOR NORMALITY 
* Base code from Elgin Perry 02/16/2010 
* additional code written by Jackie Johnson 02/17/2010 
* ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ • 
t 
libname ALGAE "G:\LR\OTHER_LR_DATA\Criteria_work\chlorophyll\2010" 
*libname ALGAE "C:\Projects\CBP\CHLCRIT\LogNormal\"; 
options Is=72; 
♦OPTIONS LS=120 PS=55 REPLACE NOCENTER; 
OPTIONS formchar = '|-|+|- +=|-/\<>*'; 
♦PROC IMPORT 0UT= ALGAE.JAMES_SPRING_CHL 
DATATABLE= "JAMES_SPRING_CHL" 
DBMS=ACCESS20 0 0 REPLACE; 
