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

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

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

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

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

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

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

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

8、驗(yàn),例1：已知正常人乙酰膽堿脂酶的平均值為1.44單位，現(xiàn)測得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。問：慢性支氣管炎患者與正常人的平均乙酰膽堿脂酶之間的差別有無顯著性意義？,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過程,TTEST過程,proc ttest data=aa h0=1.44; /*設(shè)立無效假設(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)過輻射處理后,得到一穩(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)定品系的百粒重和輻射前有無差別?,配對t檢驗(yàn),例3

11、：擬對某種藥物的降壓效果進(jìn)行評價，治療前后病人的舒張壓水平如下：問：治療前后舒張壓的差異有無統(tǒng)計學(xué)意義？,UNIVARIATE、MEANS和TTEST過程,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;,配對t檢驗(yàn)-練習(xí),例4：將同種屬的16只大白鼠按性別相同，年齡、體重接近者配成8對，將同對中的兩只大白鼠隨機(jī)分配到正常飼料組和維生素E缺乏飼料組，以兩種飼料喂養(yǎng)一段時間后，測定其肝中維生素A的含量（IU/g）,測得數(shù)據(jù)如下：試比較不同飼料組大白鼠肝中維生素A的含量,兩樣本t檢驗(yàn),例5：測得 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問克山病患者的血磷值與健康人是否不同?,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過程：,兩樣本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問：兩工種工人尿氟含量有無差別？,綜合練習(xí)（1）,將lbw.dbf數(shù)據(jù)庫轉(zhuǎn)為SAS數(shù)據(jù)集；數(shù)據(jù)中哪些變量為定量變量，哪些為定性變量？,綜合練習(xí)（2）,既往有吸煙史的母親和無吸煙史的母親，