data work1;
infile 'e:\trapping\trapping1_10.txt' delimiter='09'x;
input breed$ id size party sniff rush;
proc genmod data=work1 descending;
class breed id sniff crush;
model crush=breed size party sniff breed*sniff/dist=bin link=log type3;
repeated subject=id(breed)/type=ind covb corrw;
contrast 'breed=pm sniff=0 vs pm' sniff 0/estimate=parm;
contrast 'breed=pm sniff=1 vs pm' sniff 1/estimate=parm;
contrast 'breed=py sniff=0 vs py' sniff 1 breed*sniff 1/estimate=parm;
contrast 'breed=py sniff=1 vs py' sniff 0 1 breed*sniff 0 1/estimate=parm;
output out=work1_1 pred=breed;
run;
1、想输出自变量breed和sniff的预测概率值？
2、想进行breed*sniff互作效应间的两两比较？