非參數(shù)檢驗

基礎所流行病王艷紅流行病學資料分析2022/10/181非參數(shù)檢驗

基礎所流行病流行病學資料分析2022/10/15一、相關理論知識的回顧非參數(shù)檢驗特點(定期復習十分重要)總體分布不作嚴格規(guī)定；不依賴于總體分布類型（如正態(tài)分布）

——任意分布檢驗非參數(shù)檢驗的適用條件總體的分布不易確定偏態(tài)分布方差不齊等級資料（統(tǒng)計資料一般分為計量資料與計數(shù)資料兩大類，介于其間的一類資料稱為等級資料.不同類型的資料應采用不同的分析方法。例如用某藥治1.5指標的特異性指測量特定項目而不療若干痢疾患者,其中治愈、好轉、無效人受其它項目干擾的能力文獻來源；化驗結果分為:-、+、++、+++、++++5個有序等級）2一、相關理論知識的回顧非參數(shù)檢驗特點(定期復習十分重要)2一、相關理論知識的回顧基于秩次的

非參數(shù)檢驗

3一、相關理論知識的回顧基于秩次的

非參數(shù)檢驗

3一、相關理論知識的回顧（1）、Wilcoxon秩和檢驗應用兩組連續(xù)變量資料的秩和檢驗兩組有序變量資料的秩和檢驗基本思想：假設兩總體的分布相同（H0），兩樣本可認為是從同一總體中抽取的隨機樣本；將二者混合后由小到大編秩，然后分別計算兩組的平均秩和R1與R2，R1和R2如果近似相等，其差別可能是由隨機抽樣引起的。4一、相關理論知識的回顧（1）、Wilcoxon秩和檢驗4一、相關理論知識的回顧5一、相關理論知識的回顧5一、相關理論知識的回顧（1）、Wilcoxon秩和檢驗公式1——沒有節(jié)點（即沒有相同的觀察值）公式2——有節(jié)點6一、相關理論知識的回顧（1）、Wilcoxon秩和檢驗6一、相關理論知識的回顧（2）、Kruskal-Wallis秩和檢驗應用多組連續(xù)變量資料的秩和檢驗多組有序變量資料的秩和檢驗基本思想同Wilcoxon秩和檢驗7一、相關理論知識的回顧（2）、Kruskal-Wallis秩一、相關理論知識的回顧（2）、Kruskal-Wallis秩和檢驗，當K=3且至少有一組例數(shù)大于5時，或者k>3時服從卡方分布；若K=3且沒有一組例數(shù)大于5時，服從H分布。公式1——沒有結公式2——有結8一、相關理論知識的回顧（2）、Kruskal-Wallis秩一、相關理論知識的回顧（3）、Friedman秩和檢驗應用隨機區(qū)組設計資料基本思想各區(qū)組的觀察值按由小到大的順序進行編秩；如果各處理的效應相同，各區(qū)組內秩1，2，……，k應以相同的概率出現(xiàn)在各處理組（列）中，各處理組的秩和應該大致相等，不可能出現(xiàn)較大的差別。9一、相關理論知識的回顧（3）、Friedman秩和檢驗9（3）、Friedman秩和檢驗公式10（3）、Friedman秩和檢驗10二、相關SAS過程步（1）、NPAR1WAY過程兩組連續(xù)變量資料的秩和檢驗多組連續(xù)變量資料的秩和檢驗（2）、RANK過程+ANOVA(或GLM)過程完全隨機設計連續(xù)變量資料的秩和檢驗隨機區(qū)組設計資料的秩和檢驗多重比較（3）、Freq過程等級資料兩組連續(xù)變量資料的秩和檢驗完全隨機設計連續(xù)變量資料的秩和檢驗隨機區(qū)組設計資料的秩和檢驗單因素單因素、雙因素、兩兩比較單因素、雙因素、等級資料11二、相關SAS過程步（1）、NPAR1WAY過程單因素單因素二、相關SAS過程步（1）、NPAR1WAY過程的基本結構procnpar1waywilcoxon;classgroup（分組變量）;varx（連續(xù)變量）;run;Wilcoxon是指資料的分析基于Wilcoxon得分兩樣本資料：Wilcoxon秩和檢驗多樣本作Kruskal-Wallis秩和檢驗12二、相關SAS過程步（1）、NPAR1WAY過程的基本結構1二、相關SAS過程步（2）、RANK過程+ANOVA(或GLM)過程基本結構完全隨機設計（Kruskal-Wallis秩和檢驗）procrankdata=sampleout=a;varx;ranksr;

/*產生秩次*/run;procglmdata=a;classgroup;modelr=group;meansgroup/snk;run;procsortdata=a;byblock;run;procrankdata=aout=b;varx;byblock;ranksr;run;procglmdata=b;classtreatblock;modelr=treatblock;meanstreat/snk;run;隨機區(qū)組設計（Friedman秩和檢驗）13二、相關SAS過程步（2）、RANK過程+ANOVA(或GL二、相關SAS過程步（3）、Freq過程的基本結構（CMH解法）procfreqdata=aa;tablesgroup*x/scores=rank

cmh;run;14二、相關SAS過程步（3）、Freq過程的基本結構（CMH解三、常見的非參數(shù)檢驗實例（一）、配對設計和單樣本資料的非參數(shù)檢驗（二）、完全隨機設計兩組獨立樣本的秩和檢驗（三）、完全隨機設計多組獨立樣本的秩和檢驗（四）、隨機區(qū)組設計資料的秩和檢驗15三、常見的非參數(shù)檢驗實例（一）、配對設計和單樣本資料的非參數(shù)（一）、配對資料的非參數(shù)檢驗擬對某種藥物的降壓效果進行評價，治療前后病人的舒張壓水平如下：試用非參數(shù)檢驗的方法比較治療前后舒張壓有無差異？datadbp;inputbeforeafter@@;d=before-after;cards;11712312710814112010710711010011498115102138152127104122107;run;procunivariatedata=dbpnormal;vard;run;16（一）、配對資料的非參數(shù)檢驗datadbp;procun（二）、完全隨機設計兩組獨立樣本的秩和檢驗實例某缺碘地區(qū)測定9名12歲男孩和8名女孩的尿碘含量（μg/l），結果如下：男孩（9名）：102.00、45.70、103.00、105.00、98.30、113.00、64.00、113.00、54.00女孩（8名）：61.00、49.70、45.80、33.90、116.00、118.00、120.00、72.00試用非參數(shù)檢驗的方法比較不同性別的兒童尿碘含量是否不同?17（二）、完全隨機設計兩組獨立樣本的秩和檢驗17SAS數(shù)據(jù)集的建立dataaa;if_n_<10thensex=1;if_n_>=10thensex=2;inputx@@;cards;102.0045.70103.00105.0098.30113.0064.00113.0054.0061.0049.7045.8033.90116.00118.00120.0072.00;run;SAS過程步procnpar1waywilcoxondata=aa;varx;classsex;run;procfreqdata=aa;tablessex*x/scores=rankcmh;run;行均值得分的差值Rowmeanscoresdiffer18SAS數(shù)據(jù)集的建立dataaa;SAS過程步procnp練習1應用免疫酶法對鼻咽癌患者和非鼻咽癌者分別測定13人的血清病毒VCA-LOGA抗體滴度（倍），結果如下鼻咽癌患者：5204080808016016032032064012801280非鼻咽癌患者：510102020204040808080160640試用非參數(shù)檢驗法比較兩組患者的血清病毒的VCA-LOGA抗體滴度是否有差異？19練習119（三）、完全隨機設計多組獨立樣本的秩和檢驗實例某研究者欲研究A、B兩個菌種對小鼠巨噬細胞吞噬功能的激活作用，將60只小鼠隨機分為三組，其中生理鹽水為對照組，用常規(guī)巨噬細胞吞噬功能的監(jiān)測方法，獲得三組吞噬指數(shù)如下：對照組：1.301.401.501.501.601.601.701.701.701.701.702.102.302.302.302.402.40A菌組：1.801.802.202.202.202.202.302.302.302.662.662.682.682.702.802.803.003.103.103.103.103.204.304.30B菌組：1.501.801.802.002.302.302.402.402.402.402.502.502.602.602.602.702.703.103.10試用非參數(shù)檢驗的方法比較三組吞噬指數(shù)有無差別？20（三）、完全隨機設計多組獨立樣本的秩和檢驗20SAS數(shù)據(jù)集的建立dataaa;if_n_<18thengroup=1;if_n_>=18and_n_<42thengroup=2;if_n_>=42thengroup=3;inputx@@;cards;1.301.401.501.501.601.601.701.701.701.701.702.102.302.302.302.402.401.801.802.202.202.202.202.302.302.302.662.662.682.682.702.802.803.003.103.103.103.103.204.304.301.501.801.802.002.302.302.402.402.402.402.502.502.602.602.602.702.703.103.10;run;SAS過程步procnpar1waywilcoxondata=aa;varx;classgroup;run;procrankdata=aaout=bb;varx;ranksr;run;procglmdata=bb;classgroup;modelr=group;meansgroup/dunnett("1");run;procfreqdata=aa;tablesgroup*x/scores=rankcmh;run;Kruskal-Wallis檢驗結果：R-Square×Df總行均值得分的差值Rowmeanscoresdiffer21SAS數(shù)據(jù)集的建立dataaa;SAS過程步procnp練習2正常人、單純型肥胖和皮質醇增多癥病人的血漿總皮質醇測定值如下（μg/L）：正常人：1.41.92.22.52.83.13.73.94.67.0單純型肥胖：1.61.22.02.43.14.15.07.410.813.6皮質醇增多癥：9.810.210.613.014.014.815.615.621.624.0試用非參數(shù)的方法比較三組人群的血漿總皮質醇含量有無差別？22練習222（四）、隨機區(qū)組設計資料的秩和檢驗實例用某新藥治療血吸蟲病患者，采用三天療法，在治療前及治療后測定7名患者的血清谷丙轉氨酶（SGPT）的變化，以觀察該藥對肝功能的影響，測定結果如下。試用非參數(shù)檢驗的方法比較四個階段SGPT有無差別？23（四）、隨機區(qū)組設計資料的秩和檢驗23SAS數(shù)據(jù)集的建立dataaa;doblock=1to7;dotreat=1to4;inputx@@;output;end;end;cards;63188138549023822014454300839245140213100541751503672300163906420718587;run;SAS過程步procsortdata=aa;byblock;run;procrankdata=aaout=b;varx;byblock;ranksr;run;procglmdata=b;classtreatblock;modelr=treatblock;run;procfreqdata=aa;tablesblock*treat*x/scores=rankcmh;run;行均值得分的差值RowmeanscoresdifferFriedman檢驗結果：R-Square×（Df總－Df區(qū)組）24SAS數(shù)據(jù)集的建立dataaa;SAS過程步procso小結（1）配對設計和單樣本資料的非參數(shù)檢驗Wilcoxon符號檢驗和Wilcoxon符號秩和檢驗SAS中Univariate過程：“位置檢驗:Mu0=0”中的符