醫(yī)學(xué)統(tǒng)計(jì)學(xué)：非參數(shù)檢驗(yàn)SPSS

非參數(shù)檢驗(yàn)的SPSS實(shí)現(xiàn)一、兩獨(dú)立樣本的檢驗(yàn)

[例1]調(diào)查某廠的鉛作業(yè)工人7人和非鉛作業(yè)工人10人的血鉛值（μg/100g）如下，問(wèn)兩組工人的血鉛值有無(wú)差別？

非鉛作業(yè)組556791213151821鉛作業(yè)組17182025344344數(shù)據(jù)準(zhǔn)備

定義2個(gè)變量：

1.group（非鉛作業(yè)組為1，鉛作業(yè)組為2），2.Pb。按順序輸入數(shù)據(jù)。

統(tǒng)計(jì)分析

分析→非參數(shù)統(tǒng)計(jì)→兩獨(dú)立樣本→Ｐb選入檢驗(yàn)變量列表框，group選入分組變量框→定義組別，定義組1為1，組2為2→檢驗(yàn)類型：Mann-WhitneyU→確定結(jié)果解釋

結(jié)果表明，第1組的平均秩次（MeanRank）為5.95，第2組的平均秩次為13.36，Mann-Whitney

U=4.5，W=59.5，z=-2.98,精確雙側(cè)概率P=0.001，可認(rèn)為鉛作業(yè)組工人的血鉛值高于非鉛作業(yè)組。

二、多個(gè)獨(dú)立樣本的檢驗(yàn)

[例2]隨機(jī)抽樣得右表三組人的血槳總皮質(zhì)醇測(cè)定值（μg/L），試比較三組人有無(wú)差異？正常人單純性肥胖皮質(zhì)醇增多癥0.40.69.81.91.210.22.22.010.62.52.413.02.83.114.03.14.114.83.75.015.63.95.915.64.67.421.67.013.624.0數(shù)據(jù)準(zhǔn)備

定義2個(gè)變量：1.group（正常人為1，單純性肥胖為2，皮質(zhì)醇增多癥為3），2.x(總皮質(zhì)醇測(cè)定值)

按順序輸入數(shù)據(jù)。統(tǒng)計(jì)分析

分析→非參數(shù)統(tǒng)計(jì)→Ｋ個(gè)獨(dú)立樣本→x選入檢驗(yàn)變量框，組別選入分組變量框→定義組別，最小值為1，最大值為3→確定結(jié)果解釋

結(jié)果表明，1至3組的平均秩次（MeanRank）分別為9.65、11.75、25.10，χ2值（即H值）為18.1219，P=0.0001；可認(rèn)為三組人的血槳總皮質(zhì)醇測(cè)定值有差異，根據(jù)本例情況可看出皮質(zhì)醇增多癥組高于其他兩組人。大鼠對(duì)別正常飼料組VitE缺乏飼料組137.225.7220.925.1331.418.8441.433.5539.834.0639.328.3736.126.2831.918.3表2正常飼料組、VitE缺乏飼料組大鼠肝中VitA含量三、配對(duì)資料秩和檢驗(yàn)[例3]數(shù)據(jù)準(zhǔn)備

兩個(gè)變量:

1.正常飼料2.VitE

缺乏

按順序輸入數(shù)據(jù)。

統(tǒng)計(jì)分析

分析→非參數(shù)統(tǒng)計(jì)→兩個(gè)相關(guān)樣本檢驗(yàn)→將正常飼料和VitE

選入檢驗(yàn)對(duì)列框→檢驗(yàn)類型:

選擇前兩種(Wilcoxon和符號(hào)檢驗(yàn))→確定結(jié)果解釋

首先顯示兩變量正常飼料和VitE缺乏的例數(shù)、均數(shù)、標(biāo)準(zhǔn)差、最大值和最小值；

配對(duì)符號(hào)秩和檢驗(yàn)（WilcoxonMatched-PairsSigned-RanksTest）結(jié)果，其平均秩分別為5.00和1.00，Z=-2.380，雙側(cè)P=0.017，可認(rèn)為兩組大鼠肝中VitA含量有差別，飼料中缺乏VitE會(huì)使大鼠肝中VitA含量降低；

符號(hào)檢驗(yàn)（SignTest）的結(jié)果，雙側(cè)P=0.0703，則認(rèn)為兩組大鼠肝中VitA含量無(wú)差別。

在這種情況下，應(yīng)取配對(duì)符號(hào)秩和檢驗(yàn)（Wilcoxon）結(jié)果，因兩法比較之下，配對(duì)符號(hào)秩和檢驗(yàn)較為敏感，效率較高。四、隨機(jī)區(qū)組設(shè)計(jì)多個(gè)樣本的秩和檢驗(yàn)

8名受試對(duì)象在相同條件下分別接受4種不同頻率聲音的刺激，他們的反應(yīng)率（%）見(jiàn)下表。問(wèn)4種不同頻率聲音的刺激的反應(yīng)率（%）是否有差別？

[例4]8名受事對(duì)象對(duì)4種不同頻率聲音刺激的反應(yīng)率（%）比較

受試號(hào)頻率A頻率B頻率C頻率D

18.49.69.811.7211.612.711.81239.49.110.49.849.88.79.91258.388.68.668.69.89.610.678.9910.611.4

87.88.28.510.8

數(shù)據(jù)準(zhǔn)備

定義4個(gè)變量：

1.頻率A2.頻率B3.頻率C4.頻率D

統(tǒng)計(jì)分析

分析→非參數(shù)統(tǒng)計(jì)→K個(gè)相關(guān)樣本檢驗(yàn)→頻率A、頻率B、頻率C、頻率D均進(jìn)入檢驗(yàn)對(duì)列框→檢驗(yàn)類型:選擇Friedman→確定描述性統(tǒng)計(jì)量

N均值標(biāo)準(zhǔn)差極小值極大值頻率A89.101.197.811.6頻率B89.391.47812.7頻率C89.901.088.511.8頻率D810.861.198.612秩

秩均值頻率A1.38頻率B2.00頻率C2.94頻率D3.69檢驗(yàn)統(tǒng)計(jì)量(a)N8卡方15.152df3漸近顯著性0.00169表６兩種療效對(duì)腎小球腎病的療效比較五、等級(jí)資料秩和檢驗(yàn)療效

消炎痛

合劑

無(wú)效

154部分緩解

69基本緩解

45完全緩解

219

[例5]等級(jí)資料的SPSS操作要點(diǎn)：1.定義3個(gè)變量，即療效、療法（或疾?。┖皖l數(shù)。2.在進(jìn)行統(tǒng)計(jì)分析前，需對(duì)頻數(shù)加權(quán)。3.采用2個(gè)和K個(gè)獨(dú)立樣本的非參數(shù)檢驗(yàn)。

[例6]某村發(fā)生一種地方病，其住戶沿一條河排列，調(diào)查時(shí)對(duì)發(fā)病的住戶標(biāo)記為“1”，共9例，對(duì)非發(fā)病的住戶標(biāo)記為“0”，共17戶：?jiǎn)柌舻姆植寂帕惺浅示奂厔?shì)，還是隨機(jī)分布？01100010010000110010000101六、Runs過(guò)程（游程檢驗(yàn)）主要功能

依時(shí)間或其他順序排列的有序數(shù)列中，具有相同的事件或符號(hào)的連續(xù)部分稱為一個(gè)游程。調(diào)用Runs過(guò)程可進(jìn)行游程檢驗(yàn)，即用于檢驗(yàn)序列中事件發(fā)生過(guò)程的隨機(jī)性分析。假定下面是由0和1組成的一個(gè)這種變量的樣本：

01100010010000110010000101

其中相同的0（或相同的1）在一起稱為一個(gè)游程（單獨(dú)的0或1也算）。這個(gè)數(shù)據(jù)中有７個(gè)0組成的游程和７個(gè)1組成的游程。數(shù)據(jù)準(zhǔn)備

定義1個(gè)變量為epi。

按住戶順序輸入數(shù)據(jù)，發(fā)病的住戶為1,非發(fā)病的住戶為0。統(tǒng)計(jì)分析

分析→非參數(shù)統(tǒng)計(jì)→游程檢驗(yàn)→epi選入檢驗(yàn)變量列框

→選擇割點(diǎn)→確定

1、中位數(shù)：中位數(shù)作割點(diǎn)，小于臨界割點(diǎn)的為一類，大于或等于臨界割點(diǎn)的為另一類；

2、眾數(shù)：眾數(shù)作割點(diǎn)，小于臨界割點(diǎn)的為一類，大于或等于臨界割點(diǎn)的為另一類；

3、均數(shù)：均數(shù)作臨界割點(diǎn)，小于臨界割點(diǎn)的為一類，大于或等于臨界割點(diǎn)的為另一類；

4、定制：用戶指定割點(diǎn)，小于定制割點(diǎn)的為一類，大于或等于割點(diǎn)的為另一類；

本例選定制，在其方框中鍵入1（根據(jù)需要選項(xiàng)，本例是0、1二分變量，故臨界割點(diǎn)值用1），再點(diǎn)擊OK鈕即可。結(jié)果解釋

檢驗(yàn)結(jié)果可見(jiàn)本例游程個(gè)數(shù)為14，檢驗(yàn)臨界割點(diǎn)值（Testvalue）=1.00，

Z=0.325，雙側(cè)P=0.746。所以認(rèn)為此地方病的病戶沿河分布的情況無(wú)聚集性，而是呈隨機(jī)分布。1-SampleK-S過(guò)程

調(diào)用此過(guò)程可對(duì)單樣本進(jìn)行Kolmogorov-SmirnovZ檢驗(yàn)，它將一個(gè)變量的實(shí)際頻數(shù)分布與正態(tài)分布（Normal）、均勻分布（Uniform）、泊松分布（Poisson）進(jìn)行比較。數(shù)據(jù)準(zhǔn)備

以SPSS自帶文件cars．Sav為例，試分析馬力（horse）和加速度（accel）的數(shù)據(jù)是否服從于正態(tài)分布。統(tǒng)計(jì)分析

分析→非參數(shù)統(tǒng)計(jì)→單樣本

Kolmogorov-Smirnov檢驗(yàn)→馬力和加速度選入檢驗(yàn)變量列框

→檢驗(yàn)分布：正態(tài)分布→確定結(jié)果解釋練習(xí)一

某醫(yī)生欲比較某新療法與傳統(tǒng)療法治療腎綜合征出血熱患者的降溫效果，將病人隨機(jī)分為兩組，分別以新療法與傳統(tǒng)療法治療，以用藥開(kāi)始的體溫降至正常值所用的時(shí)間（小時(shí)）為療效指標(biāo)，結(jié)果見(jiàn)表1，試比較兩種療法的退熱時(shí)間有無(wú)差別？新療法傳統(tǒng)療法25363040324435483750395639594260466448195240表1新療法與傳統(tǒng)療法治療結(jié)果(小時(shí))練習(xí)二

批次測(cè)定方法甲乙丙丁12.722.463.953.8622.322.424.313.9532.482.224.523.30表2四種方法測(cè)定甘藍(lán)葉核黃素濃度練習(xí)三

某研究者欲研究保健食品對(duì)小鼠抗疲勞作用，將同種屬的小鼠按性別和年齡相同、體重相近配成對(duì)子，共10對(duì)，并將每對(duì)中的兩只小鼠隨機(jī)分到保健食品兩個(gè)不同的劑量組，測(cè)量小鼠負(fù)重游泳時(shí)間（min，負(fù)重5%體重），結(jié)果見(jiàn)表3，問(wèn)不同劑量組的小鼠負(fù)重游泳時(shí)間有無(wú)差別？表3不同劑量組小鼠負(fù)重游泳時(shí)間(min)

小鼠對(duì)子號(hào)中劑量組高劑量組11415.22135.5315144176.55135.561813.5717.510810.210.2910101010.59.51113.86.8123.033.481315.25.51416.59練習(xí)四

某地監(jiān)測(cè)大氣中SO2的日均濃度，按不同功能區(qū)設(shè)置采樣點(diǎn)，結(jié)果見(jiàn)下表，問(wèn)各功能區(qū)SO2的日均濃度有無(wú)差別？對(duì)照區(qū)工業(yè)區(qū)商業(yè)區(qū)居民區(qū)1046723133830665501352307096304854080266951151851677630表4

不同功能區(qū)大氣中SO2的日均濃度練習(xí)五

四種病人痰液內(nèi)嗜酸性粒細(xì)胞的檢查結(jié)果見(jiàn)表7。問(wèn)四種疾病患者痰液內(nèi)嗜酸粒細(xì)胞的等級(jí)有無(wú)差別?白細(xì)胞等級(jí)支氣管擴(kuò)張肺水腫肺癌病毒性呼吸道感染－1353＋