編輯問你要SNK兩兩比較的P值，快來看這裡

如何追尋SNK的統計量和P值？

前幾天收到統計聯盟小夥伴的問題，說做了方差分析並進行了SNK兩兩比較，但被雜誌社修回，編輯讓提供Q檢驗的統計量和P值。小夥伴很無奈，說SAS中SNK兩兩比較沒有提供具體的統計量和P值。那麼面對如此挑剔的雜誌社編輯，我們該怎麼辦？有的小夥伴很實在，直接回編輯：「軟體不提供，我也不能提供」。這樣回答估計你的文章就要悲劇了。下面就讓小編教你在SAS中算出兩兩比較的統計量和P值。

1

SAS原方差分析結果

問題來源：作者想比較3組不同手術方式處理組（group）間出血量(blood)的比較？

實施順序：正態性檢驗(滿足)—方差齊性檢驗(滿足)—方差分析—SNK

SAS 代碼：

/*****正態性檢驗**********/

proc univariate data=shun normal;

var blood;

by group;

run;

/****方差齊性檢驗***SNK比較***/

proc anova data=shun;

class group;

model blood=group;

means group/snk hovtest;

run;

quit;

SNK結果展示：

很遺憾確實沒有提供統計量和P值，不能滿足雜誌社要求。

小編最愛的SNK英雄

2

解決辦法

可用GLM模型中CONTRAST語句來進行兩兩比較並提供P值。由於是單因素3水平，contrast的語句構成如下：contrast 「 標籤哪兩組相比」 比較的變數名 （1和-1代表比較組，0不參與比較）。

SAS代碼如下：

proc glm data=shun;

class group;

model blood=group;

means group/hovtest;

means group/snk bon dunnett("1");

contrast "A VS B" GROUP 1 -1 0 ;

contrast "A VS C" GROUP 1 0 -1;

contrast "B VS C" GROUP 0 1 -1;

run;

quit;

結果展示如下：

小編最愛的SNK美女

3

小編寄語

大家可以比對上面的SNK分析結果，是完全一致的，所以這既可以清楚的反映各組兩兩比較的結果，又可以提供相應的P值滿足雜誌社的要求，小夥伴們可以放心使用！

為方便聯盟內部成員交流，特建立微信交流群，由高校、科研院所、臨床醫生構成，由於群人數超過100，需要由工作人員拉你才能入群，工作人員二維碼如下，添加時請備註「申請入群」：

科研統計交流群

由高校、醫院、科研院所等科研人員組成，旨在搭建交流平台，添加，請備註「申請入群」

喜歡這篇文章嗎？立刻分享出去讓更多人知道吧！