data bioassaysurvivalcounty;
infile "Bioassay - survival data.csv" dlm=',' firstobs=2;
input pop $ fieldtype $ year corn $ survival infested;
propsurv=(survival/infested);
arsinpropsurv=arsin(propsurv);
asr_propsurv=arsin(sqrt(propsurv));
run;
proc print;
run;

proc means mean stderr;
class fieldtype corn;
var propsurv;
run;


proc sort;
by year;
run;

proc glm;
class corn pop;
model asr_propsurv=pop|corn;
where year=2015;
lsmeans pop pop*corn;
contrast "PC5 vs controls (BT)" pop -1 1 pop*corn -1 1 0 0; 
contrast "PC5 vs controls (ISO)" pop -1 1 pop*corn 0 0 -1 1;
run;

proc glm;
class corn pop year;
model asr_propsurv=pop|corn;
where year=2016;
lsmeans pop pop*corn;
contrast	"2016 - PC11 vs controls (Bt)"	pop	0	2	0	0	-1	-1	pop*corn	0	2	0	0	-1	-1	0	0	0	0	0	0	;
contrast	"2016 - NPC6 vs controls (Bt)"	pop	0	0	2	0	-1	-1	pop*corn	0	0	2	0	-1	-1	0	0	0	0	0	0	;
contrast	"2016 - PC7 vs controls (Bt)"	pop	0	0	0	2	-1	-1	pop*corn	0	0	0	2	-1	-1	0	0	0	0	0	0	;
contrast	"2016 - PC11 vs controls (Iso)"	pop	0	2	0	0	-1	-1	pop*corn	0	0	0	0	0	0	0	2	0	0	-1	-1	;
contrast	"2016 - NPC6 vs controls (Iso)"	pop	0	0	2	0	-1	-1	pop*corn	0	0	0	0	0	0	0	0	2	0	-1	-1	;
contrast	"2016 - PC7 vs controls (Iso)"	pop	0	0	0	2	-1	-1	pop*corn	0	0	0	0	0	0	0	0	0	2	-1	-1	;
run;

proc glm;
class corn pop year;
model asr_propsurv=pop|corn;
where year=2017;
lsmeans pop pop*corn;
contrast	"2017 - NPC14 vs controls (Bt)"	pop	4	0	0	0	0	0	-1	-1	-1	-1	pop*corn	4	0	0	0	0	0	-1	-1	-1	-1	0	0	0	0	0	0	0	0	0	0	;
contrast	"2017 - PC19 vs controls (Bt)"	pop	0	4	0	0	0	0	-1	-1	-1	-1	pop*corn	0	4	0	0	0	0	-1	-1	-1	-1	0	0	0	0	0	0	0	0	0	0	;
contrast	"2017 - PC20 vs controls (Bt)"	pop	0	0	4	0	0	0	-1	-1	-1	-1	pop*corn	0	0	4	0	0	0	-1	-1	-1	-1	0	0	0	0	0	0	0	0	0	0	;
contrast	"2017 - NPC11 vs controls (Bt)"	pop	0	0	0	4	0	0	-1	-1	-1	-1	pop*corn	0	0	0	4	0	0	-1	-1	-1	-1	0	0	0	0	0	0	0	0	0	0	;
contrast	"2017 - NPC12 vs controls (Bt)"	pop	0	0	0	0	4	0	-1	-1	-1	-1	pop*corn	0	0	0	0	4	0	-1	-1	-1	-1	0	0	0	0	0	0	0	0	0	0	;
contrast	"2017 - NPC13 vs controls (Bt)"	pop	0	0	0	0	0	4	-1	-1	-1	-1	pop*corn	0	0	0	0	0	4	-1	-1	-1	-1	0	0	0	0	0	0	0	0	0	0	;
contrast	"2017 - NPC14 vs controls (Iso)"	pop	4	0	0	0	0	0	-1	-1	-1	-1	pop*corn	0	0	0	0	0	0	0	0	0	0	4	0	0	0	0	0	-1	-1	-1	-1	;
contrast	"2017 - PC19 vs controls (Iso)"	pop	0	4	0	0	0	0	-1	-1	-1	-1	pop*corn	0	0	0	0	0	0	0	0	0	0	0	4	0	0	0	0	-1	-1	-1	-1	;
contrast	"2017 - PC20 vs controls (Iso)"	pop	0	0	4	0	0	0	-1	-1	-1	-1	pop*corn	0	0	0	0	0	0	0	0	0	0	0	0	4	0	0	0	-1	-1	-1	-1	;
contrast	"2017 - NPC11 vs controls (Iso)"	pop	0	0	0	4	0	0	-1	-1	-1	-1	pop*corn	0	0	0	0	0	0	0	0	0	0	0	0	0	4	0	0	-1	-1	-1	-1	;
contrast	"2017 - NPC12 vs controls (Iso)"	pop	0	0	0	0	4	0	-1	-1	-1	-1	pop*corn	0	0	0	0	0	0	0	0	0	0	0	0	0	0	4	0	-1	-1	-1	-1	;
contrast	"2017 - NPC13 vs controls (Iso)"	pop	0	0	0	0	0	4	-1	-1	-1	-1	pop*corn	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	4	-1	-1	-1	-1	;

run;




/** within population difference in survival and development on Cry3Bb1 vs isoline**/
proc sort data=bioassaysurvivalcounty;
by pop;
run;
proc ttest data=bioassaysurvivalcounty sides=L;
class corn;
var asr_propsurv;
by pop;
run;