統(tǒng)計(jì)第13次作業(yè)

1、統(tǒng)計(jì)第13次作業(yè) 郭曉蘭 微生物學(xué) 12213641三、計(jì)算題1. 為了評(píng)價(jià)某試驗(yàn)藥物與對(duì)照藥物對(duì)提高心率變異性（HRV）的療效，按時(shí)域的HRV低于50 ms的標(biāo)準(zhǔn)收治患者，共收治20名HRV異常的患者，隨機(jī)分為試驗(yàn)組和對(duì)照組，試驗(yàn)組服用試驗(yàn)藥，對(duì)照組服用對(duì)照藥。對(duì)每一患者在治療前以及治療后1周、2周、3周分別測(cè)量一次時(shí)域的HRV水平，測(cè)量結(jié)果如教材表17-19，請(qǐng)作統(tǒng)計(jì)分析。教材表17-19 試驗(yàn)組與對(duì)照組治療前后HRV水平測(cè)量結(jié)果試驗(yàn)組HRV測(cè)量值/ms對(duì)照組HRV測(cè)量值/ms編號(hào)治療前治療后編號(hào)治療前治療后第1周第2周第3周第1周第2周第3周145.352.461.070.4136.14

2、1.848.363.0240.348.754.966.4242.046.853.867.4341.048.854.868.8343.750.553.269.2438.044.952.763.0441.446.254.364.7542.750.856.770.2541.550.053.665.7639.145.155.364.9645.052.457.570.6743.152.356.569.1742.247.853.270.3839.346.256.967.8842.545.455.868.7941.948.057.365.4937.042.748.161.91037.844.754.165.

3、01035.540.845.959.5分析：本題的處理因素有兩個(gè)水平，即試驗(yàn)組和對(duì)照組；本題的時(shí)間因素有4個(gè)水平，分別是治療前、治療后第1周、第2周和第3周，屬于典型的重復(fù)測(cè)量設(shè)計(jì)，應(yīng)該用重復(fù)測(cè)量設(shè)計(jì)的統(tǒng)計(jì)分析方法。解：（1）統(tǒng)計(jì)描述表1 試驗(yàn)組與對(duì)照組HRV的平均水平組別基線第1周第2周第3周試驗(yàn)組40.85±2.4248.19±2.9656.02±2.2667.10±2.52對(duì)照組40.69±3.2946.44±3.8752.37±3.7066.10±3.78由SPSS軟件可得輪廓圖，如圖1：圖1 輪廓圖（2）

4、統(tǒng)計(jì)推斷建立檢驗(yàn)假設(shè)，確定檢驗(yàn)水準(zhǔn)：對(duì)于處理因素：處理主效應(yīng)相等。：處理主效應(yīng)不相等。對(duì)于時(shí)間因素：時(shí)間主效應(yīng)相等。：時(shí)間主效應(yīng)不全相等。對(duì)于處理因素與時(shí)間因素的交互效應(yīng)：在4個(gè)時(shí)間點(diǎn)兩種處理的總體均數(shù)之差是相等的，即處理因素與時(shí)間因素?zé)o交互效應(yīng)。：在4個(gè)時(shí)間點(diǎn)兩種處理的總體均數(shù)之差不全相等，即處理因素與時(shí)間因素存在交互效應(yīng)。均??；假設(shè)資料符合重復(fù)測(cè)量設(shè)計(jì)分析的前提條件，用SPSS軟件進(jìn)行分析，步驟如下：AnalyzeGeneral Linear modelRepeated Measurename: time Number of Levels:4 ADDDefine： Within-subje

5、cts variables: 治療前（1），治療后1周（2），治療后2周（3），治療后3周（4）Model ：full factorialGroup:Between-Subjects Factor(s) OK借助spss對(duì)本例兩組各個(gè)時(shí)點(diǎn)的HRV觀察資料進(jìn)行球形檢驗(yàn)（Mauchly's test of sphericity），=0.1，得到球形檢驗(yàn)的=0.005<0.1，所以認(rèn)為資料不滿足球形對(duì)稱性條件，故選用Greenhouse-Geisser。方差分析結(jié)果如表2：表2 兩組4次HRV觀察資料的方差分析表方差來(lái)源處理主效應(yīng)53.7921537921.5110.235時(shí)間主效應(yīng)7

6、316.9312.1973330.1401675.9960.000處理×時(shí)間33.2612.19715.1387.6190.001處理間誤差640.9211835.607重復(fù)測(cè)量誤差78.58339.5491.987合計(jì)8123.48862.9433436.664由表可知：檢驗(yàn)處理因素與時(shí)間因素之間的交互效應(yīng)的0.05，按水準(zhǔn)拒絕，接受，可認(rèn)為處理因素與時(shí)間因素存在交互效應(yīng)。處理主效應(yīng)的0.05，不拒絕,所以處理間的差異沒(méi)有統(tǒng)計(jì)學(xué)意義。時(shí)間主效應(yīng)0.05，拒絕，接受，可以認(rèn)為各個(gè)時(shí)間點(diǎn)的總體均數(shù)不全相等。由于存在交互效應(yīng)，各個(gè)時(shí)間點(diǎn)的兩兩比較比較復(fù)雜，所以可采用作圖作趨勢(shì)分析，得到大

7、致的統(tǒng)計(jì)推斷。2. 為了比較某治療婦女痛經(jīng)藥物的療效與對(duì)照藥效應(yīng)的差異，采用量表Vas對(duì)服藥后60 min時(shí)的痛經(jīng)進(jìn)行評(píng)分。痛經(jīng)評(píng)分越高說(shuō)明痛經(jīng)越嚴(yán)重，反之痛經(jīng)評(píng)分越低說(shuō)明痛經(jīng)程度越低，0分為無(wú)痛經(jīng)。某醫(yī)生共收治30名經(jīng)常痛經(jīng)的女性，隨機(jī)分為兩組，第1組受試者在接受治療第1個(gè)月經(jīng)期感到痛經(jīng)時(shí)服用試驗(yàn)藥，在接受治療的第2個(gè)月經(jīng)期感到痛經(jīng)時(shí)服用對(duì)照藥；第2組受試者在接受治療第1個(gè)月經(jīng)期感到痛經(jīng)時(shí)服用對(duì)照藥，在接受治療的第2個(gè)月經(jīng)期感到痛經(jīng)時(shí)服用試驗(yàn)藥。用雙盲的方法測(cè)定痛經(jīng)評(píng)分，其測(cè)量結(jié)果如教材表17-20，請(qǐng)作統(tǒng)計(jì)分析。教材表17-20 雙盲法測(cè)定疼經(jīng)的評(píng)分結(jié)果第1組痛經(jīng)評(píng)分第2組痛經(jīng)評(píng)分受試者編

8、號(hào)第1個(gè)痛經(jīng)期第2個(gè)痛經(jīng)期受試者編號(hào)第1個(gè)痛經(jīng)期第2個(gè)痛經(jīng)期143341284723431237313483132936434274343952513526366432563342737307354683422834519393093443104222104737113836112242124126124043132832133046143527144039154229153831分析：本題關(guān)注的不僅是處理效應(yīng)，還有階段效應(yīng)和順序效應(yīng)，應(yīng)該采用交叉設(shè)計(jì)分析方法。解：（1）統(tǒng)計(jì)描述表3 Vas分統(tǒng)計(jì)描述（mean±S） 順序132.6±7.65階段135.67±6.

9、36試驗(yàn)藥39.07±6.03順序237.2±6.87階段234.13±8.66對(duì)照藥30.73±6.66表4 Vas分統(tǒng)計(jì)描述（mean±S）時(shí)期1時(shí)期2順序1試驗(yàn)藥：37.53±6.07對(duì)照藥：27.67±5.69順序2對(duì)照藥：33.80±6.28試驗(yàn)藥：40.60±5.78（2）統(tǒng)計(jì)推斷建立檢驗(yàn)假設(shè)，確定檢驗(yàn)水準(zhǔn)：無(wú)效假設(shè)：兩種藥物治療痛經(jīng)的平均痛經(jīng)評(píng)分相同。備擇假設(shè)：兩種藥物治療痛經(jīng)的平均痛經(jīng)評(píng)分不同。無(wú)效假設(shè)：受試者兩個(gè)階段的平均痛經(jīng)評(píng)分相同。備擇假設(shè)：受試者兩個(gè)階段的平均痛經(jīng)評(píng)分不同。無(wú)效假設(shè)

10、：兩種順序受試者的平均痛經(jīng)評(píng)分相同。備擇假設(shè)：兩種順序受試者的平均痛經(jīng)評(píng)分不同。=0.05用SPSS軟件進(jìn)行交叉設(shè)計(jì)方差分析，步驟如下：Analyze General Linear Model Univariate y：Dependent Variable ，period：Fixed Factor(s) ，seq：Fixed Factor(s)，treat：Fixed Factor(s) ，id：Random Factor(s) Model：Custom， Build term(s):Main Effect ，period：Model ，seq：Model，treat：ModelBuild t

11、erms(s):interactionid,seq：modelcontinueoptionsdescriptive statisticsOK表5 交叉設(shè)計(jì)的方差分析變異來(lái)源順序 317.400 1 317.400 8.360.007 3時(shí)期 35.267 1 35.267 1.070.310 6藥物1 041.667 11 041.66731.500.000 0個(gè)體間變異1 063.00028 37.964處理內(nèi)誤差 926.06728 33.074合計(jì)3 383.4005957.346從表可知：對(duì)于無(wú)效假設(shè)（兩種藥物治療痛經(jīng)的平均痛經(jīng)評(píng)分相同）而言，檢驗(yàn)處理效應(yīng)的0.05，故可拒絕，接受，差異有統(tǒng)計(jì)學(xué)意義，可以認(rèn)為兩種藥物治療痛經(jīng)的療效有不同。的95%置信區(qū)間為（5.292,11.375）分，其下限明顯大于0，所以可認(rèn)為試驗(yàn)藥物治療痛經(jīng)效果差于對(duì)照藥物。對(duì)無(wú)效假設(shè)（受試者兩個(gè)階段的平均痛經(jīng)評(píng)分相同）而言，檢驗(yàn)階段效應(yīng)的=0.3110.05，故不拒絕，尚不能認(rèn)為受試者兩個(gè)階段的平均痛經(jīng)評(píng)分有不同。對(duì)于無(wú)效假設(shè)（兩種順序受試者的平均痛經(jīng)評(píng)分相同）而言，檢驗(yàn)順序效應(yīng)的=0.0070.05，故可拒絕，接受，差異有統(tǒng)計(jì)學(xué)意義，可以認(rèn)為受試