醫(yī)學(xué)統(tǒng)計學(xué)實習(xí)指導(dǎo)1

1、實習(xí)一計量資料的統(tǒng)計描述、統(tǒng)計圖表一、目的與要求：1、了解計量資料的頻數(shù)分布表的編制方法和分布規(guī)律；2、掌握計量資料的統(tǒng)計描述方法；3、掌握正態(tài)分布的概念和特征，標(biāo)準(zhǔn)正態(tài)分布的概念和標(biāo)準(zhǔn)化變換，正態(tài)分布 的應(yīng)用；4、理解統(tǒng)計表和統(tǒng)計圖在統(tǒng)計描述中的作用；5、掌握統(tǒng)計表的制表原則和基本要求；6、掌握制作統(tǒng)計圖的基本要求和各種統(tǒng)計圖的適用條件。二、學(xué)時：3學(xué)時。三、案例：1某市1974年留駐該市一年以上，無明顯肝、腎疾病，無汞作業(yè)接觸史的 居民238人的發(fā)汞值（pmol/kg ）檢測結(jié)果見表1-1的第1和3列所示: 表1-1某市1974年238人的發(fā)汞檢測結(jié)果組段組中值x 0人數(shù)f頻率累計頻率fx

2、 0fx o 21. 52. 52 00. 0840. 08 45012>3. 54. 5660. 2770. 3612971336. 55. 5 6. 5600.2520. 61339(.)25357. 5 & 5480. 2020. 81540834689. 5 10. 5180. 0760. 8911891984. 511. 5 12. 5160. 0670. 9582 00250013. 5 14. 560. 0250. 983871261. 515. 5 16. 510. 0040. 98716. 5272. 2517. 5 1& 500. 0000. 9870

3、019. 5 21. 520. 530. 0131. 00061. 51260. 75合計2381. 000169914743. 5問題1 ：簡述發(fā)汞的分布規(guī)律。t= = 7.1387問題2 :某人按照加權(quán)法計算的均數(shù)為wmmol/kg ）,是否可以用該均數(shù)說明其發(fā)汞的平均水平？為什么？ 問題3 :計算發(fā)汞的95%醫(yī)學(xué)參考值范圍。2 .根據(jù)1999年某地某單位的體檢資料，116名正常成年女子的血清甘油三 酯（mmol/l ）測量結(jié)果見表1-2 :表1-2某單位1999年正常成年女子血清甘油三酯（mmol/l ）測量結(jié)果組段頻數(shù)0. 6 10. 7 30. 8 90. 9 131. 0 191.

4、 1251. 2 181. 3 131.4 91. 5 51.6 1.71合計116問題1 :描述集中趨勢應(yīng)選擇何指標(biāo)？并計算。問題2 :描述離散趨勢應(yīng)選擇何指標(biāo)？并計算。問題3 :求該地正常成年女子血清甘油三酯的95%參考值范圍。問題4 :試估計該地正常成年女子血清甘油三酯在0.8 mrnol/l以下者及1.5 mmol/l以下者各占正常女子總?cè)藬?shù)的百分比。問題5 :該地90%正常成年女子血清甘油三酯集中在哪個范圍。3 .某單位1990-1998年普通感冒與支氣管炎發(fā)病率見表1-3 :表1-3某單位1990-1998年普通感冒與支氣管炎發(fā)病率（%o ）年份感冒支氣管炎1990126. 276

5、. 37199192. 196. 371992107. 595. 901993101. 935. 69199492. 605. 49199573. 204. 32199651. 403. 04199742. 392. 42199833. 922. 27問題1 :說明兩種疾病隨時間的變化趨勢應(yīng)該選用何種統(tǒng)計圖形？并繪制。 問題2 :說明兩種疾病隨時間的變化速度應(yīng)選擇何種統(tǒng)計圖形？并繪制。4某人制作了表1-4 :表1-4某地eb病毒抗體陽性人群藥物干預(yù)治療后eb病毒抗體陽性例數(shù)第一療程第二療程1月2月3月1月2月3月血清學(xué)指檢血清學(xué)指檢血清學(xué)指檢血清學(xué)指檢血清學(xué)指檢血清學(xué)指標(biāo)查 標(biāo)查 標(biāo)查 標(biāo)查

6、標(biāo)查 標(biāo)v e d人v e d人v e d人v e d人v e d人v e dcan數(shù)can數(shù)can數(shù)can數(shù)can數(shù)canaaaaaaaaaaaasesesesesese328252845231523778463062296228283192562227216844202661443818問題1 :該統(tǒng)計分析表存在哪些問題？問題2 :根據(jù)該表所反映的內(nèi)容制作正確的統(tǒng)計分析表。附：spss程序描述統(tǒng)計量以實習(xí)一第2題為例，求描述統(tǒng)計量，如均數(shù)，中位數(shù)，標(biāo)準(zhǔn) 差，標(biāo)準(zhǔn)誤，最大值，最小值，第2.5、25、50、75、97.5百分位數(shù)；求偏 度和峰度系數(shù)及其標(biāo)準(zhǔn)誤；繪制直方圖。數(shù)據(jù)文件：“例1-2.

7、 sav”。數(shù)據(jù)格式：1列116行，1個反應(yīng)變量，變量名為“甘油三酯”。程序：frequenciesvariables=甘油三酯/ntiles= 4/percentiles= 2. 5 97. 5/statistics=stddev semean mean skewness seskew kurtosis sekurt /histogram normal/order= analysis .附：sas程序描述統(tǒng)計量以下的幾個程序均以實習(xí)一第2題為例。/*制作頻數(shù)表*/data exl-2;/*建立數(shù)據(jù)集*/input x q)q);/*確定變量名稱x*/low=0. 6;/*確定頻數(shù)表最低下限為

8、0.6*/dis=0. 1;/*確定頻數(shù)表組距為0.1*/z=x-mod (x-low,dis);/*新建變量z,該變量的值就是將變量x轉(zhuǎn)化為 該數(shù)值所在組段的下限值*/cards;/*給變量賦值*/略proc freq;/*調(diào)用freq過程*/tables z;/*計算變量z各數(shù)值的頻數(shù)*/run;/*用means過程計算例數(shù)、均數(shù)、標(biāo)準(zhǔn)差、標(biāo)準(zhǔn)誤和95%可信區(qū)間*/proc means data=exl-2/*調(diào)用means it程，處理數(shù)據(jù)集 exl 2. sd2*/n mean std stderr elm;/*計算的統(tǒng)計量為例數(shù)、均數(shù)、標(biāo)準(zhǔn)差、標(biāo) 準(zhǔn)誤和95%可信區(qū)間*/var x;

9、/*分析變量為x*/run;/*用means過程對頻數(shù)表計算例數(shù)、均數(shù)、標(biāo)準(zhǔn)差、最小值和最大值*/data exl-2;/*建立數(shù)據(jù)集*/input x f 6)6);/*確定變量名稱，x的值為各組段的組中值，f的值為各組段的頻數(shù)*/cards;/*變量賦值*/略proc means;/*調(diào)用means過程*/freq f;/*定義頻數(shù)變量為"/var x;/*定義分析變量為x*/run;/*用univariate過程進行描述性統(tǒng)計*/proc univariate data=exl-2;/*調(diào)用 univariate 過程*/freq f/*定義頻數(shù)變量為f*/var x;/*定義

10、分析變量為x （組中值）*/run;/*計算2. 5 %和97. 5 %分位數(shù)*/proc univariate data=exl-2;/*調(diào)用 univariate 過程*/var x;/*定義分析變量為x*/output out=pct/*將結(jié)果輸出到外部數(shù)據(jù)集pct*/pctlpre=p/*數(shù)據(jù)集pct中百分位數(shù)變量的前綴*/pctlpts=2. 5 97. 5;/*需要計算2. 5%和97. 5%百分位數(shù)*/run;proc print data=pct;/*顯示數(shù)據(jù)集pct*/run;八作正態(tài)性檢驗、莖葉圖、箱式圖和正態(tài)概率圖*/proc univariate data=exl.2/

11、*調(diào)用 univariate 過程*/normal/*作正態(tài)性檢驗*/plot;/*作莖葉圖、箱式圖和正態(tài)概率圖*/var x;/*定義分析變量為x*/run;實習(xí)二計量資料的統(tǒng)計分析方法一、目的與要求：1、熟悉抽樣誤差的概念；熟悉標(biāo)準(zhǔn)誤的意義及其應(yīng)用；標(biāo)準(zhǔn)差與均數(shù)標(biāo)準(zhǔn)誤的 區(qū)別；2、掌握總體均數(shù)95 %的可信區(qū)間計算及適用條件；3、掌握常用的t檢驗方法；4、了解方差齊性檢驗和變量變換；5、掌握方差分析的基本思想；各種設(shè)計方案（完全隨機設(shè)計、隨機區(qū)組設(shè)計） 變異和自由度的分解方法；6、掌握不同設(shè)計非參數(shù)統(tǒng)計假設(shè)檢驗方法和步驟。二、學(xué)時：3學(xué)時。三、案例：1某地隨機抽樣調(diào)查了部分健康成人的紅細胞

12、數(shù)和血紅蛋白量，結(jié)果見表2-1 :表2-1某年某地健康成年人的紅細胞數(shù)和血紅蛋白含量指標(biāo)性別例數(shù)均數(shù)標(biāo)準(zhǔn)差標(biāo)準(zhǔn)值*紅細胞數(shù)男3604. 660. 584. 84(10 /l )女2匸4. 180. 294. 33血紅蛋白男360134. 57. 1140. 2(g/l )女255117. 610. 2124. 7*實用內(nèi)科學(xué)（1976年）所載均數(shù)（轉(zhuǎn)為法定單位）問題1 :說明女性的紅細胞數(shù)與血紅蛋白的變異程度何者為大？問題2 :計算男性兩項指標(biāo)的抽樣誤差。問題3 :估計該地健康成年女性紅細胞數(shù)的均數(shù)。問題4 :該地健康成年男、女血紅蛋白含量是否不同？問題5 :該地男性兩項血液指標(biāo)是否均低于上表

13、的標(biāo)準(zhǔn)值（若測定方法相 同）？2. 某醫(yī)生應(yīng)用潑尼松、轉(zhuǎn)移因子和胸腺肽治療系統(tǒng)性紅斑狼瘡（sle ）患者14 人。治療前后血清sil-2r （ u/ml ）數(shù)據(jù)見表2-2 :表2-2治療前后血清sil-2r （ u/ml ）數(shù)據(jù)no.1234567治療前1410. 37893. 541569. 45936. 51529. 94477. 23999. 4治療后1353. 57876. 881534.42879. 9846& 66427. 23971. 56no.891011121314治療前474. 8587 3. 04252. 611227. 20595. 4 0359. xi1097

14、. 99治療后446. 67825. 06175. 301110. 19470. 83337. 751022. 31問題1 :該醫(yī)生對此數(shù)據(jù)應(yīng)用兩組獨立樣本的f檢驗，結(jié)果為t =0.3737 , 自由度為26 ,7=0.7116。于是，該醫(yī)生的結(jié)論是治療前后血清sil-2r的差異沒有統(tǒng)計學(xué)意義。你是否同意這種分析結(jié)果？問題2 :有人提議做配對資料的t檢驗，如果治療前后的差異有統(tǒng)計學(xué)意義就 可以說明治療有效。你是否同意這樣做？3. 腦積液磷酸己糖異構(gòu)酶檢測用于腦膜炎診斷的探討一文為比較三組患兒 cst中phi值是否不同，數(shù)據(jù)及分析結(jié)果見表2-3 o表2-3三組患兒cst中phi值的比較組別刀tp

15、pm15407. 0 ± 294. 75. 34<0. 01wm、 ve1415. 0± 13. 16. 47<0. 01對照組237. 0±4. 8問題1 ：該資料采用的是何種統(tǒng)計分析方法?問題2 :使用的統(tǒng)計分析方法是否正確？若不正確，可以采用何種正確的統(tǒng)計 分析方法。問題3 :采用該統(tǒng)計分析方法應(yīng)滿足什么條件？該資料是否滿足？4. 某研究者對惡性滋養(yǎng)細胞腫瘤患者灌注治療，測得其治療前后的hcg （pmol/l ）值，其結(jié)果見表2-4 :表2-4惡性滋養(yǎng)細胞腫瘤患者灌注前后hcg值（pmol/l ）病例號治療前（*1 ）治療后（xi ）log x

16、1log x 211280000210000. 06. 107215. 322222755003300. 04.877953. 518513124502210. 04.095173. 34439415000006. 96. 176090. 838855100002500. 04. 000003. 39794697001203. 03. 986773. 080277155884825. 04. 192793. 6835084223914. 03. 625622. 96095對于該資料研究者作了如下分析：分別對 禺和 乂取對數(shù)，血-比竝衍 斥=-2. 3733以 一7查£界值表， j 皇

17、茨*gsm<005。得出治療有效的結(jié)論。問題1 ：該資料的分析方法是否合適？為什么？問題2 :應(yīng)該如何分析該資料？問題3 :從該研究者的分析結(jié)果您還應(yīng)該注意什么問題？ 附：spss程序1. 配對檢驗，以實習(xí)二第2題為例。數(shù)據(jù)文件：“例2-2. sav”。數(shù)據(jù)格式：2列14行，2個反應(yīng)變量，變量名為“ xl ”和“ x2 ”。 程序：t-testpairs = xl with x2 （paired）/criteria = ci （. 95）/missing = analysis.2兩樣本均數(shù)比較的（檢驗，以醫(yī)學(xué)統(tǒng)計學(xué)教材例題3-7為例。數(shù)據(jù)文件：“例 3-7. sav ”。數(shù)據(jù)格式：2歹寸

18、40行，1個反應(yīng)變量，變量名為“” ；1個分組變量，變量名為“ group ”，有2個水平。 程序：t-testgroups = group (1 2)/missing = analysis/variables = x/criteria = ci (. 95)3完全隨機設(shè)計資料的方差分析，以實習(xí)二第3題為例。 數(shù)據(jù)文件：“例2-3. sav ”。數(shù)據(jù)格式：2列52行，1個反應(yīng)變量，變量名為“phi值”；1個分組變量, 變量名為“ group ”，有3個水平。程序：onewayphi 值 by group/statistics descriptives homogeneity/plot mean

19、s/missing analysis/posthoc = snk lsd alpha (. 05).附：sas程序1. 配對資料兩樣本均數(shù)比較的(檢驗以實習(xí)二第2題為例。 /*用means作配對資料兩個樣本均數(shù)比較的t檢驗*/data ex2 2;/*建立數(shù)據(jù)集*/input xl x2 5)5);/*確定變量名稱，xl和x2分別為兩組個體值*/d=xl-x2;/*計算兩個變量的差值作為變量d的值*/cards;/*變量賦值*/略proc means t prt;/*調(diào)用means it程，要求輸出1值和所對應(yīng)的/ 值*/var d;/*定義分析變量為d*/run;/*用univariate過

20、程作配對資料兩樣本均數(shù)比較的t檢驗*/proc univariate data=ex2 2 ;/*調(diào)用 univariate 過程*/var d ;/*定義分析變量為d*/run ;2兩樣本均數(shù)比較的/檢驗以醫(yī)學(xué)統(tǒng)計學(xué)教材例題3-7為例data ex3-7;/*建立數(shù)據(jù)集*/input x cfflffl;/*確定變量名稱，x為分析變量，c為分組變量*/cards;/*變量賦值*/略proc (test;/*調(diào)用ties!過程*/var x;/*定義分析變量為x*/class c;/*定義分組變量為c*/run;3完全隨機設(shè)計資料的方差分析以實習(xí)二第3題為例。data ex2 3;/*建立數(shù)據(jù)

21、集*/input x c;/*確定變量名稱，x為分析變量，c為分組變量*/cards;/*變量賦值*/略proc anova;/*調(diào)用anova過程*/class c;/*定義分組變量為c*/model x=c;/*定義模型，分析g對x的影響*/means c/lsd;/*用lsd法對多組均數(shù)進行兩兩比較*/means c/hovtest;/*作方差齊性檢驗，默認levene法*/run;實習(xí)三計數(shù)資料與等級資料的統(tǒng)計分析方法一、目的與要求：1、掌握常用相對數(shù)的區(qū)別和易犯的錯誤；2、掌握使用標(biāo)準(zhǔn)化法的目的，標(biāo)準(zhǔn)化率的計算方法和應(yīng)用條件；3、掌握等級資料秩和檢驗的方法；4、了解poisson分布

22、的基本概念及方法；5、掌握尸檢驗的基本思想及方法。二、學(xué)時：3學(xué)時。三、案例：1 某醫(yī)生在一個有5萬人口的社區(qū)進行肺癌調(diào)查，通過隨機抽樣共調(diào)查2000人，全部調(diào)查 工作在10天內(nèi)完成，調(diào)查內(nèi)容包括流行病學(xué)資料和臨床實驗室檢查資料。調(diào)查結(jié)果見表3-1 :表3-1某社區(qū)不同性別人群肺癌情況性別檢查人數(shù)有病人數(shù)死亡人數(shù)死亡率（%發(fā)病率（% ）男10506350. 00. 6女9503266. 70. 3合計20009558. 40.45該醫(yī)生對表中資料進行了統(tǒng)計分析，認為男性肺癌的發(fā)病率高于女性，而死亡情況則完全相反。 問題1 :該醫(yī)生所選擇的統(tǒng)計指標(biāo)正確嗎？問題2 :該醫(yī)生對指標(biāo)的計算方法恰當(dāng)嗎？

23、問題3 :應(yīng)該如何作適當(dāng)?shù)慕y(tǒng)計分析？2 .某醫(yī)生統(tǒng)計了本醫(yī)院從i960年至1999年的住院病例，得到表3-2和表3-3的數(shù)據(jù)。表3-2 1960-99年住院患者中心臟病患者的頻率年度住院患者中10種心臟病患者的頻 率（%）1960-4. 61970-7. 31980-8. 91990-200010. 1表3-3 1960-99年住院心臟病患者中冠心病患者的頻率年度10種心臟病住院患者中冠心病患者的頻率 (%)1960-9. 381 970-30. 101980-29. 211990-200051. 45根據(jù)表3-2 ,該醫(yī)生認為“40年來，10種心臟病總的發(fā)病率逐年代增高”;根據(jù)表3-3 ,

24、該醫(yī)生認為“在10種心臟病中，冠心病發(fā)病率迅速增高，從二十世紀(jì)70年代起居心臟病中 第一位”。問題1 :如何評價這位醫(yī)生的分析結(jié)論？ 問題2 :這兩份資料各有什么用處？3a、 b兩個條件相近的縣級醫(yī)院院長在一起談?wù)摳髯葬t(yī)院的治愈率時得到見表3-4 : 表3 -4 a、 b兩個縣級醫(yī)院的治愈率科室a醫(yī)院b醫(yī)院治療人數(shù)治愈率治療人數(shù)治愈率內(nèi)科8306511()80外科5506316368兒科3153952639婦產(chǎn)科4293339535中醫(yī)科1862837530其他科室36 3626810合計267346183737a醫(yī)院院長自豪地說：“盡管都是縣級醫(yī)院，硬件條件也差不多，但本院的治愈率46%明顯

25、 高于貴醫(yī)院的37%，可見我們的醫(yī)療質(zhì)量和管理水平較高”。b醫(yī)院院長不服氣，說道：“我 們醫(yī)院的各個治療科室的治愈率均比你們高，應(yīng)該是我們醫(yī)院的醫(yī)療質(zhì)量和管理水平高”。二 人各不相讓，爭執(zhí)不下。問題1 :為什么同樣的資料會產(chǎn)生完全不同的看法，根本原因何在？如何才能避免類似的爭執(zhí)？ 問題2 :你認為哪個的結(jié)論正確，或者都不正確？要怎樣才能做出合理的評判？4 .2000年某地愛滋病病毒的感染率為十萬分之七，該地10萬人口， 2001年感染了愛滋病病毒的人數(shù)為17人，有人說，該地2001年總體上愛滋病病毒感染率與2000年持平。如果 是這樣的話，該地2001年感染了愛滋病病毒的人數(shù)為17人這種情況發(fā)

26、生的概率為pcjc = 17) = <? = 0.0006x '171因為發(fā)生的概率太小了，所以說該地2001年總體上愛滋病病毒感染率與2000年持平的說法是 不成立的。問題1 :該分析是否正確？問題2 :如果有問題，出在哪里？5. 某地對區(qū)級醫(yī)院20012002年醫(yī)療質(zhì)量進行總體評價與比較，按分層抽樣方法抽取兩年內(nèi)某 病患者1250例，患者年齡構(gòu)成與病情兩年間差別沒有統(tǒng)計學(xué)意義，觀察三項指標(biāo)分別為療效、住院日、費用。規(guī)定很好、好、一般、差的標(biāo)準(zhǔn)見表3-5 ,病人醫(yī)療質(zhì)量各等級頻數(shù)分布見表3-6 : 表3-5很好、好、一般、差的標(biāo)準(zhǔn)指標(biāo)很好好一般差療效治愈顯效好轉(zhuǎn)無效住院日1 6

27、2 02廠25> 25費用（元）<1400140018001800 2200> 2200表3-6兩年病人按醫(yī)療質(zhì)量等級的頻數(shù)分配表標(biāo) 指好 量 質(zhì)差效 療年年1 2 o oo803o o2 160o日 院 住年年1 2 o ooo2o8131050212013020o4用年年1 2 o oo1111o31o23o27o2o3o o11o7對療效、住院日、費用三項指標(biāo)分別采用尸檢驗，結(jié)果為療效 尸=6. 786尸=0. 079住院日=12. 568尸=0006費用尸=9. 613尸=0. 020故不能認為兩年療效不同；而兩年的住院日和費用均有差別，根據(jù)調(diào)查所得平均住院日和平均費

28、 用可以認為平均住院日2001年比2002年長，而費用2001年低于2002年。問題1 :以上f檢驗方法是否正確？為什么？問題2 :如果不正確，問題出在什么地方？該選擇何種假設(shè)檢驗方法？6. 某研究者欲比較某藥物與對照藥的療效，將類似病情的患者隨機分成兩組，分別接受該藥物 和對照藥，結(jié)果見表2-7 :表2-7某藥物治療某病人療效結(jié)果療效治療組對照組合計惡化112無變化51419進步131023顯著進步9514基本痊愈202合計303060對于該資料研究者作了列聯(lián)表尸檢驗:5. 79 ,u-4 ,尸a at在m = q.«的檢驗水準(zhǔn)上，不能認為某藥物對治療組與對照組的療效不同。問題1

29、:該資料的分析方法是否合適？為什么？問題2 :應(yīng)該如何分析該資料？附：spss程序1行列表資料的x 2檢驗以實習(xí)三第5題為例。數(shù)據(jù)文件：“例3-5. sav”。數(shù)據(jù)格式：3列8行。2個分類變量，“年份”和“醫(yī)療質(zhì)量” ；1個頻數(shù)變量“廠”。 程序：weightby fcrosstabs/tables=年份by醫(yī)療質(zhì)量/format= avalue tables/statistic=chisq/cells= count row/count round cell .2. 兩組等級資料比較的秩和檢驗以實習(xí)三第5題為例。數(shù)據(jù)文件：“例3-5. sav”。數(shù)據(jù)格式：3列8行。1個分組變量“年份”，1個反

30、應(yīng)變量“醫(yī)療質(zhì)量”，1個頻數(shù)變量 “ f ”丄0程序：weightby f .npar tests/m-w=醫(yī)療質(zhì)量by年份(1 2)/missing analysis.附：sas程序1poisson分布的樣本均數(shù)與總體均數(shù)比較(直接法)以實習(xí)三第4題為例。data ex3 4;/*建立數(shù)據(jù)集*/n=l 00000;/ *確定樣本例數(shù)*/pai=0. 000007;/*確定總體率*/lam=n*pai;/*計算總體均數(shù)*/x=17;/*確定實際發(fā)生數(shù)*/p=l-poisson (lam,x-1);/*計算實際發(fā)生數(shù)所對應(yīng)的概率*/proc print;/*調(diào)用print過程*/var lam

31、p;/*顯示變量lam和p */run;2行x列表資料的x 2檢驗以實習(xí)三第5題為例。data ex3-5;/*建立數(shù)據(jù)集*/input r c f a)a);/*確定變量名稱，r為行變量，c為列變量, f為頻數(shù)變量*/cards;/*變量賦值*/略proc freq;/*調(diào)用freq過程*/weight f;/*定義f為頻數(shù)變量*/tables r*c/*作r*c的列聯(lián)表*/chisq;/*對列聯(lián)表x 2檢驗*/run;3.兩組等級資料比較的秩和檢驗以實習(xí)三第5題為例。data ex3 5;/*建立數(shù)據(jù)集*/input c g fw;/*確定變量名稱*/cards;/*變量賦值*/略proc

32、 nparlway wilcoxon;/*調(diào)用nparlway過程,進行wi lcoxon分析 /freq f;/*確定頻數(shù)變量為f*/var g;/*定義分析變量g*/class c;/*定義分組變量c*/run;實習(xí)四雙變量直線相關(guān)與回歸一、目的與要求：1、理解直線回歸的意義，掌握直線回歸方程的求法、回歸系數(shù)的意義及其假設(shè) 檢驗的方法（方差分析與t檢驗的基本思想）；2、掌握直線回歸的區(qū)間估計方法；3、掌握直線相關(guān)的概念、相關(guān)系數(shù)的意義，相關(guān)系數(shù)的計算方法、假設(shè)檢驗方 法（查表法與t檢驗法）；4、掌握秩相關(guān)分析適用的條件、spearman等級相關(guān)系數(shù)的意義，分析的步驟 和等級相關(guān)系數(shù)r s的

33、校正。二、學(xué)時：3學(xué)時。三、案例：1為了探討鎘對機體免疫功能的影響，分別對每組20只lace小鼠以劑量為0.3、1.2 .2. 4mg/kg/da卩的氯化鎘灌胃染毒14天，以20只未施染毒的小鼠為對照，分別測定小鼠脾淋巴細胞內(nèi)鈣調(diào)素含量（10 5ng/kg ）見表 4-1 :表4-1學(xué)毒計量與鈣調(diào)素含量的牛組信息染毒劑量（mg/kg/da卩）|樣本含量憐丐調(diào)素含量（10 5ng/kg j對照組（0 ）205. 38 ±2. 860. 32 04. 68 ±2. 721.2204. 32 ±2. 262.4203. 70 ±2. 67對這份資料，研究人員作

34、了如下的統(tǒng)計分析：四個劑量組間比較的方差分析、 鈣調(diào)素均值與染毒劑量的相關(guān)系數(shù)、鈣調(diào)素均值關(guān)于染毒劑量的線性回歸。結(jié)論 為：（1）經(jīng)方差分析得f=1.417,p>0. 05 ,認為為三個實驗組的鈣調(diào)素含量差異無統(tǒng) 計學(xué)意義；（2）必各2且的染毒劑量和鈣調(diào)素的組均值計算的相關(guān)系數(shù)（r=-0. 9996,p<0. 05 ）,故又結(jié)論為：染毒劑量與鈣調(diào)素含量呈負相關(guān)；（3）通過染毒劑量預(yù)測鈣調(diào)素含量效果很好。問題1 :對該研究數(shù)據(jù)進行方差分析的目的是什么？問題2 :染毒劑量和鈣調(diào)素的相關(guān)分析應(yīng)該怎么做？問題3 :為了探討小鼠脾淋巴細胞內(nèi)鈣調(diào)素含量與氯化鎘染毒劑量的劑量一 反應(yīng)關(guān)系，應(yīng)采用

35、何種統(tǒng)計分析方法？問題4 :研究人員上述做法存在何種問題？2 .某地10名一年級女大學(xué)生的胸圍（cm ）與肺活量（l ）數(shù)據(jù)見表4-2所 示。表4-2 10名一年級女大學(xué)生的胸圍（cm ）與肺活量（l ）學(xué)生編號12891034567胸圍/72. 583. 974. 873. 779. 478. 388.477. 181. 778. 3肺活 y2. 513. 111. 912. 983. 282. 723. 382. 832. 862. 72問題1 :按此資料繪制散點圖？問題2 :求直線回歸方程并對回歸系數(shù)作假設(shè)檢驗。問題3 :求直線相關(guān)系數(shù)，相關(guān)系數(shù)的假設(shè)檢驗結(jié)果與直線回歸系數(shù)會一致嗎？ 為

36、什么？問題4 :試估計胸圍為75cm時的平均肺活量，計算其95 %的可信區(qū)間，并 說明其含義。問題5 :求胸圍為75cm時，某地一年級大學(xué)生肺活量的95 %的散布范圍， 并解釋其含義。3 .對某省不同地區(qū)水質(zhì)的碘含量及其甲狀腺腫的患病率作了調(diào)查，結(jié)果見表4-3 :表4-3某省不同地區(qū)水質(zhì)碘含量與甲狀腺腫患病率地區(qū)碘含量 (mg/l )患病率（%）地區(qū)碘含量(pg/l )患病率 （%）11. 040. 5107. 76. 322. 037. 7118. 07. 132. 539. 0128. 09. 043. 520. 0138. 34. 053. 522. 0148. 54. 064. 037

37、.4158. 55.474. 431. 516& 84. 784. 515. 61724. 50. 094. 621. 0研究者發(fā)現(xiàn)不同地區(qū)的甲狀腺腫的患病率的高低與本地區(qū)水質(zhì)的碘含量有關(guān)，于 是利用pearson積差相關(guān)的計算公式，把碘含量視為變量把甲狀腺腫的患 病率視為因變量卩，計算出相關(guān)系數(shù)，得r =-0.712 ,經(jīng)檢驗p <0.002,據(jù) 此認為甲狀腺腫的患病率與水質(zhì)的碘含量之間有負相關(guān)關(guān)系。問題1 :這是否正確？為什么？問題2 :應(yīng)當(dāng)搜集哪些資料？如何分析？附：spss程序直線相關(guān)回歸分析(含散點圖)以實習(xí)四第2題為例。數(shù)據(jù)文件：“例4-2sav”。數(shù)據(jù)格式：2列10

38、行。1個自變量“無”，1個因變量“ y ”。程序：regression/missing listwise/statistics coeff outs ci r anova/criteria=pin(. 05) poutc 10)/noorigin/dependent y/method=enter x graph/scatterplot (bivar)=x with y/missing=listwise 附：sas程序1 兩個變量的直線回歸分析：以實習(xí)四第2題為例。data ex4 2;/*建立數(shù)據(jù)集*/input x y;/*確定變量名稱*/cards;/*變量賦值*/略proc reg;/*

39、調(diào)用reg過程*/model y=x;/*定義模型，以y為應(yīng)變量，以x為自變量 /run;在model語句后面可以加上選項，得到一些有用的統(tǒng)計量。常用選項包括: stb :輸出標(biāo)準(zhǔn)化偏回歸系數(shù)。p :輸出每個觀測的實際值、預(yù)測值和殘差。cli :輸出每個觀測預(yù)測值均數(shù)的雙側(cè)95%置信區(qū)間， elm :輸出每個觀測預(yù)測值的雙側(cè)95%置信范圍。2.兩個變量的直線相關(guān)分析：以實習(xí)四第2題為例。data ex4 2;/*建立數(shù)據(jù)集*/input x y;/*確定變量名稱*/cards;/*變量賦值*/略proc corr;"調(diào)用co"過程*/var x y/*確定作相關(guān)分析的變量*

40、/run;3兩個變量的秩相關(guān)分析：以實習(xí)四第3題為例。data ex4.3;/*建立數(shù)據(jù)集*/input x y;/*確定變量名稱*/cards;/*變量賦值*/略proc corr spearman;/*調(diào)用corr過程，要求作spearman相關(guān)分 析*/var x y;/*確定作相關(guān)分析的變量*/run;實習(xí)五多元線性回歸分析、logistic回歸分析一、目的與要求：1、掌握多元回歸分析的概念和偏回歸系數(shù)、復(fù)相關(guān)系數(shù)，校正復(fù)相關(guān)系數(shù)，剩 余標(biāo)準(zhǔn)差的意義；2、了解多元回歸分析和logistic回歸分析 的步驟；3、了解多元回歸分析和logistic回歸分析中自變量選擇的方法。二、學(xué)時：3學(xué)

41、時。三、案例：1有學(xué)者認為，血清中低密度脂蛋白增高和高密度脂蛋白降低是引起動脈硬化 的一個重要原因。現(xiàn)測量了 30名被懷疑患有動脈硬化的就診患者的載脂蛋白a i 、載脂蛋白b、載脂蛋白e、載脂蛋白c、低密度脂蛋白中的膽固醇、高 密度脂蛋白中的膽固醇含量，資料見表5-1 :表5-1 30名就診患者血清中低、高密度脂蛋白中的膽固醇含量及載脂蛋白的測量結(jié)果序號1載脂蛋白a i(mg/dl)x 1載脂蛋白b (mg/dl) xi載脂蛋白e (mg/dl)x 3載脂蛋白c (mg/dl)x 4低密度脂蛋白(mg/dl)y 1高密度脂蛋白(mg/dl)y 211731067. 014. 713762213

42、91326.417. 81624331981126. 916. 71348141181387. 115. 7188395139948. 613.613851617516012. 120. 321565713115411. 221. 51714081581419. 729. 61484291581377. 41& 219756101321517. 517. 2113371116211u().()15. 9145701214411310. 142. 88141131621377. 22 0. 718?5()141691298. 516. 715758151291386. 310. 1197

43、471616614811. 533.415649171851186. 017. 515669181551216. 120.415457191751114. 127. 214474201361109. 426. 09039211531338. 516. 921565221101499. 524. 71844023160865. 310. 81185724112123& 016. 61273425147110& 5413754262041226. 121. 012672271311026. 613.413 05128170127& 424. 71356229173123&a

44、mp; 719. 0188853013213113. 829. 212238用spss或sas統(tǒng)計軟件完成：問題1 :分別求孑1,孑2關(guān)于才1,才2,才3,才4的線性回歸方 程，并做分析。問題2 :分別用向前法、向后法和逐步回歸法選擇變量，看結(jié)果是否一致。 問題3 :作7 2/71關(guān)于zl , xi ,/3 ,*4的逐步回歸分析，并與前面的分析結(jié)果進行比較。問題4 :進一步做殘差分析，看是否滿足回歸分析的條件和有無離群值。問題5 :根據(jù)上面的分析，從專業(yè)上用文字扼要地敘述最后的分析結(jié)果。問題6 :有人做同類研究時，得到載體蛋白/i與高密度脂蛋白/ 2的回歸 系數(shù)為負值，試分析產(chǎn)生這一錯誤結(jié)論的

45、原因。2某醫(yī)院骨科在研究某種藥物治療骨折效果時，收集了 516例病例資料，對 每一患者采用相同的標(biāo)準(zhǔn)按照“好、差”做療效評價。在評價時需要同時考 慮骨折的類型、是否手術(shù)、是否服藥和治療時間4個因素對結(jié)果的影響，各因 素的說明及部分計算結(jié)果見表5-2 :表5-2各因素與說明因素變量名賦值說明骨折類型x 1閉合=0 ,開放=1治療方法j 2非手術(shù)=0 ,手術(shù)=1服藥情況x 3未服藥=0 ,服藥=1治療周數(shù)x 41 =1 ,11 =2 ,21 =3療效評價y差=0 ,好=1將治療周數(shù)化為啞變量水平z 4 -1才4 - 2100210301模型1 : 2p立尿辰屁£4屈磯夙elogistic

46、回歸參數(shù)估計（模型1）估計x 1x 2x 3x 4 1x 4 2回歸系數(shù)（勺）-0. 148- 0. 0062. 2202. 6544 . 745標(biāo)準(zhǔn)誤（鬲）0. 2190. 1890. 2570. 2430. 376變量標(biāo)準(zhǔn)差（%0. 7010. 6680. 7760. 7680. 666模型對數(shù)似然函數(shù)值：in l 1 =- 423.033模型 2 :a *a-ti + ajti* axi*+ ajvh.logistic回歸參數(shù)估計（模型2）估計x 1/ 2/ 3x 4.1x 4.2xxxi回歸系數(shù)（坊）- 1. 092-0. 3172. 3512. 7384. 9201. 430標(biāo)準(zhǔn)誤（

47、）0. 3840. 2150. 2630. 2470. 3840. 474模型對數(shù)似然函數(shù)值：in z 2 = -418. 292問題1 ：根據(jù)擬合模型1的結(jié)果說明，在相同的治療時間條件下，藥物及手術(shù) 的作用。問題2 :根據(jù)擬合模型1的結(jié)果說明，各因素對結(jié)果的影響大小順序及不同治 療時間的優(yōu)勢比值。問題3 :根據(jù)擬合模型2的結(jié)果說明，*1*2是否有必要考慮？手術(shù)究竟 有無作用？問題4 :根據(jù)擬合模型2的結(jié)果說明，在開放性骨折中，藥物和手術(shù)同時使用相對于單純手術(shù)作用有多大?3某醫(yī)生在河南平頂山煤礦區(qū)人群糖尿病（dm）現(xiàn)況調(diào)查基礎(chǔ)上，對篩選出的 174例糖尿病病例和3 066例糖耐量正常者進行以人

48、群為基礎(chǔ)的病例-對照 研究。調(diào)查內(nèi)容包括性別、年齡、糖尿病史、肥胖、體力勞動、飲酒和飲食等因 素，分析目的主要是糖尿病患病是否與肥胖有關(guān)。該醫(yī)生應(yīng)用非條件logistic 回歸分析糖尿病發(fā)生與各種危險因素的關(guān)聯(lián)性，先做糖尿病患病與逐個因素的單 因素logistic回歸結(jié)果摘要在表5-3中。然后將單因素分析有統(tǒng)計學(xué)意義的 危險因素引入多因素logistic回歸模型進行逐步篩選，篩選結(jié)果列在表5-4中。結(jié)果表明，年齡大、母親有糖尿病史、同胞有糖尿病史、最重時體質(zhì)指數(shù)(bmi) 高、腰臀比值(whr)高、舒張壓高、多食高梁和豆類可能是糖尿病患病的獨立 危險因子；職業(yè)性體力活動強度高和多食淺色蔬菜可能是糖尿病患病的獨立保護 因子。表5-3非條件logistic回歸單因素分析結(jié)果變量oror 的 95%ci變量oror 的 95%ci性別0. 520. 380. 71年齡2. 442. 08 2. 87居住年限1. 581. 24 2. 03母親糖尿病史5. 933. 25 10. 82子女有糖尿病8. 952. 2236. 08同胞糖尿病史4. 992. 44 10. 21現(xiàn)時bmi3. 352. 41 4. 65最重時bmi4. 443. 01 6. 56現(xiàn)時whr5. 914. 05 8. 63飲酒指數(shù)0. 670. 47 0. 96職業(yè)性體力活動0. 710. 630. x