****
proc mixed data=root;
class site pit sdepth rdia;
model ca =site pit sdepth rdia sdepth*rdia/solution;
random site pit(site);
store  sasuser.rootanalysis;
run;
proc plm source=sasuser.rootanalysis;
run;
proc plm source=sasuser.rootanalysis;
   show fit parms;
   test sdepth rdia sdepth*rdia;
run;
*****
proc mixed data=root;
class site pit sdepth rdia;
model ca =site pit sdepth rdia sdepth*rdia/solution;
random site pit(site);
store  sasuser.rootanalysis / label='PLM: Getting Started';
run;
proc plm restore=sasuser.rootanalysis;
score data=scoreX out=pred;
run;
****
proc mixed data=root;
class site pit sdepth rdia;
model ca =site pit sdepth rdia sdepth*rdia/solution;
random pit(site);
store  sasuser.rootanalysis / label='PLM: Getting Started';
run;
proc plm restore=sasuser.rootanalysis;
score data=rootnew out=out / ilink;
run; 
***
proc plm restore=sasuser.rootanalysis;
score data=root out=rootmodel predicted/ilink;
run;
ods graphics on;
proc plm restore=sasuser.rootanalysis;
effectplot;
run;
ods graphics off;