在SAS Proc Freq中使用数组批量执行Fisher精确检验的问题咨询

你好呀！看起来你踩了SAS里一个常见的小坑——Proc Freq是无法直接识别DATA步中定义的数组的，所以你直接写comp_vars(22)才会抛出那两条语法错误，因为Proc Freq根本不知道这个数组指代的是哪些变量~

针对你要批量跑22个变量和switch_var的Fisher精确检验的需求，给你两个实用的解决方案：

方案一：如果变量名有规律（比如comp_var1到comp_var22）——用宏循环 #

这是最简洁的方式，通过宏循环自动生成每一组检验的语句：

%macro run_fisher_tests();
    proc freq data=dataset;
        %do i=1 %to 22;
            tables switch_var*comp_var&i / fisher;
        %end;
    run;
%mend run_fisher_tests;

/* 调用宏执行所有检验 */
%run_fisher_tests();

这个宏会循环22次，每次生成一个switch_var*comp_varX的交叉表语句，自动完成所有组合的Fisher检验。

方案二：如果变量名无规律——先把数组变量名转成宏变量列表 #

如果你已经在DATA步里定义了数组comp_vars，可以先把数组里的变量名提取出来拼成一个宏变量，再传给Proc Freq：

/* 第一步：把数组中的变量名提取到宏变量里 */
data _null_;
    set dataset;
    /* 这里要和你定义的数组类型匹配，数值型就去掉$ */
    array comp_vars[22] $;
    /* 定义足够长的变量来存所有变量名，22个变量的话1000长度完全够用 */
    length var_list $ 1000;
    do i=1 to dim(comp_vars);
        var_list = catx(' ', var_list, vname(comp_vars[i]));
    end;
    /* 把拼接好的变量列表存入宏变量comp_var_list */
    call symput('comp_var_list', var_list);
    stop; /* 只需要读取第一行数据即可，不用遍历全表 */
run;

/* 第二步：用宏变量批量执行检验 */
proc freq data=dataset;
    tables switch_var*(&comp_var_list) / fisher;
run;

这个方法先把数组里的所有变量名转成空格分隔的字符串，Proc Freq就能识别这个变量列表，自动生成所有交叉表并执行Fisher检验。如果你的变量是数值型的，记得把数组定义里的$去掉哦~

备注：内容来源于stack exchange，提问作者Ryan