data countyfield history;
infile "Field history analysis data.csv" dlm="," firstobs=2;
input field $ farm year corn nonbt cry3 cry3435 pyramid FH1 FH2 FH3 FH4 FH5 FH6 FH7 FH8 FH9 FH10 FH11 FH12 cat root wcr ncr both;
FH1a=log(FH1);
FH2b=arsin(sqrt(FH2));
FH5b=arsin(sqrt(FH5));
FH6b=arsin(sqrt(FH6));
FH7b=arsin(sqrt(FH7));
FH8b=arsin(sqrt(FH8));
FH9b=arsin(sqrt(FH9));
FH10b=arsin(sqrt(FH10));
FH11b=arsin(sqrt(FH11));
FH12b=arsin(sqrt(FH12));
sqrt_root=sqrt(root);
sqrt_wcr=sqrt(wcr);

run;
proc print;
run;

proc corr;
var nonbt cry3 cry3435 pyramid FH4 FH1 FH2 FH5 FH6 FH7 FH8 FH10 FH11 FH12;
run;

proc freq;
tables NonBt*cat /chisq;
tables Cry3*cat / chisq;
tables Cry3435*cat / chisq;
tables Pyramid*cat / chisq;
tables FH4*cat / chisq;
run;

proc sort;
by cat year;
run;
ods output summary=with_stackods;
proc means stackodsoutput n mean stderr;
var FH1 FH2 FH5 FH6 FH7 FH8 FH10 FH11 FH12;
by cat;
run; 

/*total years continuous corn*/
proc ttest sides=L;
	class Cat;
	var FH1a;
	run;
/*Proportion corn planted in 6 yr*/
proc ttest sides=L;
	class Cat;
	var FH2b;
	run;
/* Proportion S.I. use in 6 years */
proc ttest sides=L;
	class Cat;
	var FH5b;
	run;
/* Proportion non Bt planted */
proc ttest sides=U;
	class Cat;
	var FH6b;
	run;
/* Proportion Cry3Bb1 or mCry3a */
proc ttest;
	class Cat;
	var FH7b;
	run;
/* Proportion Cry34/35Ab1 */
proc ttest;
	class Cat;
	var FH8b;
	run;
/*Prop pyramided corn*/
proc ttest;
	class Cat;
	var FH10;
	run;
/*Prop SI on non-BT*/
proc ttest;
	class Cat;
	var FH12;
	run;
/*Prop Bt+S.I.*/
proc ttest;
	class cat;
	var FH11b;
	run;



/* DID ANY FIELD HISTORY METRICS AFFECT ROOT INJURY OR WCR ABUNDANCE? */

proc reg plots (only)=ALL;															/*USE THIS ONE*/
	Model sqrt_wcr= nonbt cry3435 pyramid FH1a FH2 FH5 FH6 FH7 FH8 FH10 FH12/ selection=stepwise
	SLE=0.25 SLS=0.15 include=0;
	run;

proc reg plots (only)=ALL;															/*USE THIS ONE*/
	Model sqrt_root= nonbt cry3435 pyramid FH1a FH2 FH5 FH6 FH7 FH8 FH10 FH12  / selection=stepwise
	SLE=0.25 SLS=0.15 include=0;
	run;