衛(wèi)生統(tǒng)計(jì)學(xué)（教學(xué)講義）第十一章無(wú)序分類資料的統(tǒng)計(jì)分析的sas實(shí)現(xiàn)

1、第十一章 無(wú)序分類資料的統(tǒng)計(jì)分析的SAS實(shí)現(xiàn)例11-2 某研究小組為研究慢支口服液II號(hào)對(duì)慢性支氣管炎治療效果，以口服消咳喘為對(duì)照進(jìn)行了臨床試驗(yàn)，試驗(yàn)組120人、對(duì)照組117人（兩組受試者病程、病情等均衡），療程2周，療效見(jiàn)表11-3。問(wèn)慢支口服液II號(hào)與消咳喘治療慢性支氣管炎的療效是否相同？表11-3 試驗(yàn)組與對(duì)照組療效組別有效無(wú)效合計(jì)有效率（%）試驗(yàn)組116412096.67對(duì)照組823511770.09合計(jì)1983923783.54SAS分析程序data ch11_2; 定義數(shù)據(jù)集名input r c f; 輸入行變量，列變量以及頻數(shù)變量cards;1 1 116 1 2 4 2 1 8

2、2 2 2 35;proc freq; 利用freq過(guò)程進(jìn)行Pearson檢驗(yàn)weight f; 指定頻數(shù)變量tables r*c/chisq expected; 指定行變量，列變量并輸出理論頻數(shù)run;SAS軟件輸出結(jié)果The FREQ Procedure Statistics for Table of r by c Statistic DF Value Prob Chi-Square 1 30.4463 .0001 Likelihood Ratio Chi-Square 1 34.1018 .0001 Continuity Adj. Chi-Square 1 28.5435 .0001 M

3、antel-Haenszel Chi-Square 1 30.3178 .0001 Phi Coefficient 0.3584 Contingency Coefficient 0.3374 Cramers V 0.3584 Fishers Exact Test Cell (1,1) Frequency (F) 116 Left-sided Pr = F 9.747E-09 Table Probability (P) 9.006E-09 Two-sided Pr = P 1.178E-08 Sample Size = 237SAS軟件輸出結(jié)果解釋該結(jié)果包含兩個(gè)部分，第一部分為檢驗(yàn)結(jié)果，列出了各

4、種檢驗(yàn)方法，本例考慮Pearson （Chi-Square）檢驗(yàn)方法，統(tǒng)計(jì)量值為30.4463，p值.0001，因此按=0.05水準(zhǔn)拒絕，差別有統(tǒng)計(jì)學(xué)意義。另外，教材中提到校正檢驗(yàn)，輸出結(jié)果中的Continuity Adj. Chi-Square就是校正檢驗(yàn)；第二部分給出的是Fisher精確概率檢驗(yàn)方法（Fishers Exact Test）。教材中的說(shuō)明查附表8， ， ，按=0.05水準(zhǔn)拒絕，差別有統(tǒng)計(jì)學(xué)意義，可認(rèn)為慢支口服液II號(hào)治療慢性支氣管炎有效率高于消咳喘。教材中關(guān)于校正檢驗(yàn)，F(xiàn)isher精確檢驗(yàn)方法的使用有如下說(shuō)明：分布是服從標(biāo)準(zhǔn)正態(tài)分布的隨機(jī)變量的平方和所形成的統(tǒng)計(jì)量的概率分布，

5、是一連續(xù)型分布，。而四格表資料是二分類變量資料，因而統(tǒng)計(jì)量的抽樣分布是非連續(xù)的。在大樣本情況下（一般指且各個(gè)格子的理論頻數(shù)），二項(xiàng)分布逼近正態(tài)分布，統(tǒng)計(jì)量近似地服從于自由度為1的分布。當(dāng)不滿足以上條件時(shí)，非校正的值偏大，因此不宜采用基于連續(xù)分布的檢驗(yàn)，而應(yīng)采用校正的檢驗(yàn)或改用精確概率法。一般認(rèn)為，若，采用校正檢驗(yàn)或Fisher精確概率檢驗(yàn) ；若或，采用Fisher精確概率檢驗(yàn)。例11-3 為評(píng)價(jià)中西結(jié)合治療抑郁發(fā)作的療效。將187例患者隨機(jī)分為2組，兩組患者均選用阿咪替林西醫(yī)綜合治療，中西醫(yī)結(jié)合組在上述治療的同時(shí)，再配合中醫(yī)辨證治療，根據(jù)中醫(yī)辨證分型采用不同的方劑，治療結(jié)果見(jiàn)表11-5，問(wèn)兩種

6、治療方案的療效有無(wú)差別？ 表11-5 試驗(yàn)組與對(duì)照組療效組別有效無(wú)效合計(jì)有效率（%）中西醫(yī)結(jié)合組92（88.973）2（5.027）9497.87西醫(yī)組85（88.027）8（4.973）9391.40注 ：括號(hào)內(nèi)為理論頻數(shù)SAS分析程序data ch11_3; 定義數(shù)據(jù)集名input r c f; 輸入行變量，列變量以及頻數(shù)變量cards;1 1 92 1 2 2 2 1 85 2 2 8;proc freq; 利用freq過(guò)程進(jìn)行Pearson檢驗(yàn)weight f; 指定頻數(shù)變量tables r*c/chisq expected; 指定行變量，列變量并輸出理論頻數(shù)run;SAS軟件輸出結(jié)果

7、The FREQ Procedure Statistics for Table of r by c Statistic DF Value Prob Chi-Square 1 3.8716 0.0491 Likelihood Ratio Chi-Square 1 4.1265 0.0422 Continuity Adj. Chi-Square 1 2.6981 0.1005 Mantel-Haenszel Chi-Square 1 3.8509 0.0497 Phi Coefficient 0.1439 Contingency Coefficient 0.1424 Cramers V 0.143

8、9 WARNING: 25% of the cells have expected counts less than 5. Chi-Square may not be a valid test. Fishers Exact Test Cell (1,1) Frequency (F) 92 Left-sided Pr = F 0.0482 Table Probability (P) 0.0395 Two-sided Pr = P 0.0578 Sample Size = 187SAS軟件輸出結(jié)果解釋本例考慮校正Chi-Square檢驗(yàn)方法，統(tǒng)計(jì)量值為2.6981，p值.0.1005，因此按=0.

9、05水準(zhǔn)不能拒絕，差別沒(méi)有統(tǒng)計(jì)學(xué)意義教材中的說(shuō)明由于有一理論數(shù)為4.973，需采用校正檢驗(yàn)或Fisher精確概率檢驗(yàn)。查附表8， ，按=0.05水準(zhǔn)不拒絕，差別無(wú)統(tǒng)計(jì)學(xué)意義，還不能認(rèn)為中西醫(yī)結(jié)合治療的有效率不同于西醫(yī)組。例11-4 為了解國(guó)產(chǎn)紫外線瞬間消毒器與進(jìn)口高壓蒸汽消毒機(jī)對(duì)牙科手機(jī)消毒滅菌的效果，將剛?cè)ジ铱ㄟ^(guò)垢等待處理的牙科手機(jī)29個(gè)隨機(jī)分為A、B兩組，A組為紫外線消毒組，B組為高壓蒸汽組。消毒前細(xì)菌培養(yǎng)均為陽(yáng)性，消毒后細(xì)菌培養(yǎng)結(jié)果見(jiàn)下表。問(wèn)兩種消毒法消毒后細(xì)菌培養(yǎng)陽(yáng)性率有無(wú)差別？表11-6 兩種方法消毒后細(xì)菌培養(yǎng)結(jié)果組別陽(yáng)性陰性合計(jì)A10515B11314合計(jì)111829SAS分析程

10、序data ch11_4; 定義數(shù)據(jù)集名input r c f; 輸入行變量，列變量以及頻數(shù)變量cards;1 1 10 1 2 5 2 1 1 2 2 13;proc freq; 利用freq過(guò)程進(jìn)行Pearson檢驗(yàn)weight f; 指定頻數(shù)變量tables r*c/chisq expected; run;SAS軟件輸出結(jié)果The FREQ Procedure Statistics for Table of r by c Statistic DF Value Prob Chi-Square 1 10.8977 0.0010 Likelihood Ratio Chi-Square 1 12

11、.1957 0.0005 Continuity Adj. Chi-Square 1 8.5160 0.0035 Mantel-Haenszel Chi-Square 1 10.5219 0.0012 Phi Coefficient 0.6130 Contingency Coefficient 0.5226 Cramers V 0.6130 Fishers Exact Test Cell (1,1) Frequency (F) 10 Left-sided Pr = F 0.0013 Table Probability (P) 0.0012 Two-sided Pr = P 0.0017 Samp

12、le Size = 29SAS軟件輸出結(jié)果解釋本例考慮Fisher精確概率檢驗(yàn)，雙邊檢驗(yàn)p值=0.0017，因此按=0.05水準(zhǔn)拒絕，差別有統(tǒng)計(jì)學(xué)意義。教材中的說(shuō)明本例為雙側(cè)檢驗(yàn)，小于或等于的概率有、，F(xiàn)isher精確概率檢驗(yàn)的，差別有統(tǒng)計(jì)學(xué)意義，按=0.05水準(zhǔn)拒絕，可認(rèn)為紫外線消毒法細(xì)菌培養(yǎng)陽(yáng)性率高于高壓蒸汽法。例11-5 為探討埃茲蛋白（Ezrin）在胃癌組織中的表達(dá)情況，采用免疫組化法檢測(cè)50例胃癌組織、25例胃粘膜不典型增生和25例正常胃粘膜中Ezrin的表達(dá)，結(jié)果見(jiàn)表11-9。問(wèn)不同胃組織Ezrin表達(dá)陽(yáng)性率是否相同？表11-9 Ezrin在不同胃組織中的表達(dá)組別觀測(cè)例數(shù)陽(yáng)性例數(shù)陰

13、性例數(shù)陽(yáng)性率（%）正常胃粘膜257（15.250）18（9.750）28.0不典型增生2511（15.250）14（9.750）44.0胃癌組織5043（30.500）7（19.500）86.0合計(jì)100613961.0注 ：括號(hào)內(nèi)為理論頻數(shù)SAS分析程序data ch11_5; 定義數(shù)據(jù)集名do r=1 to 3; 定義行變量do c=1 to 2; 定義列變量input f;output;end;end;cards;7 18 11 14 43 7;proc freq; 利用freq過(guò)程進(jìn)行Pearson檢驗(yàn)weight f; 指定頻數(shù)變量tables r*c /chisq; 指定行變量，列

14、變量run;SAS軟件輸出結(jié)果Statistics for Table of r by cStatistic DF Value ProbChi-Square 2 27.6166 .0001Likelihood Ratio Chi-Square 2 29.3091 .0001Mantel-Haenszel Chi-Square 1 26.0618 .0001Phi Coefficient 0.5255Contingency Coefficient 0.4652Cramers V 0.5255Sample Size = 100SAS軟件輸出結(jié)果解釋本例考慮Pearson檢驗(yàn)，p值 S 0.0719

15、 Simple Kappa Coefficient Kappa 0.4589 ASE 0.0907 95% Lower Conf Limit 0.2811 95% Upper Conf Limit 0.6366 Sample Size = 100PART IIStatistics for Table of r by c Statistic DF Value Prob Chi-Square 1 21.8872 .0001 Likelihood Ratio Chi-Square 1 22.1508 .0001 Continuity Adj. Chi-Square 1 19.8710 .0001 M

16、antel-Haenszel Chi-Square 1 21.6683 .0001 Phi Coefficient 0.4678 Contingency Coefficient 0.4238 Cramers V 0.4678 Fishers Exact Test Cell (1,1) Frequency (F) 52 Left-sided Pr = F 3.821E-06 Table Probability (P) 3.423E-06 Two-sided Pr = P 5.240E-06 Sample Size = 100SAS軟件輸出結(jié)果解釋PART I是配對(duì)四格表資料兩種處理的差別分析，因此考慮McNemars 檢驗(yàn)（McNemars Test），p值=0.0719，因此按=0.05水準(zhǔn)不能拒絕，差別沒(méi)有統(tǒng)計(jì)學(xué)意義