close

當兩組樣本資料屬性為連續變項,相互配對或重覆且

(1) 其差距服從常態分佈時,則採Paired t test檢定兩組差距之平均值是否相等;

(2) 其差距未服從常態分佈時,則採Wilcoxon Signed Rank Test檢定兩組差距之中位數是否相等。

 

目標: 檢測before & after 是否有顯著差異

Step 1: 產生10個1~50隨機配對資料;

data sample;
*產生10個1~50隨機整數;
do i=1 to 10;
  before=int(50*ranuni(123));
  after=int(50*ranuni(123));
  output;
end;
run;
proc print;run;
 
data a;
do i=1 to 10;
  before=int(50*ranuni(123));
  after=int(50*ranuni(123));
  output;
end;
run;

 

Step 2: 將前後資料相減

 
data a1;set a;
  diff=after-before;
proc print;run;

 

兩組配對資料檢定-Wilcoxon Signed Rank

 

Step 3. 用 PROC UNIVARIATE 語法擷取相關檢定

 
proc univariate normal data=a1; 
var diff;
run;

 

報表解讀

兩組配對資料檢定-Wilcoxon Signed Rank

樣本數小,採Shapiro-Wilk檢定,p=0.3684 >0.05,通過常態檢定,採Paired t test:

 

兩組配對資料檢定-Wilcoxon Signed Rank

Paired t test 的目的是比較前後兩組的差距是否為0,因此可以前後相減,再檢驗此差距的常態性,

(1) 若通過常態分佈檢定,則檢驗此差距平均值是否為0

在SAS Proc Univariate 中,可由位置檢定中的 Student's t 得到該差距數據的統計量t為-0.18077,p=0.8606

即,該數據前後沒有統計上的顯著差異(p=0.8606)

 

(2) 若該差距未通過常態檢定,則檢驗此差距中位數是否為0

可由位置檢定中的 符號秩 得到該差距數據的統計量S為-3,p=0.7910

 

參考:

https://communities.sas.com/t5/Base-SAS-Programming/wilcoxon-signed-rank-test-significance/td-p/234753
 
.....
Step 1 & Step 2 可以再簡化~試試看?
 
arrow
arrow
    創作者介紹
    創作者 Vivian 的頭像
    Vivian

    統計散記-Vivian 經驗分享

    Vivian 發表在 痞客邦 留言(0) 人氣()