10=chap9非參數(shù)檢驗ridit分析

1、Chap9非參數(shù)檢驗、Ridit分析 P139155 教學(xué)目的與要求： 3掌握：非參數(shù)統(tǒng)計概念、優(yōu)缺點，配對、兩組、多組資料的秩和檢驗及兩兩比較。熟悉：Ridit分析的意義和分析方法。了解：游程檢驗。 教學(xué)內(nèi)容提要 重點講解：配對、兩組、多組秩和檢驗及兩兩比較。講解：Ridit分析的意義和分析方法。介紹：游程檢驗。 重點：配對、兩組、多組資料的秩和檢驗及兩兩比較。難點：Ridit分析，游程檢驗。教學(xué)時間: 理論3學(xué)時。13 , 45, 67 2022/7/2711非參數(shù)檢驗的意義 P1391. 比較 1.1 參數(shù)檢驗（parametric test）：t、Z、F，相關(guān)、回歸等。 總體分布的類型

2、已知，對未知參數(shù)進(jìn)行檢驗。優(yōu)點：信息利用充分；檢驗效能較高。缺點：資料要求較高；適用范圍有限。1.2非參數(shù)檢驗（nonparametric test）任意分布檢驗（distribution-free test）：秩和（rank sum test）、 2。與總體分布無關(guān)，不比較參數(shù)，比較分布的位置（中位數(shù)M） 。優(yōu)點：適用范圍 。缺點：信息利用 ,效能,型錯誤,假(-)率。2022/7/2722Wilcoxon符號秩和檢驗（配對及單樣本比較用） P139一、配對設(shè)計資料 【例9-1】對12份血清分別用原法（檢測時間20）和新法（10）檢測谷-丙轉(zhuǎn)氨酶，結(jié)果見表9-1，兩種檢測結(jié)果是否相同？表9-

3、1 原法和新法檢測12份血清谷-丙轉(zhuǎn)氨酶（nmolS1/L）資料對子編號123456789101112原 法6023695802422121902522038142195新 法80200100822402432053822044152243差值d-2036-5-22-31-15-130-6-10-48編 秩-810-3-1.51.5-9-7-6-4-5-11對d進(jìn)行正態(tài)性檢驗，P0.05 配對t test P0.05 Wilcoxon符號秩和檢驗 (d =0)2022/7/273以原法、新法為變量名配對計量的數(shù)據(jù)文件L9-1.sav 配對符號秩和檢驗。Analyze Nonparametric

4、 Tests（非參數(shù)檢驗）2 Related Samples 新法 、原法同時送入Test Pairs List框中OK。P0.05，不拒絕H0，不能認(rèn)為兩法檢測結(jié)果不同。 2022/7/2742022/7/275二、單樣本比較的秩和檢驗（樣總比較） 【例9-2】已知某地正常人尿氟的中位數(shù)為0.86mg/L，今在該地某廠隨機(jī)抽取12名工人，測得尿氟含量（mg/L）為：0.84，0.86，0.88，0.94，0.97，1.01，1.05，1.09，1.20，1.28，1.35，1.83。推斷該廠工人的尿氟含量是否高于當(dāng)?shù)卣Ｈ恕?與配對不同：與每一個樣本觀測值配對的是中位數(shù)，建立數(shù)據(jù)文件L9-2

5、.sav。操作過程同例9-1。 Analyze Nonparametric Tests2 Related Samples 尿氟含量、中位數(shù)同時送入Test Pairs ListOK。2022/7/2762022/7/2773 Wilcoxon秩和& Mann-Whitney U檢驗(兩組) 141一、兩連續(xù)型資料比較的秩和檢驗 原始數(shù)據(jù) 兩小n Wilcoxon秩和檢驗 兩大n 近似正態(tài)Z檢驗 頻數(shù)表資料(近似等級) 兩小n Mann-Whitney U 兩大n 近似正態(tài)Z檢驗2022/7/2781.Wilcoxon秩和 ：2小n原始數(shù)據(jù)：n110，(n2n1) 0.10。2.Mann-Whi

6、tney U檢驗（Mann-Whitney U test）：兩小n頻數(shù)表n1n2，將n1中的每個變量值與n2個變量值逐個比較，小記1，等記0.5，大記0，求其和得U。2022/7/2710【例9-3】以葛根劑量和抑制率為變量名，在葛根劑量的Values中標(biāo)簽1=3g/100ml組，2=5g/100ml組，建立13行2列的完全隨機(jī)設(shè)計計量資料數(shù)據(jù)文件L9-3.sav以后，進(jìn)行兩樣本比較的Wilcoxon秩和檢驗和Mann-Whitney U檢驗：Analyze Nonparametric Tests2 Independent Samples，抑制率 入Test variable list，葛根劑

7、量入grouping variables Define groups，group1和group2中輸1和2 Continue OK 2022/7/27112022/7/2712二、等級資料兩樣本比較的Mann-Whitney U檢驗【例9-4】用祖?zhèn)麽樉姆ㄅc一般針灸療法各治療哮喘病人46例、28例，以痰液嗜酸性白細(xì)胞為指標(biāo)，結(jié)果見表9-3，判斷祖?zhèn)麽樉姆ǖ寞熜欠窀摺?2022/7/2713以組別、療效和例數(shù)為變量名，建立8行3列的數(shù)據(jù)文件L9-4.savAnalyze Nonparametric Tests2 Independent Samples，療效 入Test variable lis

8、t，組別入grouping variables Define groups，在group1和group2中輸入1和2 Continue OK。 2022/7/27142022/7/27154 Kruskal-wallis H 檢驗（多組比較用） 一、多組連續(xù)型資料1. Kruskal-wallis H檢驗： 計量：完全隨機(jī)設(shè)計兩或多樣本比較 等級：推斷雙向有序資料的相關(guān)性 【例9-5】比較A、B、C中藥的滅釘螺作用，各觀察5批，每批對200只活釘螺用藥后清點死亡數(shù)，計算死亡率（），結(jié)果見表9-4的、行。2022/7/2716以藥物和死亡率為變量名，藥物Values標(biāo)簽1=A藥，2=B藥，3C藥

9、，建立2列15行的數(shù)據(jù)文件L9-5.savAnalyze Nonparametric Test K Independent Samples，死亡率入Test variable list，藥物入grouping variables Define Range ，在Minimun和Maximun中填入1和3（第1至第3組間的比較） Continue。在Test Type復(fù)選框中有Kruskal-Wallis H、Median（中位數(shù)）、Jonckheere-Terpstra，系統(tǒng)默認(rèn)或全選OK2022/7/27172022/7/27182.兩兩比較的Nemenyi test2022/7/2719DP

10、S操作 （1）正態(tài)性檢驗：分別選定各組數(shù)據(jù)塊數(shù)據(jù)分析正態(tài)性檢驗。 若 P0.05，可認(rèn)為來自正態(tài)總體。 可用ANOVA。（2） Kruskal-wallis H檢驗：選定數(shù)據(jù)塊分類數(shù)據(jù)統(tǒng)計非參數(shù)檢驗 Kruskal-wallis H檢驗。 2022/7/27202022/7/27212022/7/27222022/7/2723二、等級資料多樣本比較的秩和檢驗【例9-6】某中醫(yī)院以蛞蝓膠囊為主綜合治療（簡稱蛞蝓綜合法）中晚期肺癌，并與中西醫(yī)結(jié)合治療及聯(lián)合化療作比較觀察，其近期療效分為部分緩解，穩(wěn)定，擴(kuò)展三級，資料見表9-5，比較三組的療效。 2022/7/2724以療法、療效、例數(shù)為變量名，對療

11、法、療效標(biāo)簽其水平，建立3列9行的RC表文件L9-6.sav。Dataweight casesweight cases by例數(shù)AnalyzeNonparametricTestKIndependent Samples，療效入Test variable list，療法入grouping variables Define Range ，在Minimun和Maximun填入1和3（第1至第3組間比較） ContinueOK。結(jié)果： 2 =4.030，P=0.133。 2022/7/27252022/7/27265 Friedman M檢驗、Cochrans Q檢驗（配伍資料比較用） 一、Friedm

12、an M 秩和檢驗：配伍組設(shè)計連續(xù)資料【例9-7】某醫(yī)師按中醫(yī)辨證，把肺癌病人分成五類，研究辨證分型的療效，由于療效又受病期的影響，所以又按病期分為、三個配伍組（b3），資料見表9-6。分析不同辨證分型肺癌病人的一年生存率是否不同。 2022/7/2727以病期、陽虛、氣陰兩虛、氣虛、陰陽兩虛、氣滯血瘀為變量名，建立3行6列配伍組設(shè)計連續(xù)型資料數(shù)據(jù)文件L9-7.sav。Analyze Nonparametric Test K Related Samples，陽虛、氣陰兩虛、氣虛、陰陽兩虛、氣滯血瘀進(jìn)Test VariablesOK。 結(jié)果： 2 =5.867， P=0.2092022/7/27

13、282022/7/2729二、Cochrans Q檢驗 ：配伍組設(shè)計兩分類資料。 【例9-8】 在嬰兒哭鬧時，可以采用不同的哄法，如A（喂水）、B（輕輕搖晃）、C（哄以橡皮奶頭）、D（說話逗他）。對12名未滿一月的嬰兒用這4種方法進(jìn)行試驗，以1表示哄法有效、0表示無效，結(jié)果如表9-8。判斷這4種哄嬰兒方法的效果是否不同。2022/7/2730以編號、哄法A、B、C、D為名，建立5列12行L9-8.sav. Analyze Nonparametric Test K Related Samples，A、B、C、D入Test Variable List，Test Type只選 Cochrans Q

14、OKCochrans Q =3.692，P=0.297 2022/7/27312022/7/27326 Ridit分析 一、Ridit分析的意義、R值均數(shù)的估計 Ridit：與特定分布相對應(yīng)（relative to an identified distribution，Rid）的單位（unit）。適用于有序分類變量比較（各組n50），說明各等級是否相同及優(yōu)劣。等級數(shù)據(jù)定量（連續(xù)型） 按正態(tài)分布的理論推斷基本思想：選標(biāo)準(zhǔn)組（大n, 基準(zhǔn)） 用標(biāo)準(zhǔn)組資料計算各等級的Ridit值（即參照單位）用參照單位計算各對比組的平均Ridit值推斷。2022/7/2733二、樣本與總體比較的Ridit分析【例9

15、-9】 某醫(yī)院用中草藥制劑甲、乙治療慢支病人結(jié)果如表9-9，乙與甲制劑相比，療效如何？表9-9 兩制劑對慢性支氣管炎的療效制劑無效好轉(zhuǎn)顯效痊愈合計甲7761926685333420乙115524151052022/7/27341.建立數(shù)據(jù)文件 ：療效、甲組、乙組為名 L9-9.sav。2.計算標(biāo)準(zhǔn)組各等級的Ridit值（1）將標(biāo)準(zhǔn)組的例數(shù)加權(quán)：甲。DataWeight Cases Weight Cases by甲組例數(shù) Frequency VariableOK（2）計算標(biāo)準(zhǔn)組療效各等級的Ridit值：TransformRank Case（觀察單位排秩） 療效入variableRank Type

16、s （排秩方式）， Proportion estimate（百分位數(shù)估計） Continue OK2022/7/27353.計算對比組總體R均值及其置信區(qū)間（1）將對比組的例數(shù)變量進(jìn)行加權(quán)：乙DataWeight Cases，F(xiàn)requency Variable甲出乙進(jìn)OK（2）求總體R均值及其95置信區(qū)間：analyzedescriptive Statistics Explore，將PROPORTION of 療效 using BLOM p療效送入Dependent List框OK2022/7/27362022/7/27372022/7/2738三、兩樣本比較的Ridit分析【例9-10】觀

17、察針刺與藥物治療內(nèi)分泌性突眼癥的療效，資料見表9-11，比較兩療法的療效有無差別。表9-11 針刺與藥物治療的效果比較組 別痊愈顯效好轉(zhuǎn)無效合計針刺組193281372藥療組3642942合 計223812421142022/7/2739（1）建數(shù)據(jù)文件：療效、針刺、藥療為名 L9-10.sav（2）確定標(biāo)準(zhǔn)組：Transform Compute，Target Variable輸標(biāo)準(zhǔn)例數(shù)針刺 “+” 藥療 OK（3）計算標(biāo)準(zhǔn)組療效各等級R：A. 加權(quán):DataWeight Cases 標(biāo)準(zhǔn)例數(shù)進(jìn)OKB.求標(biāo)準(zhǔn)組R:TransformRank Case,療效variableRank Types

18、Proportion estimate Continue OK。p療效（4）求對比組：DataWeight Cases，標(biāo)出針進(jìn)OK analyzedescriptive Statistics Explore，PROPORTION of 療效 入Dependent ListOK DataWeight Cases，針刺出藥療進(jìn)OK analyzedescriptive Statistics ExploreOK2022/7/27402022/7/27412022/7/2742針刺=0.3969，總體R針刺值95的CI為（0.3391，0.4548）；藥療0.6767，總體R藥療值95的CI為（0.6043，0.7491）。兩組95CI無重疊，拒絕H02022/7/2743四、多樣本比較的Ridit分析 【例9-11】某醫(yī)院在進(jìn)行黃連、黃柏雙黃氣霧劑治療急性呼吸道感染的課題研究中，對三種不同劑型的雙黃制劑的療效進(jìn)行了比較，其結(jié)果見