Question

我正在尝试计算双侧Kolmogorov-Smirnov检验的临界值（ PROC NPAR1WAY 不输出这些！）。计算公式为c（a）* sqrt（（n + m）/（nm））其中n和m是每个数据集中的观测数，c（a）= 1.36，置信水平a = 0.05。

无论哪种，

A）SAS中有一个例程可以为我计算这些吗？（我已经搜索了一段时间）或，

B）自己计算统计数据的最佳方法是什么？我最初的方法是从每个数据集中选择行数到宏变量然后计算统计量，但这感觉很难看。

提前致谢

Answer 1

A）如果您搜索了所有相关文档，可能不会。

B）该方法听起来不错，但如果您愿意，可以使用数据步骤，例如

data example1 example2;
    set sashelp.class;
    if _n_ < 6 then output example1;
    else output example2;
run;

data _null_;
    if 0 then set example1 nobs = n;
    if 0 then set example2 nobs = m;
    call symput('Kolmogorov_Smirnov_05',1.36 * sqrt((n+m)/(n*m)));
run;

%put &=Kolmogorov_Smirnov_05;

SAS - Kolmogorov-Smirnov双边临界值

1 个答案: