各位大神，我从网上找了个卡方两两比较的SAS程序，但是只能比较到2和3，我的是7个组，求帮忙看一下程序哪里有问题，谢谢！
/* libname:The libname of your dataset */
/* filename:The name of your dataset */
/* fvar: The WEIGHT variable */
/* gvar: The variable for grouping */
/* outcomevar: The outcome variable*/
%macro multicompare(libname=,filename=,fvar=,gvar=,outcomevar=);
options nodate nonumber formdlim='*';
title;
/*Conduct Chi-sqare test for comparison of all groups*/
proc freq data=test;
tables &gvar*&outcomevar / nocol norow nopercent chisq;
weight f;
run;
/*Create symbols to distinguish different subgroups*/
proc sql noprint;
select n(distinct &gvar) into :gn from &libname..&filename.;
select distinct a.&gvar,b.&gvar into :ag separated by '-',:bg separated by '-'
from &libname..&filename. a,&libname..&filename. b where a.&gvar<b.&gvar;
quit;
%let gn1=%eval(&gn);
/*Combine subgroups*/
%Macro chisq;
%do t=1 %to &gn1;
%let m=%scan(&ag,&t);
%let n=%scan(&bg,&t);
proc sql ;
create table test as select * from &libname..&filename. where &gvar=&m or &gvar=&n;
quit;
/*Conduct Chi-sqare test for multicomparison of subgroups*/
proc freq data=test;
tables &gvar*&outcomevar / nocol norow nopercent chisq;
weight f;
run;
/*Calculate the adjusted alpha*/
data alpha;
alpha=0.05/(&gn1*(&gn1-1)/2+1);
proc print data=alpha;
run;
%end;
%mend chisq;
%chisq;
%mend multicompare;
data a;
do group=1 to 7;
do factor=1 to 2;
input f @@;output;
end;
end;
cards;
6 36 9 11 17 18 19 31 15 6 7 7 9 14
;
run;
%multicompare(libname=work,filename=a,fvar=f,gvar=group,outcomevar=factor);
