data county;
infile "Root injury and abundance data.csv" dlm=',' firstobs=2;
input Field $ Farm Year Cat $ Root WCR corrected;
logroot=log(1+root);
inverseroot=(1/root); 
logwcr=log(1+wcr); 
sqrt_wcr=sqrt(wcr);
sqrt_root=sqrt(root);
run;
proc print;
run;
/*MEANS FOR GRAPHS*/
proc sort;
by Year Cat;
run;
proc univariate;
var root;
by Year Cat;
run;
proc univariate;
var WCR;
by Year Cat;
run;
proc sort;
by Cat;
run;
proc univariate;
var Root;
by Cat;
run;
proc univariate;
var WCR;
by Cat;
run;

proc glm;
class Year Cat;
model sqrt_wcr = Year Cat Year*Cat / solution;
run;
proc glm;
class Cat;
model sqrt_wcr = Cat / solution;
run;

proc glm;
class Year Cat;
model sqrt_root = Year Cat Year*Cat / solution;
lsmeans year*cat/ pdiff adjust=bon lines;
run;

proc mixed method=type3;
class Cat;
model logboth = Cat / residual solution;
run;

/* CORRECTED SURVIVAL */
proc glm;
class cat;
model corrected = cat;
run;