SAS統(tǒng)計分析-T檢驗課件

1、t 檢驗鄧 偉2012.121t 檢驗鄧 偉1 t檢驗簡介均數(shù)的抽樣分布規(guī)律:來自正態(tài)總體或來自偏態(tài)總體但樣本量比較大的樣本其均數(shù)服從正態(tài)分布標準變換后u統(tǒng)計量:uN(0,1),以樣本標準差s代替總體標準差,此時變換后的統(tǒng)計量稱為t統(tǒng)計量, t統(tǒng)計量服從自由度為-1的t分布t檢驗的種類單樣本均數(shù)的比較：樣本所代表的總體均數(shù)和某已知總體均數(shù)比較配對設(shè)計兩樣本的比較：每對差值的總體均數(shù)是否為0完全隨機設(shè)計兩樣本：兩樣本所代表的總體均數(shù)是否相等2 t檢驗簡介均數(shù)的抽樣分布規(guī)律:2配對t檢驗原理:檢驗差值d的總體均數(shù)是否為0假設(shè): 公式:常見配對的情況同一對象治療（或處理）前后的比較同一樣本接受不同處

2、理的比較配對的兩個受試對象分別給予兩種處理3配對t檢驗原理:檢驗差值d的總體均數(shù)是否為03 data ttest1; input x1 x2; /*輸入治療前、后的血紅蛋白值*/ d=x1-x2; /*計算治療前后的差值 */ cards; 113 140 150 138 150 140 135 135 128 135 100 120 110 147 120 114 130 138 123 120 ; proc means mean stderr t prt; /*輸出對均數(shù)是否為0的t檢驗及其p值*/ var d; run;配對t檢驗4 data ttest1;配對t檢驗4 Analysis

3、 Variable : D Mean Std Error T Prob|T| - -6.8000000 5.2042717 -1.3066189 0.2237 -配對t檢驗5 Analysis Variable : D配對t檢驗5成組t檢驗 原理:兩樣本所代表的總體均數(shù)是否相等,檢驗兩個樣本均數(shù)的差是否為0假設(shè):公式6成組t檢驗 原理:兩樣本所代表的總體均數(shù)是否相等,檢驗兩個樣成組t檢驗 應(yīng)用條件獨立性正態(tài)性univariate 正態(tài)性檢驗方差齊性齊性檢驗t檢驗7成組t檢驗 應(yīng)用條件7data ttest3; do c=1 to 2; /* c循環(huán)開始,循環(huán)次數(shù)為2 */ input n; /

4、* 輸入n值 */ do i=1 to n; /* i循環(huán)開始,循環(huán)次數(shù)為n */ input x ; /* 輸入x 值*/ output; /* 將x輸出 */ end; /* i循環(huán)結(jié)束 */ end; /* c循環(huán)結(jié)束 */ cards; 5 279 334 303 338 198 3 229 274 310;proc print;proc ttest;class c; var x;run;成組t檢驗 8data ttest3;cards;成組t檢驗 8循環(huán)語句的應(yīng)用 do ; /* 循環(huán)1開始*/ do ; /* 循環(huán)2開始*/ input ; output; end; /* 循環(huán)2結(jié)

5、束 */ end; /* 循環(huán)1結(jié)束 */注意：1）的合理應(yīng)用 2）output的作用：9循環(huán)語句的應(yīng)用 do ; /* 循環(huán)1開始/* 循環(huán)的結(jié)果 * / OBS C N I X 1 1 5 1 279 2 1 5 2 334 3 1 5 3 303 4 1 5 4 338 5 1 5 5 198 6 2 3 1 229 7 2 3 2 274 8 2 3 3 310成組t檢驗 10/* 循環(huán)的結(jié)果 * /成組t檢驗 10 TTEST PROCEDUREVariable: XC N Mean Std Dev Std Error-1 5 290.40000000 56.99385932 25.

6、488428752 3 271.00000000 40.58324778 23.43074903Variances T DF Prob|T|-Unequal 0.5603 5.6 0.5971Equal 0.5099 6.0 0.6284For H0: Variances are equal, F = 1.97 DF = (4,2) ProbF = 0.7272成組t檢驗 11 TTEST PROCEDURE成組t檢驗 11單樣本t test原理:檢驗?zāi)硺颖舅鶃碜缘目傮w均數(shù)與已知總體均數(shù)是否相等假設(shè): 公式:例:某疾病的平均住院天數(shù)是否是80天？12單樣本t test原理:檢驗?zāi)硺颖舅鶃碜缘目?/p>

7、體均數(shù)與已知總體/*程序1*/data time;input time ;cards;43 90 84 87 116 95 86 99 93 92121 71 66 98 79 102 60 112 105 98;proc ttest h0=80;var time;run;/* 程序2 */Data time;Set time;D=time-80;Proc means stderr t prt;Run;單樣本t test13/*程序1*/proc ttest h0=80;單樣本t tSAS程序 :proc ttestProc ttest options; /* 過程名 */Class ; /*

8、 分組變量名,決定哪兩組進行比較 */Var ; /* 反應(yīng)變量列表 */PAIRED var1*var2; /*專用于配對變量檢驗, class和var語句在使用該語句時無效*/By variables; /*指明分層變量,按該變量分別分析,需先排序 */Freq variable; /* 指明頻數(shù)變量 */Run;14SAS程序 :proc ttestProc ttest optionsAlpha 定義可信區(qū)間的alpha水平，默認為0.05H0 定義特定的總體均數(shù)值，默認為0by：根據(jù)by變量分別計算分析結(jié)果，要求事先根據(jù)by變量排序class：定義分組變量，決定進行哪兩組比較。在tte

9、st中只允許有兩組SAS程序 :proc ttest15optionsSAS程序 :proc ttest15Freqweight非常接近Freq 頻數(shù)，只能是整數(shù)Weight可以是小數(shù)SAS程序 :proc ttest16FreqweightSAS程序 :proc ttest16 Paired定義配對變量SAS程序 :proc ttest17 Paired定義配對變量SAS程序 :proc tte例1的分析用ttest過程解決: Proc ttest data=ttest1; Paired x1*x2; Run; Lower CL Upper CL Lower CL Upper CLDiffe

10、rence N Mean Mean Mean Std Dev Std Dev Std Dev Std Err x1 - x2 10 -18.57 -6.8 4.9729 11.32 16.457 30.045 5.2043 T-Tests Difference DF t Value Pr |t| x1 - x2 9 -1.31 0.223718例1的分析用ttest過程解決:18SAS的循環(huán)語句介紹記數(shù)循環(huán)語法 do 記數(shù)變量起始值 to 終止值 by 增減值 .end；條件循環(huán)語句 Do while （條件）；End； Do until （條件）；End； 19SAS的循環(huán)語句介紹記數(shù)循環(huán)語法 19記數(shù)循環(huán)語法data a;amount=1000;rate=0.08;year_int=amount*rate;do i=1 to 12; mont_int+(mont_int+amount)*rate/12;end;run;proc print;run; 20記數(shù)循環(huán)