data Landscape; infile "Root injury, abundance ANOVA and correlations DATA.csv" dlm=',' firstobs=2; input Field $ Year Location Cat Dis Root WCR; sqrt_root=sqrt(Root); sqrt_WCR=sqrt(WCR); run; proc print; run; proc sort; by cat; run; proc means mean stderr; var root_a WCR; by Cat; run; /*Adult abundance, full model*/ proc mixed; class year location cat; model sqrt_WCR=cat year cat*year; random location location*year location*cat location*year*cat; run; proc mixed; class year location cat; model sqrt_WCR=cat year cat*year; random location location*year location*cat; run; proc mixed; class year location cat; model sqrt_WCR=cat year cat*year; random location location*year; run; /*FINAL model for WCR abundance*/ proc mixed; class year location cat; model sqrt_WCR=cat year cat*year; random location; run; proc mixed; class year location cat; model sqrt_wcr=cat year cat*year; run; /*Model selection for Root injury*/ proc mixed; class year location cat; model sqrt_root=cat year cat*year; random location location*year location*cat location*year*cat; run; proc mixed; class year location cat; model sqrt_root=cat year cat*year; random location location*year location*cat; run; proc mixed; class year location cat; model sqrt_root=cat year cat*year; random location location*year; run; proc mixed; class year location cat; model sqrt_root=cat year cat*year; random location; run; /*FINAL MODEL for root injury --> No random terms retained*/ proc mixed; class year location cat; model sqrt_root=cat year cat*year; run; proc corr pearson; var root wcr dis; run;