data LandscapeCorrectedSurvival;
infile "Corrected survival DATA.csv" dlm=',' firstobs=2;
input Pop $ Type $ cat $ Year Corrected;
run;
proc print;
run;
proc sort;
by year;
run;

/* Corrected survival contrasts with controls for 2015 */
proc glm;
class year pop;
model corrected=pop;
lsmeans pop;
contrast "1c vs controls" pop 20 0 0 0 -5 -5 0 -5 -5;
contrast "1d vs controls" pop 0 20 0 0 -5 -5 0 -5 -5;
contrast "1 vs controls" pop 0 0 20 0 -5 -5 0 -5 -5;
contrast "5a vs controls" pop 0 0 0 20 -5 -5 0 -5 -5;
contrast "2 vs controls" pop 0 0 0 0 -5 -5 20 -5 -5;
where year=2015;
run;

/* Corrected survival contrasts with controls for 2016 */
proc glm;
class year pop;
model corrected=pop;
lsmeans pop;
contrast "10a vs controls" pop 25 0 0 0 0 -5 -5 -5 -5 -5 0;
contrast "9 vs controls" pop 0 25 0 0 0 -5 -5 -5 -5 -5 0;
contrast "9a vs controls" pop 0 0 25 0 0 -5 -5 -5 -5 -5 0;
contrast "11 vs controls" pop 0 0 0 25 0 -5 -5 -5 -5 -5 0;
contrast "11a vs controls" pop 0 0 0 0 25 -5 -5 -5 -5 -5 0;
contrast "9b vs controls" pop 0 0 0 0 0 -5 -5 -5 -5 -5 25;
where year=2016;
run;

/* Corrected survival contrasts with controls for 2017 */
proc glm;
class year pop;
model corrected=pop;
lsmeans pop;
contrast "20b vs controls" pop 20 0 0 -5 -5 -5 -5;
contrast "14c vs controls" pop 0 20 0 -5 -5 -5 -5;
contrast "16a vs controls" pop 0 0 20 -5 -5 -5 -5;
where year=2017;
run;

/* No difference between FF and SURR (sig effect of year but not type or interaction) */
proc sort;
by year;
run;

Proc mixed;
class type year;
model corrected=year type year*type;
where type not eq "CONTROL";
run;

proc mixed;
class type;
model corrected=type;
where type not eq "CONTROL";
run;