t檢驗(yàn)與u檢驗(yàn)

1、數(shù)據(jù)處理技術(shù)數(shù)據(jù)處理技術(shù) 2015-2016-12015-2016-1三三 t t檢驗(yàn)檢驗(yàn) u u檢驗(yàn)檢驗(yàn) 秩和檢驗(yàn)秩和檢驗(yàn)東南大學(xué)公共衛(wèi)生學(xué)院東南大學(xué)公共衛(wèi)生學(xué)院 孫金芳孫金芳 Email: Passwords: sas2015review a POWERFUL learning tool and makes it easier to find relevant guidance： /sas/v8/sashtml/main.htm The MEANS Procedure The UNIVARIATE Procedure The CHART Proce

2、dure The PLOT Procedure The GCHART Procedure The GPLOT ProcedureBase SASSAS/GRAPH本節(jié)主要內(nèi)容 t檢驗(yàn)(t test)、u(z)檢驗(yàn)、秩和檢驗(yàn)方法回顧與適用資料特點(diǎn) UNIVARIATE 過程與t 檢驗(yàn) TTest 過程 NPAR1WAY 過程t檢驗(yàn)與u檢驗(yàn) 正態(tài)性(Normality)：正態(tài)性檢驗(yàn) 方差相等(Homogenity, Equal variance)：方差齊性檢驗(yàn) 獨(dú)立性(Indenpendenct)：如果不滿足前兩個條件，非參數(shù)檢驗(yàn)非參數(shù)檢驗(yàn)。若兩樣本比較，僅方差不齊，可考慮t檢驗(yàn)大樣本時可選u(Z

3、)檢驗(yàn)近似，但軟件不能直接求得，需寫程序。Univariate 過程與t 檢驗(yàn) 單樣本形式(one sample)資料檢驗(yàn)： 已知樣本均數(shù)與總體均數(shù)的比較 配對資料的t檢驗(yàn)例3-1正常成年男子的血液中血小板計(jì)數(shù)的均數(shù)為225109/L，今測得20名男性油漆作業(yè)工人的血小板計(jì)數(shù)值（單位109/L）：220，188，162，230，145，160，238，188，247，113 126，245，164，231，256，183，190，158，224，175。問油漆作業(yè)工人的血小板計(jì)數(shù)與正常成年男子有無顯著差異？ 程序： DATA EX3_1; INPUT X; CARDS; 220 188 162

4、 230 145 160 238 188 247 113 126 245 164 231 256 183 190 158 224 175 ; PROC UNIVARIATE NORMAL mu0=225; VAR x; RUN;配對資料的t檢驗(yàn) 程序Data ex3_2;Input x1 x2;Y=x1-x2;/*賦值語句*/Cards;113 138 120 116 138 125 120 136 100 110 118 132 138 130 123 110;Proc univariate normal;Var y;Run;TTest過程：所有t檢驗(yàn)資料格式：*proc ttest opt

5、ions ( H0=mu/*單樣本t-test*/) ;class variable;/*兩樣本t檢驗(yàn)不可少*/Paired variables;/*配對t檢驗(yàn)*/ var variables;By variables; freq variable;Weight variable;run;Data exp4_2; Input x; Cards; 504 498 496 487 509 476 482 510 469 472 ; proc ttest h0=500;/*指定總體均數(shù)*/ run; The TTEST Procedure Statistics Lower CL Upper CL L

6、ower CL Upper CL Variable N Mean Mean Mean Std Dev Std Dev Std Dev Std Err x 10 479.4 490.3 501.2 10.482 15.239 27.821 4.8191 T-Tests Variable DF t Value Pr |t| x 9 -2.01 0.07502 配對設(shè)計(jì)資料的t檢驗(yàn)Data a;input x1 x2;Cards;113 138 120 116 138 125 120 136 100 110118 132 138 130 123 110;proc ttest;Paired x1*x2

7、;run;另:Paired x1*x2Paired x1*x2 a*b x1-x2 a-bPaired (x1-x2)*(b1-b2) x1-b1 x1-b2 x2-b1 x2-b2Paired (a b)(b c) a-b a-c b-cThe TTEST Procedure Statistics Lower CL Upper CL Lower CL Upper CL Difference N Mean Mean Mean Std Dev Std Dev Std Dev Std Err x1 - x2 8 -15.63 -3.375 8.8789 9.6911 14.657 29.832 5

8、.1822 T-Tests Difference DF t Value Pr |t| x1 - x2 7 -0.65 0.53573 成組設(shè)計(jì)(completely random design)兩樣本均數(shù)比較的t檢驗(yàn)例3-3 為研究國產(chǎn)四類新藥阿卡波糖膠囊的降血糖效果，某醫(yī)院用40名II型糖尿病病人進(jìn)行同期隨機(jī)對照試驗(yàn)。試驗(yàn)者將這些病人隨機(jī)等分到試驗(yàn)組（阿卡波糖膠囊組）和對照組（拜唐蘋膠囊組），分別測得試驗(yàn)開始前和8周后的空腹血糖，算得空腹血糖下降值如表3-2，能否認(rèn)為該國產(chǎn)四類新藥阿卡波糖膠囊與拜唐蘋膠囊對空腹血糖的降糖效果不同？DATA EX3_3;DO C=1 TO 2;DO I=1 T

9、O 20;INPUT X;OUTPUT;END;END;CARDS;-0.70 -5.60 2.00 2.80 0.70 3.50 4.00 5.80 7.10 -0.50 2.50 -1.60 1.70 3.00 0.40 4.50 4.60 2.50 6.00 -1.40 3.70 6.50 5.00 5.20 0.80 0.20 0.60 3.40 6.60 -1.10 6.00 3.80 2.00 1.60 2.00 2.20 1.20 3.10 1.70 -2.00 ; PROC UNIVARIATE NORMAL;VAR X;BY C;PROC TTEST COCHRAN;VAR

10、X;CLASS C;RUN; 循環(huán)語句“DO-END”對變量進(jìn)行賦值，并通過輸出語句OUTPUT將所賦的值輸出到數(shù)據(jù)集，以方便數(shù)據(jù)集的建立 . 此處為一個雙循環(huán)嵌套語句，由一個外循環(huán)嵌套一個內(nèi)循環(huán)。 從“DO I=1 TO 20”到第一個END為內(nèi)循環(huán)； 從“DO C=1 TO 2”到第二個END為外循環(huán)。P358數(shù)據(jù)21-4的錄入形式 Data ex21_4; /*不是ex21-4*/ Do a=1 to 2; Do b=1 to 2; do id= 1 to 6; Input freq;output; End;end;end; Cards; 2.1 2.2 2.4 2.1 2.3 1.9

11、1.3 1.0 0.9 1.0 1.2 1.1 0.8 1.0 1.1 0.9 1.2 0.9 0.6 0.4 0.7 0.6 0.8 0.5 ;The TTEST ProcedureStatistics Lower CL Upper CL Lower CL Upper CLVariable C N Mean Mean Mean Std Dev Std Dev Std Dev Std Err X 1 20 0.6328 2.065 3.4972 2.3272 3.0601 4.4696 0.6843 X 2 20 1.4922 2.625 3.7578 1.8408 2.4205 3.5353

12、 0.5412 X Diff (1-2) -2.326 -0.56 1.2062 2.2547 2.7589 3.5556 0.8724(第四部分)T-TestsVariable Method Variances DF t Value Pr |t|X Pooled Equal 38 -0.64 0.5248X Satterthwaite Unequal 36.1 -0.64 0.5250X Cochran Unequal 19 -0.64 0.5286（第五部分）Equality of VariancesVariable Method Num DF Den DF F Value Pr F X

13、Folded F 19 19 1.60 0.3153沒有原始觀測值的樣本均數(shù)與總體均數(shù)比較的t檢驗(yàn) 例3-4 某藥廠制劑車間用自動裝瓶機(jī)封裝藥液，在裝瓶機(jī)工作正常時，每瓶藥液凈重500克。某日隨機(jī)抽取10瓶成品，得凈重量的均數(shù)為490.3克，標(biāo)準(zhǔn)差為15.2392，問此時的裝瓶機(jī)工作是否正常? DATA EX3_4; N=10; MEAN=490.3; SD=15.2392; MU=500; T=(MEAN-MU)/(SD/SQRT(N); DF=N-1; P=(1-PROBT(ABS(T),DF)*2; PROC PRINT; VAR T P; RUN;U(Z)檢驗(yàn) 例3-5 青海省20名男

14、子153人，體重均數(shù)為57.38kg。標(biāo)準(zhǔn)差為5.78kg。內(nèi)蒙古自治區(qū)同年齡的男子206人，體重均數(shù)為55.92kg，標(biāo)準(zhǔn)差為5.17kg，問內(nèi)蒙古與青海省20歲成年男子的體重有無顯著差異？ DATA EX3_5; INPUT N1 N2 MEAN1 MEAN2 SD1 SD2; U=(MEAN1-MEAN2)/SQRT(SD1*2/N1+SD2*2/N2); P=(1-PROBNORM(ABS(U)*2; CARDS; 153 206 57.38 55.95 5.78 5.17 ; PROC PRINT; VAR U P; RUN; 導(dǎo)入excel數(shù)據(jù)的t檢驗(yàn) 臨床新藥數(shù)據(jù) 入組年齡基線水

15、平是否一致？SAS Online SAS/STAT - SAS/STAT Users Guide The TTEST Procedure Syntax Example秩和檢驗(yàn) 非正態(tài)性資料 等級資料 開口資料 方差不齊 NPAR1WAY過程 FREQ過程主要內(nèi)容 配對資料秩和檢驗(yàn)與UNIVARIATE過程 兩個獨(dú)立樣本秩和檢驗(yàn)與NPAR1WAY過程 多個獨(dú)立樣本秩和檢驗(yàn)與NPAR1WAY過程 配伍組設(shè)計(jì)資料秩和檢驗(yàn)與FREQ過程配對資料秩和檢驗(yàn)與UNIVARIATE過程 符號檢驗(yàn)與Wilcoxon符號秩和檢驗(yàn) UNIVARIATE 過程 Proc univariate normal; 例5-1

16、 利用兩種不同的方法測得10名學(xué)習(xí)困難兒童的智商(IQ) 見表5-1，試分析兩種方法的測量結(jié)果有無差異？ 編號12345678910方法A61605663566359566261方法B55544759516753656554Data ex;Input x y;D=x-y;Cards;61 55 60 54 56 47 63 59 56 51 63 67 59 53 56 65 62 65 61 54;Proc univariate normal;Var d;Run;結(jié)果一 先看正態(tài)性檢驗(yàn)結(jié)果 tests for normality 正態(tài)性檢驗(yàn) 檢驗(yàn) -統(tǒng)計(jì)量- -P 值- Shapiro-Wi

17、lk W 0.836773 Pr D 0.0197 Cramer-von Mises W-Sq 0.16066 Pr W-Sq 0.0147 Anderson-Darling A-Sq 0.822584 Pr A-Sq 0.0222結(jié)果二 再看位置檢驗(yàn)結(jié)果：tests for location:Mu=0 位置檢驗(yàn): Mu0=0 檢驗(yàn) -統(tǒng)計(jì)量- -P 值- 學(xué)生 t t 1.4499 Pr |t| 0.1810 符號 M 2 Pr = |M| 0.3438 符號秩 S 14.5 Pr = |S| 0.1543第二節(jié) 兩樣本秩和檢驗(yàn)NPAR1WAY過程 *PROC NPAR1WAY option

18、s; *CLASS 變量名; Freq 變量名; Var 變量名； Exact 統(tǒng)計(jì)量; Run;NPAR1WAY options; Data=sasdata; WILCOXON|MEDIAN|ANOVA WILCOXON檢驗(yàn)|中位數(shù)檢驗(yàn)|方差分析其余語句說明 Class 指定分組變量，必須語句 Freq 指定頻數(shù)變量，可選語句 Var 指定分析變量 Exact 指定精確檢驗(yàn)統(tǒng)計(jì)量例5-2 兩組工人的血鉛值比較非鉛作業(yè)組242629344358637287101鉛作業(yè)組828797121164208213Data ex;G=1;If _n_=11 then g=2;Input x;Cards;

19、24 26 2934 43 58 63 72 87 101 82 87 97 121 164 208 213;Proc npar1way wilcoxon;Class g;var x;Run;等級資料 例5-3 某醫(yī)院用某藥治療兩型慢性支氣管炎患者共216例，療效見表5-3。問該藥對兩型慢性支氣管炎的治療效果是否不同？ DATA EX5_3; DO J=1 TO 2; DO G=1 TO 4; INPUT N; OUTPUT; END; END; CARDS; 62411411 20371615 ; PROC NPAR1WAY Wilcoxon; VAR G; CLASS J; freq n; RUN;第三節(jié) 多個樣本秩和檢驗(yàn) Kruskal-Wallis H檢驗(yàn) 兩兩比較 : 檢驗(yàn)水準(zhǔn)需作調(diào)整：=/m，m為需比較的次數(shù) 例5-4 以小白鼠研究正常肝核糖核酸(RNA)對癌細(xì)胞的生物學(xué)作用，試驗(yàn)分為對照組(生理鹽水)、水層RNA組和酚層RNA組，分別用此三種不同處理誘導(dǎo)肝癌細(xì)胞的果糖二磷酸酯酶(FDP酶)活力，得數(shù)據(jù)見表5-4。問這三種不同處理的誘導(dǎo)作用是否不同？ DATA EX5_4;DO i=1 TO 8; 以列為樣本組輸入數(shù)據(jù)DO j=1 TO 3;INPUT X; OUTPUT;END;END;CARDS;2.79 3.83 5.41