05t檢驗(yàn)與方差分析

數(shù)據(jù)的一般描述性統(tǒng)計(jì)分析描述性統(tǒng)計(jì)---DescriptiveStatistics一、Frequencies過程Frequencies過程專門為產(chǎn)生頻數(shù)表而設(shè)計(jì)的.

某地101例健康男子血清總膽固醇值測(cè)定結(jié)果如下4.773.376.143.953.564.234.314.715.694.124.564.375.396.305.217.225.543.935.214.125.185.774.795.125.205.104.704.743.504.694.384.896.255.324.504.633.614.444.434.254.035.854.093.354.084.795.304.973.183.975.165.105.864.795.344.244.324.776.366.384.885.553.044.553.354.874.175.855.165.094.524.384.314.585.726.554.764.614.174.034.473.403.912.704.604.095.965.484.404.555.383.894.604.473.644.345.186.143.244.903.05

Frequencies過程Format對(duì)話框用于定義輸出頻數(shù)表的格式根據(jù)頻數(shù)的多少作頻數(shù)分布選擇兩個(gè)以上變量作頻數(shù)表限制頻數(shù)表的分組數(shù)

輸出結(jié)果（略）Descriptives過程標(biāo)準(zhǔn)化值Options對(duì)話框輸出結(jié)果Explore過程對(duì)數(shù)據(jù)進(jìn)行初步檢查,判斷有無(wú)離群點(diǎn)或極端值;對(duì)前提條件假定,如正態(tài)分布和方差齊性進(jìn)行檢驗(yàn),不滿足正態(tài)分布和方差齊性時(shí),提示數(shù)據(jù)轉(zhuǎn)換方法,然后決定使用參數(shù)方法或非參數(shù)方法;了解組間差異的特征.Explore過程選入分組變量Statistics對(duì)話框顯示集中趨勢(shì)、離散趨勢(shì)及分布形狀的統(tǒng)計(jì)量顯示總體均數(shù)的可信區(qū)間顯示5個(gè)最高和最低的變量值計(jì)算M估計(jì)量Plots對(duì)話框方差齊性檢驗(yàn)Options對(duì)話框剔除含有缺失值的全部記錄只剔除此分析中含有缺失值的記錄報(bào)告缺失值輸出結(jié)果（略）T檢驗(yàn)一、統(tǒng)計(jì)方法復(fù)習(xí)t檢驗(yàn)的應(yīng)用條件獨(dú)立的隨機(jī)樣本資料服從正態(tài)分布方差齊性可比性

資料形式單樣本資料：樣本均數(shù)與總體均數(shù)比較成組設(shè)計(jì)的兩樣本資料：兩樣本均數(shù)比較配對(duì)資料:差值均數(shù)d與總體均數(shù)0比較二、CompareMeans---均數(shù)比較Means過程Options對(duì)話框?qū)Ψ纸M變量進(jìn)行單因素方差分析，并計(jì)算Eta統(tǒng)計(jì)量線性趨勢(shì)檢驗(yàn)結(jié)果輸出TTest

例正常成年男子血小板計(jì)數(shù)的均數(shù)為225109/L，今測(cè)得20名男性油漆作業(yè)工人的血小板計(jì)數(shù)值(109/L)：220，188，162，230，145，160，238，188，247，113，126，245，164，231，256，183，190，158，224，175。問油漆作業(yè)工人的血小板計(jì)數(shù)與正常成年男子有無(wú)差異？One-SampleTTest過程填入已知總體均數(shù)統(tǒng)計(jì)學(xué)結(jié)論與專業(yè)推斷p>α，差別沒有統(tǒng)計(jì)學(xué)意義（或差異無(wú)顯著性），不拒絕H0（不等于接受H0

），現(xiàn)有樣本暫(尚）不能認(rèn)為油漆作業(yè)工人的血小板計(jì)數(shù)與正常成年男子有差異。p≤α，差別有統(tǒng)計(jì)學(xué)意義（或差異有顯著性），拒絕H0，接受H1，認(rèn)為油漆作業(yè)工人的血小板計(jì)數(shù)與正常成年男子有差異。。。。。。例某克山病區(qū)測(cè)得11例克山病患者與13名健康人的血磷值(mmol/L)如下,問該地急性克山病患者與健康人的血磷值是否不同？患者:0.841.051.201.201.391.531.671.801.872.072.11健康人:0.540.640.640.750.760.811.161.201.341.351.481.561.87TTestIndependent-SampleTTest過程

結(jié)果輸出(略)例1下表是某疾病在醫(yī)改前后12家醫(yī)院的平均費(fèi)用，問改革前后醫(yī)療費(fèi)用有無(wú)差別?

TTest表某疾病醫(yī)改前和醫(yī)改后12家醫(yī)院的平均費(fèi)用被測(cè)者號(hào)(1)醫(yī)改后(2)醫(yī)改前(3)差值d(4)=(3)-(2)d2(5)11490152535122521397141518324315121508-41641401144443184951470150030900614151460452025714311390-41168181429143239914201420001012751227-482304111165126810310609121421144322484合計(jì)20621426練習(xí)

某醫(yī)院用某新藥與常規(guī)藥物治療嬰幼兒貧血,將20名貧血患兒隨機(jī)等分為2組,分別接受兩種藥物治療，測(cè)得血紅蛋白增加量(g/L)如下，問新藥與常規(guī)藥物的療效有無(wú)差別?新藥組24

36

25

14

26

34

23

20

15

19

常規(guī)藥物組14

1820

15

22

24

21

25

27

23

不同地區(qū)農(nóng)村合作醫(yī)療自負(fù)部分（元）蘇中蘇北5040201065355525202015158030353010525757030656045455060353015問該法測(cè)得的CaCO3含量是否偏高？20.9920.4120.1020.0020.9122.60已知某水樣中含CaCO3的真值為20.7mg/L，現(xiàn)用某法重復(fù)測(cè)定該水樣11次，CaCO3的含量（mg/L）為：20.4120.0023.0022.0020.99

例20例某病患者隨機(jī)分為兩組,分別用甲、乙兩藥治療，測(cè)得治療前后（治療后1個(gè)月）的血沉（mm/h），問甲、乙兩藥是否有效？·拆分文件Paired-SampleTTest過程輸出結(jié)果（略）三個(gè)不同地區(qū)農(nóng)村合作醫(yī)療自負(fù)部分（元）蘇南蘇中蘇北605040302010100653585552520202055151545803030353075105251057570306560454550603530151.問題的提出t檢驗(yàn)的局限性單因素兩水平問題的提出t檢驗(yàn)實(shí)例南京中醫(yī)藥大學(xué)的學(xué)工部門想知道2008年南京市進(jìn)入南中醫(yī)就讀學(xué)生中，市區(qū)的學(xué)生和郊區(qū)的學(xué)生成績(jī)是否有差異。因素：學(xué)生所來(lái)自的區(qū)域水平：市區(qū)，郊區(qū)單因素兩水平問題的提出江蘇、安徽、浙江三省的平均入學(xué)成績(jī)？單因素三水平安徽X=571.23江蘇X=592.79浙江X=569.83方差分析（AnalysisofVariance,ANOVA）2．ANOVA優(yōu)點(diǎn)①

不受比較的組數(shù)限制②

可同時(shí)分析多個(gè)因素的作用③

可分析因素間的交互作用方差分析的概念

方差分析（ANOVA）又稱變異數(shù)分析或F檢驗(yàn)，其目的是推斷兩組或多組資料的總體均數(shù)是否相同，檢驗(yàn)兩個(gè)或多個(gè)樣本均數(shù)的差異是否有統(tǒng)計(jì)學(xué)意義。方差分析的基本思想

某克山病區(qū)測(cè)得11例克山病患者和13名健康人的血磷值（mmol/L）如下，患者：0.841.051.201.201.391.531.671.801.872.072.11健康人：0.540.640.640.750.760.811.161.201.341.351.481.561.87問該地克山病患者與健康人的血磷值是否不同？總變異有以下兩個(gè)來(lái)源：（1）組內(nèi)變異，即由于隨機(jī)誤差的原因使得各組內(nèi)部的血磷值各不相等；（2）組間變異，即由于克山病的影響使得患者與健康人組的血磷值均數(shù)大小不等。一、統(tǒng)計(jì)方法復(fù)習(xí)單因素方差分析雙因素方差分析方差分析表變異來(lái)源SSvMSF組間SS組間k-1SS組間/v組間MS組間MS組內(nèi)組內(nèi)SS組內(nèi)N-kSS組內(nèi)/v組內(nèi)總SS總N-1一．方差分析的意義因素(factors)：將試驗(yàn)對(duì)象隨機(jī)分為若干個(gè)組，加以不同的干預(yù)，稱為處理因素。醫(yī)保對(duì)就診率的影響是否手術(shù)治療對(duì)生存率的影響醫(yī)改對(duì)醫(yī)療費(fèi)用的影響在相同的因素下的不同干預(yù)，稱為不同的水平(level)性別：男、女是否手術(shù)：是、否

例用二氧化硅(SIO2)50mg使大鼠染塵后，測(cè)量不同時(shí)期全肺濕重的變化，試比較染塵后1月、3月、6月三個(gè)時(shí)期的全肺濕重有無(wú)差別？二、單因素方差分析Means過程One-WayANOVA過程趨勢(shì)檢驗(yàn)和精確兩兩比較各組之間的多重比較Contrasts對(duì)話框PostHoc對(duì)話框滿足方差齊性要求Options對(duì)話框方差齊性檢驗(yàn)輸出結(jié)果（略）練習(xí)：某醫(yī)師將20名失眠患者隨機(jī)等分為兩組，一組服用安眠藥，另一組服用安慰劑，研究者要評(píng)價(jià)該安眠藥的催眠效果，試作統(tǒng)計(jì)分析。安眠藥組病例號(hào)12345678910療前3.53.24.32.53.83.73.04.12.83.6療后4.75.65.63.55.64.65.25.34.06.1安慰劑組

病例號(hào)12345678910療前2.42.83.24.03.92.62.93.43.84.2療后3.85.14.55.26.03.95.34.25.46.2練習(xí)：某醫(yī)師將20名失眠患者隨機(jī)等分為兩組，一組服用安眠藥，另一組服用安慰劑，研究者要評(píng)價(jià)該安眠藥的催眠效果，試作統(tǒng)計(jì)分析?！ぐ裁咚幗M病例號(hào)12345678910療前3.53.24.32.53.83.73.04.12.83.6療后4.75.65.63.55.64.65.25.34.06.1差值1.22.41.31.01.90.92.21.21.22.5安慰劑組

病例號(hào)12345678910療前2.42.83.24.03.92.62.93.43.84.2療后3.85.14.55.26.03.95.34.25.46.2差值1.42.31.31.22.11.32.40.81.62.0分析思路研究目的：觀察安眠藥是否有催眠效果主要觀察指標(biāo)：睡眠時(shí)間的改變數(shù)據(jù)類型：數(shù)值變量資料分析思路療前兩組同質(zhì)性比較，以分析可比性；各組療前療后差值比較，分別確定各自的效療；兩組療前療后差值相互比較，分析兩組的效果是否相同？同時(shí)估計(jì)差值總體均數(shù)的可信區(qū)間。兩組藥物治療失眠效果(小時(shí))的比較

療前療后差值安眠藥組

3.450.575.020.81

1.570.60

安慰劑組

3.320.644.960.83

1.640.53

療前兩組比較： t=0.481，P=0.637

可以認(rèn)為兩組具有可比性。

=n1＋n2－2=18t0.05,18=2.101t0.50,18=0.688安眠藥組服藥前后比較：t=8.248，P=0.0000

可以認(rèn)為服藥安眠藥后睡眠時(shí)間延長(zhǎng)。

=n-1=9t0.05,9=2.262t0.001,9=4.781安慰劑組服藥前后比較：t=9.745，P=0.0000

可以認(rèn)為服藥安慰劑后睡眠時(shí)間延長(zhǎng)。

=n-1=9t0.05,9=2.262t0.001,9=4.781兩組服藥前后差值相互比較：t=0.276，P=0.786

尚且不能認(rèn)為兩組延長(zhǎng)睡眠時(shí)間的效果有差別。

=n1＋n2－2=18t0.05,18=2.101t0.50,18=0.688綜合結(jié)論統(tǒng)計(jì)分析結(jié)果表明服藥前兩組具有可比性(t=0.481，P=0.637)；服用安眠藥和安慰劑均有延長(zhǎng)睡眠時(shí)間的效果(t=8.248，P<0.001；t=9.745，P<0.001)；但服用安眠藥與服用安慰劑平均延長(zhǎng)的睡眠時(shí)間無(wú)差別(t=0.276，P=0.786)。即安眠藥無(wú)效。

例20例某病患者隨機(jī)分為兩組,分別用甲、乙兩藥治療，測(cè)得治療前后（治療后1個(gè)月）的血沉（mm/h），問甲、乙兩藥是否有效？·五個(gè)地區(qū)每天發(fā)生交通事故的次數(shù)如下：

試檢驗(yàn)各地區(qū)平均每天交通事故次數(shù)是否相等東部北部中部南部西部1512101413171014912141313791117151014—1412810———79五個(gè)地區(qū)每天發(fā)生交通事故數(shù)據(jù)列表三、單變量方差分析編號(hào)工前工中工后190.53142.1287.38288.43163.1765.27347.3763.1668.434175.8166.33210.545100.01144.75194.75646.32126.3365.27773.69138.96200.028105.27126.33100.01986.32121.06105.271060.0173.6958.95不同時(shí)間段職工尿氮的含量——隨機(jī)區(qū)組設(shè)計(jì)組間因素檢驗(yàn)（Testofbetweensubjectseffects）表明：看處理組和區(qū)組的F和P值，看其是否有統(tǒng)計(jì)學(xué)意義；看處理組間的兩兩比較的結(jié)果。GeneralLinearModel---一般線性模型Univariate

過程最小二乘法權(quán)重系數(shù)選入固定因素Model：指定不同的模型，除方差分析外GeneralLinearModel可作其他統(tǒng)計(jì)分析；Contrasts：指定一種要用t檢驗(yàn)來(lái)檢驗(yàn)的priori對(duì)比；Plots：指定作某種圖；PostHoc：指定兩兩比較的方法；Save：指定將產(chǎn)生的一些指標(biāo)保存為新的變量；Options：指定要輸出的一些選項(xiàng)，如數(shù)據(jù)的描述方差齊性檢等。Model對(duì)話框分析所有的主效應(yīng)和交互作用自定義模型PostHoc對(duì)話框Contrasts對(duì)話框因素內(nèi)水平間差值比較None:不進(jìn)行因素內(nèi)不同水平之間的平均值比較;Deviation:各個(gè)因素每個(gè)水平的平均值都進(jìn)行交互比較;Simple:各個(gè)因素變量的每一水平都與指定的參考水平進(jìn)行比較;Difference:將因素變量的每一水平的平均值,除第一水平以外,都與其前面各水平的平均值進(jìn)行比較;Helmert:將因素變量的每一水平的平均值,除最后一個(gè)水平以外,都與后續(xù)的各水平的平均值進(jìn)行比較;Repeated:只對(duì)相鄰的水平進(jìn)行比較;Polynomial:多項(xiàng)式比較,比較因素變量平均值水平的線性、二次、三次以及高次影響。Plots

、Save、Options對(duì)話框輸出結(jié)果（略）乙藥甲藥用（1）不用（2）用（1）2.10.92.21.12.01.0不用（2）1.30.81.20.91.10.7治療缺鐵性貧血病人12例，分為4組給予不同治療，一個(gè)月后觀察紅細(xì)胞增加數(shù)，問甲藥、乙藥單獨(dú)使用的治療效果如果？同時(shí)使