2016-10-10 28 views
0

我正在嘗試計算雙側Kolmogorov-Smirnov測試的臨界值(PROC NPAR1WAY不會輸出這些!)。這被計算爲c(a)* sqrt((n + m)/(nm)),其中n和m是每個數據集中觀測值的數量,c(a)= 1.36,置信水平a = 0.05。SAS - Kolmogorov-Smirnov雙面臨界值

要麼,

A)有在SAS常規,將計算這些給我嗎? (我一直在尋找一段時間)或者,

B)自己計算統計量的最好方法是什麼?我最初的做法是將每個數據集中的行數選擇爲宏變量,然後計算統計量,但這種感覺很難看。

在此先感謝

回答

0

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; 
+0

這比我得到的光滑,謝謝。 – VersBersch