單因素分析一-定量資料-2013319

1、流行病學(xué)系 廖巍,單因素分析-定量資料,第一部分：背景知識(shí)復(fù)習(xí),變量（variable）和資料（data）資料的類(lèi)型定量資料：對(duì)觀察對(duì)象進(jìn)行數(shù)值上的測(cè)定，一般有計(jì)量單位，可直接比較大小。定性資料：表示觀察對(duì)象的類(lèi)別或?qū)傩?，一般無(wú)計(jì)量單位。二分類(lèi)資料 、無(wú)序多分類(lèi)資料、等級(jí)資料（有序資料）類(lèi)型轉(zhuǎn)化：定量資料?定性資料,流行病學(xué)資料分析：探究疾病結(jié)局（因變量）與病因/影響因素（自變量）之間的關(guān)系。因變量和自變量因變

2、量(Dependent variable)：反應(yīng)變量、結(jié)局變量自變量 (Independent variable): 效應(yīng)變量、處理因素根據(jù)研究假設(shè)中自變量的數(shù)目，可分為：?jiǎn)我蛩胤治龆嘁蛩胤治龈鶕?jù)研究變量的資料類(lèi)型，可分為：定量資料分析定性資料分析,單因素分析—母親是否吸煙對(duì)嬰兒出生體重的影響,定量資料,定性資料,,,定量資料的統(tǒng)計(jì)分析：統(tǒng)計(jì)描述：了解樣本觀察值的分布特征集中趨勢(shì)離散趨勢(shì)統(tǒng)計(jì)推斷：根據(jù)樣本提供的信

3、息和抽樣分布規(guī)律，以一定的概率推斷總體特征參數(shù)估計(jì)：點(diǎn)估計(jì)、區(qū)間估計(jì)假設(shè)檢驗(yàn)：參數(shù)檢驗(yàn)：t 檢驗(yàn)、F 檢驗(yàn)非參數(shù)檢驗(yàn)：秩和檢驗(yàn),t檢驗(yàn),步驟：建立檢驗(yàn)假設(shè)（H0；H1），確定檢驗(yàn)水準(zhǔn)（α=0.05）計(jì)算檢驗(yàn)統(tǒng)計(jì)量t確定p值，作出推斷結(jié)論按α=0.05 水準(zhǔn)，拒絕……，接受……，有/無(wú)統(tǒng)計(jì)學(xué)意義，可以認(rèn)為……類(lèi)型：?jiǎn)螛颖総檢驗(yàn)配對(duì)t檢驗(yàn)兩獨(dú)立樣本t檢驗(yàn),適用：樣本均數(shù) 與已知總體均數(shù) μ0的比較舉例：擬研究經(jīng)

4、常參加體育鍛煉的中學(xué)生心臟功能是否與一般的中學(xué)生相等。研究方法：已知中學(xué)一般男生的心律平均值為74次/分鐘；隨機(jī)抽取某地區(qū)中學(xué)常年參加體育鍛煉的男生16名，測(cè)量心律；比較該16名男生的心率平均值是否與中學(xué)一般男生的心律平均值相等。,單樣本t檢驗(yàn),配對(duì)t檢驗(yàn),適用：配對(duì)設(shè)計(jì)計(jì)量資料均數(shù)的比較，包括：同一受試對(duì)象處理前后；同一受試對(duì)象接受兩種不同的處理；兩個(gè)同質(zhì)受試對(duì)象分別接受兩種不同的處理。舉例：擬研究某種藥物的降壓效果

5、。研究方法：確定研究對(duì)象；研究開(kāi)始的時(shí)候，測(cè)量并記錄血壓水平；服用某種降壓藥一年；研究結(jié)束的時(shí)候，測(cè)量并記錄血壓水平；評(píng)價(jià)該藥物的降壓效果。,兩獨(dú)立樣本t檢驗(yàn),適用：完全隨機(jī)設(shè)計(jì)兩獨(dú)立樣本的均數(shù)比較舉例：擬研究某種藥物的降壓效果。研究方法：確定研究對(duì)象；隨機(jī)分組：藥物組：服用該降壓藥；對(duì)照組：不服用該藥物；隨訪(fǎng)一段時(shí)間；比較兩組血壓下降值是否相同。,t檢驗(yàn)小結(jié),第二部分：t檢驗(yàn)主要SAS程序介紹,可用于t檢

6、驗(yàn)的三個(gè)SAS過(guò)程步,UNIVARIATE、MEANS、TTEST必需結(jié)構(gòu)：相同作用：計(jì)算基本統(tǒng)計(jì)量，均值、標(biāo)準(zhǔn)差等。,UNIVARIATE 過(guò)程步,應(yīng)用：?jiǎn)螛颖総檢驗(yàn)配對(duì)t檢驗(yàn)正態(tài)性檢驗(yàn)proc univariate normal；此步驟進(jìn)行正態(tài)性檢驗(yàn)，要有normal,MEANS 過(guò)程步，做均數(shù)處理,應(yīng)用：?jiǎn)螛颖総檢驗(yàn)配對(duì)t檢驗(yàn)proc means ;常用可選擇統(tǒng)計(jì)量：,n 樣本量me

7、an 均數(shù)std 標(biāo)準(zhǔn)差stderr 標(biāo)準(zhǔn)誤,t t檢驗(yàn)的統(tǒng)計(jì)量prt 概率值 (p)lclm 可信區(qū)間下限uclm 可信區(qū)間上限,TTEST 過(guò)程步,應(yīng)用：?jiǎn)螛颖総檢驗(yàn)配對(duì)t檢驗(yàn)兩獨(dú)立樣本t檢驗(yàn)方差齊性檢驗(yàn)proc ttest ；,第三部分：SAS分析實(shí)例,學(xué)習(xí)目標(biāo),掌握SAS數(shù)據(jù)集的建立掌握常用的SAS過(guò)程步正確解讀SAS運(yùn)行結(jié)果,單樣本t檢

8、驗(yàn),例1：已知正常人乙酰膽堿脂酶的平均值為1.44單位，現(xiàn)測(cè)得13例慢性氣管炎患者的乙酰膽堿脂酶分別為：1.50，2.19，2.32，2.41，2.11，2.54，2.20，2.36，1.42，2.17，1.84，1.96，2.39。問(wèn)：慢性支氣管炎患者與正常人的平均乙酰膽堿脂酶之間的差別有無(wú)顯著性意義？,data aa;input x @@;d=x-1.44;新加入的變量要放在cards前面cards;1.50 2.19

9、 2.32 2.41 2.11 2.54 2.20 2.36 1.42 2.17 1.84 1.96 2.39;run;,UNIVARIATE、MEANS過(guò)程,TTEST過(guò)程,proc ttest data=aa h0=1.44; /*設(shè)立無(wú)效假設(shè)H0的數(shù)值 (0為零，不是英文字母o)*/var x;run;,data aa;input x @@;cards;1.50 2.19 2.32 2.41 2.11 2.54 2.

10、20 2.36 1.42 2.17 1.84 1.96 2.39;run;,單樣本t檢驗(yàn)-練習(xí),例2: 某大豆品種百粒重15g,經(jīng)過(guò)輻射處理后,得到一穩(wěn)定品系,隨機(jī)抽樣15次,百粒重(g)分別為：18.7, 19.3, 19.2, 18.4, 18.3, 18.6, 18.7, 18.3, 18.5, 18.9, l9.1, 17.7, 17.8, 17.9, 18.2。問(wèn)該穩(wěn)定品系的百粒重和輻射前有無(wú)差別?,配對(duì)t檢驗(yàn),例3

11、：擬對(duì)某種藥物的降壓效果進(jìn)行評(píng)價(jià)，治療前后病人的舒張壓水平如下：?jiǎn)枺褐委熐昂笫鎻垑旱牟町愑袩o(wú)統(tǒng)計(jì)學(xué)意義？,UNIVARIATE、MEANS和TTEST過(guò)程,data cc;input before after @@;d=before-after;cards;117 123 127 108 141 120 107 107 110 100 114 98 115 102 138 152 127 104 122 107;r

12、un;,proc ttest data=cc; var d; run;,配對(duì)t檢驗(yàn)-練習(xí),例4：將同種屬的16只大白鼠按性別相同，年齡、體重接近者配成8對(duì)，將同對(duì)中的兩只大白鼠隨機(jī)分配到正常飼料組和維生素E缺乏飼料組，以?xún)煞N飼料喂養(yǎng)一段時(shí)間后，測(cè)定其肝中維生素A的含量（IU/g）,測(cè)得數(shù)據(jù)如下：試比較不同飼料組大白鼠肝中維生素A的含量,兩樣本t檢驗(yàn),例5：測(cè)得 11 例克山病患者與 13 名健康人的血

13、磷值(mmol/L)如下： 患者：　0.84　1.05　1.20　1.20　1.39　1.53　1.67　1.80　1.87　2.07　2.11 健康者：　0.54　0.64　0.64　0.75　0.76　0.81　1.16　1.20　1.34　1.35　1.48 1.56 1.87問(wèn)克山病患者的血磷值與健康人是否不同?,data ee;input x@@;cards;0.84 1.05 1.20 1.20 1.

14、39 1.53 1.67 1.80 1.87 2.07 2.110.54 0.64 0.64 0.75 0.76 0.81 1.16 1.20 1.34 1.35 1.48 1.56 1.87;run;,data ee;input g x@@;cards;1 0.84 1 1.05 1 1.20 1 1.20 1 1.39 1 1.53 1 1.67 1 1.80 1 1.87 1 2.07 1 2.112 0.54 2

15、0.64 2 0.64 2 0.75 2 0.76 2 0.81 2 1.16 2 1.20 2 1.34 2 1.35 2 1.48 2 1.56 2 1.87;run;,data ee;input x@@;if _n_<12 then g=1;else g=2;cards;0.84 1.05 1.20 1.20 1.39 1.53 1.67 1.80 1.87 2.07 2.110.54 0.64 0.64

16、0.75 0.76 0.81 1.16 1.20 1.34 1.35 1.48 1.56 1.87;run;,X,兩樣本t檢驗(yàn)前提條件：,獨(dú)立性正態(tài)分布方差齊性,TTEST過(guò)程：,兩樣本t檢驗(yàn)-練習(xí),例6：檢查某工廠21名氟作業(yè)工人和15名非氟作業(yè)工人的尿氟含量（mg）如下：氟作業(yè)工人：3.40 3.00 2.80 2.16 2.18 2.20 2.40 2.30 2.80 1.60 1.52 1.74 1.14

17、 1.24 1.34 1.08 1.22 1.00 0.52 0.40 0.60非氟作業(yè)工人：1.76 1.44 1.08 0.15 0.40 0.52 0.42 0.62 0.46 0.40 0.72 0.84 0.20 0.85 0.48問(wèn)：兩工種工人尿氟含量有無(wú)差別？,綜合練習(xí)（1）,將lbw.dbf數(shù)據(jù)庫(kù)轉(zhuǎn)為SAS數(shù)據(jù)集；數(shù)據(jù)中哪些變量為定量變量，哪些為定性變量？,綜合練習(xí)（2）,既往有吸煙史的母親和無(wú)吸煙史的母親，