options ls=78 ps=55;
title 'Bartlett Root Chemistry Analyses';
proc import datafile='c:\stats_data\RootChem\BRootChem.xls' out=root replace;

proc mixed;
class site pit sdepth rdia;
model ca = sdepth rdia sdepth*rdia;
random site pit(site);
lsmeans sdepth*rdia /e slice=rdia;
run;
proc mixed;
class site pit sdepth rdia;
model mg = sdepth rdia sdepth*rdia;
random site pit(site);
lsmeans sdepth*rdia /e slice=rdia;
run;
proc mixed;
class site pit sdepth rdia;
model k = sdepth rdia sdepth*rdia;
random site pit(site);
lsmeans sdepth*rdia /e slice=rdia;
run;
proc mixed;
class site pit sdepth rdia;
model p = sdepth rdia sdepth*rdia;
random site pit(site);
lsmeans sdepth*rdia /e slice=rdia;
run;


options ls=78 ps=55;
title 'Bartlett Root Chemistry Analyses';

proc import datafile='c:\stats_data\RootChem\BRootChem.xls' out=root replace;
proc mixed;
class site pit sdepth rdia;
model ca = sdepth rdia ;
random site pit(site);
lsmeans rdia /e pdiff;
run;

proc mixed;
class site pit sdepth rdia;
model k = sdepth rdia ;
random site pit(site);
lsmeans rdia /e pdiff;
run;

proc mixed;
class site pit sdepth rdia;
model mg = sdepth rdia ;
random site pit(site);
lsmeans rdia /e pdiff;
run;

proc mixed;
class site pit sdepth rdia;
model p = sdepth rdia ;
random site pit(site);
lsmeans rdia /e pdiff;
run;


options ls=78 ps=55;
title 'Bartlett Root Chemistry Analyses';
proc import datafile='c:\stats_data\RootChem\BRootChem.xls' out=root replace;

proc mixed;
class age pit sdepth rdia site;
model ca = age sdepth rdia age*sdepth age*rdia age*sdepth*rdia;
random site(age) pit(site);
run;

proc mixed;
class age pit sdepth rdia site;
model k = age sdepth rdia age*sdepth age*rdia age*sdepth*rdia;
random site(age) pit(site);
run;

proc mixed;
class age pit sdepth rdia site;
model mg = age sdepth rdia age*sdepth age*rdia age*sdepth*rdia;
random site(age) pit(site);
run;

proc mixed;
class age pit sdepth rdia site;
model p = age sdepth rdia age*sdepth age*rdia age*sdepth*rdia;
random site(age) pit(site);
run;