Wilcoxon 符號秩檢定當我在 JMP 中執行 Wilcoxon Signed-Rank 選項時,我得到的答案與 SAS 中得到的答案截然不同。 有沒有辦法取得每個選項(SAS 和 JMP)中使用的程式碼或用於找出結果不同的公式?
我在這裡先向您的幫助表示感謝。回覆:Wilcoxon符號秩檢定如上所述,
大家好 - 我已在下面設定了一些 SAS 程式碼,按照前面參考文獻中的描述進行分析,然後透過 JMP 執行相同的分析...
要在 JMP 中執行相同操作,請下載 James Abbey 在www.personal.psu.edu/jda188/JMP/Signedrank.doc上的以下說明,然後按照那裡的說明對資料集的 DIFF 列進行操作。
也就是說,將 FERTILIZER1.SAS7BDAT 或 FERTILIZER3.SAS7BDAT 從 WILC 庫匯入 JMP,然後使用 Analyze | 分析 DIFF 列分配。將 DIFF 列複製到「Y,列」佔位符中,然後按一下「確定」。接下來,點擊第二個紅色三角形並選擇「檢定平均值」。在「指定假設平均值」中輸入 0,並勾選「Wilcoxon 符號秩」方塊,然後按一下「確定」。在輸出的「測試平均值=值」部分,您應該找到測試統計量= 74,以及 P 值 0.0017 - 這與 SAS 產生的答案相同。
可以嗎?
回覆:Wilcoxon符號秩檢定很有可能其中一個只是設定不正確,所以最好的方法是先找一個教科書上的例子,然後用 SAS 和 JMP 分析這些數據,找出哪個方法沒有產生預期的結果。我個人最喜歡的、簡單易懂的例子是 Sidney Siegel 的《行為科學的非參數統計》,但還有很多其他的例子。
我的第二個想法是,你執行的 SAS 或 JMP 檢定可能實際上不是配對樣本的符號秩檢定,而是獨立樣本的秩和檢定(有時也稱為 Mann-Whitney U 檢定)。這兩個檢定分別是配對 t 檢定和非配對 t 檢定的非參數類似物,並且都與 Wilcoxon 的名字有關:兩者的詳細資訊以及範例都可以在維基百科中找到,連結如下:
http://en.wikipedia.org/wiki/Wilcoxon_signed-rank_test
http://en.wikipedia.org/wiki/Mann-Whitney-Wilcoxon_test回覆:威爾科克森符號秩檢驗以下是加州大學洛杉磯分校統計學術技術服務中心的範例數據,該數據與上述觀點相符:
回覆:威爾科克森符號秩檢驗感謝二位的回饋。
David Q/RI 用的是書中的例子,所以我才對得到不同的答案感到驚訝。
有沒有辦法手動運行原始程式碼來解決這個問題,並找出程式碼在哪一步可能出現不同的行為?
再次感謝。回覆:威爾科克森符號秩檢驗如果有什麼方法可以獲取實際的源代碼,恐怕我也不知道是什麼——但我剛才快速瀏覽了一下MattF在上面帖子中提到的Wilcoxon檢驗,分別在SAS和JMP中查看了一下,它們看起來都像是對數據對之間的差異進行Mann-Whitney-Wilcoxon檢驗,這意味著你實際上需要先自己計算出數據對之間的差異,然後再比較這些差異。你實際上是這樣做的嗎?
以下提供另一個參考資料,展示瞭如何在 SAS 中執行 Wilcoxon 配對符號秩檢驗,這表明基本檢定實際上並不是以簡單程序形式提供的(否則為什麼要費這麼大勁,用這種稍微繞彎的方式呢?):
http://udel.edu/~mcdonald/statsignedrank.html
如果你願意貼出你正在運行的程式碼副本,我很樂意檢查我的判斷是否正確。