kopid Kop Stock Used bovid Bov Stock Used pwin Prediction of winner based on bacteriocin assay nemerged number of host that emerged nid number of hosts that were successfully id'd to species nbov number of hosts id'd as bov nkop number of hosts id'd as kop kophr hazard ratio of kop stock_all are relative to Kop43 bovhr hazard ratio of bov stock_all are relative to Kop43 difrel43 kophr-bovhr kgeno genotype of kop stock used bgeno genotype of bov stock used The SAS code used in the analysis: title 'Test of Predictions'; proc glimmix data=sum; class kopid bovid pwin kgeno bgeno; model nkop/nid=pwin|difrel43 /solution link=logit ; *lsmeans pwin /pdiff cl ilink at difrel43=-5; *lsmeans pwin /pdiff cl ilink at difrel43=-3; *lsmeans pwin /pdiff cl ilink at difrel43=0; lsmeans pwin /pdiff cl ilink ; *lsmeans pwin /pdiff cl ilink at difrel43=3; *lsmeans pwin /pdiff cl ilink at difrel43=5; random kgeno bgeno; run; proc glimmix data=sum ic=pq ; where pwin='either'; class kopid bovid pwin kgeno bgeno; model nkop/nid= difrel43 /solution link=logit ; random kgeno bgeno; run; title 'contrasts with just pwin'; proc glimmix data=sum; class kopid bovid pwin kgeno bgeno; model nkop/nid=pwin /solution link=logit ; contrast "Bov v. KOp" pwin -1 0 1; contrast "bov v. either" pwin -1 1 0; contrast "either v. kop" pwin 0 -1 1; contrast "bov/eith v. kop" pwin -.5 -.5 1; random kgeno bgeno; run; title 'Bov vs. Either'; proc glimmix data=sum; where pwin ne 'kop'; class kopid bovid pwin kgeno bgeno; model nkop/nid=pwin|difrel43/solution link=logit ; lsmeans pwin /pdiff cl ilink at difrel43=-5; lsmeans pwin /pdiff cl ilink at difrel43=-3; lsmeans pwin /pdiff cl ilink at difrel43=0; lsmeans pwin /pdiff cl ilink at difrel43=3; lsmeans pwin /pdiff cl ilink at difrel43=5; random kgeno bgeno; run;