data landscapebioassay1; infile "Bioassay survival and third instar DATA.csv" dlm=',' firstobs=2; input pop $ fieldtype $ year corn $ propsurv propthird; arsinpropsurv=arsin(propsurv); arsinpropthird=arsin(propthird); ASR_propsurv=arsin(sqrt(propsurv)); ASR_propthird=arsin(sqrt(propthird)); run; proc print; run; proc means n mean stderr; var propsurv propthird; where fieldtype="Control" and corn="Bt"; run; /** MEANS for survival and development **/ proc sort; by fieldtype year; run; ods output summary=with_stackods; proc means data=landscapebioassay1 stackodsoutput n mean stderr; class pop year corn; var propsurv; by fieldtype year; run; proc sort data=landscapebioassay1; by corn year; run; ods output summary=with_stackods; proc means data=landscapebioassay1 stackodsoutput n mean stderr; class pop year corn; var propthird; by corn year; run; proc glm; class pop corn; model ASR_propsurv=pop corn pop*corn; lsmeans pop pop*corn; where year=2015; contrast "2015 - 1c vs controls(Bt)" pop -4 0 0 0 1 1 0 1 1 pop*corn -4 0 0 0 0 0 0 0 1 0 1 0 0 0 1 0 1 0 ; contrast "2015 - 1c vs controls(Iso)" pop -4 0 0 0 1 1 0 1 1 pop*corn 0 -4 0 0 0 0 0 0 0 1 0 1 0 0 0 1 0 1 ; contrast "2015 - 1d vs controls(Bt)" pop 0 -4 0 0 1 1 0 1 1 pop*corn 0 0 -4 0 0 0 0 0 1 0 1 0 0 0 1 0 1 0 ; contrast "2015 - 1d vs controls(Iso)" pop 0 -4 0 0 1 1 0 1 1 pop*corn 0 0 0 -4 0 0 0 0 0 1 0 1 0 0 0 1 0 1 ; contrast "2015 - 1 vs controls(Bt)" pop 0 0 -4 0 1 1 0 1 1 pop*corn 0 0 0 0 -4 0 0 0 1 0 1 0 0 0 1 0 1 0 ; contrast "2015 - 1 vs controls(Iso)" pop 0 0 -4 0 1 1 0 1 1 pop*corn 0 0 0 0 0 -4 0 0 0 1 0 1 0 0 0 1 0 1 ; contrast "2015 - 4a vs controls(Bt)" pop 0 0 0 -4 1 1 0 1 1 pop*corn 0 0 0 0 0 0 -4 0 1 0 1 0 0 0 1 0 1 0 ; contrast "2015 - 4a vs controls(Iso)" pop 0 0 0 -4 1 1 0 1 1 pop*corn 0 0 0 0 0 0 0 -4 0 1 0 1 0 0 0 1 0 1 ; contrast "2015 - 2 vs controls(Bt)" pop 0 0 0 0 1 1 -4 1 1 pop*corn 0 0 0 0 0 0 0 0 1 0 1 0 -4 0 1 0 1 0 ; contrast "2015 - 2 vs controls(Iso)" pop 0 0 0 0 1 1 -4 1 1 pop*corn 0 0 0 0 0 0 0 0 0 1 0 1 0 -4 0 1 0 1 ; run; proc glm; class pop corn; model ASR_propsurv=pop corn pop*corn; lsmeans pop pop*corn; where year=2016; contrast "2016 - 10a vs controls(Bt)" pop -5 0 0 0 0 1 1 1 1 1 0 pop*corn -5 0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 1 0 1 0 0 0 ; contrast "2016 - 10a vs controls(Iso)" pop -5 0 0 0 0 1 1 1 1 1 0 pop*corn 0 -5 0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 1 0 1 0 0 ; contrast "2016 - 9 vs controls(Bt)" pop 0 -5 0 0 0 1 1 1 1 1 0 pop*corn 0 0 -5 0 0 0 0 0 0 0 1 0 1 0 1 0 1 0 1 0 0 0 ; contrast "2016 - 9 vs controls(Iso)" pop 0 -5 0 0 0 1 1 1 1 1 0 pop*corn 0 0 0 -5 0 0 0 0 0 0 0 1 0 1 0 1 0 1 0 1 0 0 ; contrast "2016 - 9a vs controls(Bt)" pop 0 0 -5 0 0 1 1 1 1 1 0 pop*corn 0 0 0 0 -5 0 0 0 0 0 1 0 1 0 1 0 1 0 1 0 0 0 ; contrast "2016 - 9a vs controls(Iso)" pop 0 0 -5 0 0 1 1 1 1 1 0 pop*corn 0 0 0 0 0 -5 0 0 0 0 0 1 0 1 0 1 0 1 0 1 0 0 ; contrast "2016 - 11 vs controls(Bt)" pop 0 0 0 -5 0 1 1 1 1 1 0 pop*corn 0 0 0 0 0 0 -5 0 0 0 1 0 1 0 1 0 1 0 1 0 0 0 ; contrast "2016 - 11 vs controls(Iso)" pop 0 0 0 -5 0 1 1 1 1 1 0 pop*corn 0 0 0 0 0 0 0 -5 0 0 0 1 0 1 0 1 0 1 0 1 0 0 ; contrast "2016 - 11a vs controls(Bt)" pop 0 0 0 0 -5 1 1 1 1 1 0 pop*corn 0 0 0 0 0 0 0 0 -5 0 1 0 1 0 1 0 1 0 1 0 0 0 ; contrast "2016 - 11a vs controls(Iso)" pop 0 0 0 0 -5 1 1 1 1 1 0 pop*corn 0 0 0 0 0 0 0 0 0 -5 0 1 0 1 0 1 0 1 0 1 0 0 ; contrast "2016 - 9b vs controls(Bt)" pop 0 0 0 0 0 1 1 1 1 1 -5 pop*corn 0 0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 1 0 1 0 -5 0 ; contrast "2016 - 9b vs controls(Iso)" pop 0 0 0 0 0 1 1 1 1 1 -5 pop*corn 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 1 0 1 0 -5 ; run; proc glm; class pop corn; model ASR_propsurv=pop corn pop*corn; lsmeans pop pop*corn; where year=2017; contrast "2017 - 20b vs controls(Bt)" pop -4 0 0 1 1 1 1 pop*corn -4 0 0 0 0 0 1 0 1 0 1 0 1 0 ; contrast "2017 - 20b vs controls(Iso)" pop -4 0 0 1 1 1 1 pop*corn 0 -4 0 0 0 0 0 1 0 1 0 1 0 1 ; contrast "2017 - 14c vs controls(Bt)" pop 0 -4 0 1 1 1 1 pop*corn 0 0 -4 0 0 0 1 0 1 0 1 0 1 0 ; contrast "2017 - 14c vs controls(Iso)" pop 0 -4 0 1 1 1 1 pop*corn 0 0 0 -4 0 0 0 1 0 1 0 1 0 1 ; contrast "2017 - 16a vs controls(Bt)" pop 0 0 -4 1 1 1 1 pop*corn 0 0 0 0 -4 0 1 0 1 0 1 0 1 0 ; contrast "2017 - 16a vs controls(Iso)" pop 0 0 -4 1 1 1 1 pop*corn 0 0 0 0 0 -4 0 1 0 1 0 1 0 1 ; run; /**Ttests for survival and development within each population**/ proc sort data=landscapebioassay1; by pop; run; proc ttest data=landscapebioassay1 sides=L; class corn; var ASR_propsurv; by pop; run; proc ttest data=landscapebioassay1 sides=L; class corn; var ASR_propthird; by pop; where fieldtype="Exp"; run;