t检验与u检验
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
/sas/v8/sashtml/main.htm
– – – – The MEANS Procedure The UNIVARIATE Procedure The CHART Procedure The PLOT Procedure
Base SAS
– The GCHART Procedure – The GPLOT Procedure
Std Err 0.6843 0.5412 0.8724
没有原始观测值的样本均数与总体均数 比较的t检验
• 例3-4 某药厂制剂车间用自动装瓶机封装药 液，在装瓶机工作正常时，每瓶药液净重 500克。某日随机抽取10瓶成品，得净重量 的均数为490.3克，标准差为15.2392，问 此时的装瓶机工作是否正常?
The TTEST Procedure Statistics Lower CL Upper CL Lower CL Upper CL Variable C N Mean Mean Mean Std Dev Std Dev Std Dev X 1 20 0.6328 2.065 3.4972 2.3272 3.0601 4.4696 X 2 20 1.4922 2.625 3.7578 1.8408 2.4205 3.5353 X Diff (1-2) -2.326 -0.56 1.2062 2.2547 2.7589 3.5556 (第四部分) T-Tests Variable Method Variances DF t Value Pr > |t| X Pooled Equal 38 -0.64 0.5248 X Satterthwaite Unequal 36.1 -0.64 0.5250 X Cochran Unequal 19 -0.64 0.5286 （第五部分） Equality of Variances Variable Method Num DF Den DF F Value Pr > F X Folded F 19 19 1.60 0.315验(t test)、u(z)检验、秩和检验方法 回顾与适用资料特点 • UNIVARIATE 过程与t 检验 • TTest 过程
• NPAR1WAY 过程
t检验与u检验
• 正态性(Normality)：正态性检验 • 方差相等(Homogenity, Equal variance)：方差齐性 检验 • 独立性(Indenpendenct)： 如果不满足前两个条件，非参数检验。若两样本比 较，仅方差不齐，可考虑t’检验 大样本时可选u(Z)检验近似，但软件不能直接求得， 需写程序。
表 3-2 试验组和对照组空腹血糖下降值(mmol/L) 试验组 (n1=20) 对照组 (n2=20) -0.70 2.50 3.70 6.00 -5.60 -1.60 6.50 3.80 2.00 1.70 5.00 2.00 2.80 3.00 5.20 1.60 0.70 0.40 0.80 2.00 3.50 4.50 0.20 2.20 4.00 4.60 0.60 1.20 5.80 2.50 3.40 3.10 7.10 6.00 6.60 1.70 -0.50 -1.40 -1.10 -2.00
• • • • • • • • • • •
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;
P358数据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 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 ;
配对资料的t检验
表 3-1 铁剂和饮食两种方法治疗后患者血红蛋白值（g/L） 患者对子号 1 2 3 4 5 6 7 8 铁剂治疗组 113 120 138 120 100 118 138 123 饮食治疗组 138 116 125 136 110 132 130 110
• 程序 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;
x1-x2 a-b x1-b1 x1-b2 x2-b1 x2-b2 a-b a-c b-c
The TTEST Procedure
Statistics Lower CL Difference x1 - x2 N 8 Mean Upper CL Lower CL Mean T-Tests Difference x1 - x2 DF 7 t Value -0.65 Pr > |t| 0.5357 Mean 8.8789 9.6911 Upper CL Std Dev Std Err 29.832 5.1822 14.657 Std Dev Std Dev
数据处理技术
2015-2016-1
三 t检验 u检验 秩和检验
东南大学公共卫生学院 孙金芳 Email: sas2015seu@ Passwords: sas2015
review
• a POWERFUL learning tool and makes it easier to find relevant guidance：
U(Z)检验
• 例3-5 青海省20名男子153人，体重均数为 57.38kg。标准差为5.78kg。内蒙古自治区 同年龄的男子206人，体重均数为55.92kg， 标准差为5.17kg，问内蒙古与青海省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;
Univariate 过程与t 检验
• 单样本形式(one sample)资料检验：
– 已知样本均数与总体均数的比较 – 配对资料的t检验
例3-1正常成年男子的血液中血小板计数的均数为 225×109/L，今测得20名男性油漆作业工人的血小板 计数值（单位109/L）：220，188，162，230，145， 160，238，188，247，113 126，245，164，231， 256，183，190，158，224，175。问油漆作业工人的 血小板计数与正常成年男子有无显著差异？
导入excel数据的t检验
• 临床新药数据 • 入组年龄基线水平是否一致？
SAS Online
• SAS/STAT - SAS/STAT User's Guide
– The TTEST Procedure
• Syntax • Example
Data a; input x1 x2@@; Cards;
113 138 120 116 138 125 120 136 100 110 118 132 138 130 123 110
; proc ttest; Paired x1*x2; run;
另:
Paired x1*x2 Paired x1*x2 a*b Paired (x1-x2)*(b1-b2) Paired (a b)(b c)
• • • • • • • • •
程序： DATA EX3_1; INPUT X@@; CARDS; 220 188 162 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;
Data exp4_2; Input x@@; Cards; 504 498 496 487 509 476 482 510 469 472 ; proc ttest h0=500;/*指定总体均数*/ run;
x 0 t s/ n
The TTEST Procedure
Statistics Lower CL Upper CL Lower CL Upper CL N Mean Mean Mean Std Dev Std Dev Std Dev Std Err 10 479.4 490.3 T-Tests Variable x DF t Value Pr > |t| 9 -2.01 0.0750 501.2 10.482 15.239 27.821 4.8191
DATA EX3_3; DO C=1 TO 2; DO I=1 TO 20; INPUT X@@; OUTPUT; END; END; CARDS; -0.70 -5.60 2.00 2.80 0.70 3.50 2.50 -1.60 1.70 3.00 0.40 4.50 3.70 6.50 5.00 5.20 0.80 0.20 6.00 3.80 2.00 1.60 2.00 2.20 ; PROC UNIVARIATE NORMAL; VAR X; BY C; PROC TTEST COCHRAN; VAR X; CLASS C; RUN;
-15.63 -3.375
3 成组设计(completely random design)两 样本均数比较的t检验
例3-3 为研究国产四类新药阿卡波糖胶囊的降血糖效果，某医 院用40名II型糖尿病病人进行同期随机对照试验。试验者将 这些病人随机等分到试验组（阿卡波糖胶囊组）和对照组 （拜唐苹胶囊组），分别测得试验开始前和8周后的空腹血糖， 算得空腹血糖下降值如表3-2，能否认为该国产四类新药阿卡 波糖胶囊与拜唐苹胶囊对空腹血糖的降糖效果不同？
Variable x
2 配对设计资料的t检验
表 3-1 铁剂和饮食两种方法治疗后患者血红蛋白值（g/L） 患者对子号 1 2 3 4 5 6 7 8 铁剂治疗组 113 120 138 120 100 118 138 123 饮食治疗组 138 116 125 136 110 132 130 110
TTest过程：所有t检验资料
格式：*proc ttest options ( H0=mu/*单样本t-test*/) ; class variable;/*两样本t检验不可少*/ Paired variables;/*配对t检验*/
var variables; By variables; freq variable; Weight variable; run;
4.00 4.60 0.60 1.20
5.80 2.50 3.40 3.10
7.10 6.00 6.60 1.70
-0.50 -1.40 -1.10 -2.00
• 循环语句“DO-END”对变量进行赋值，并 通过输出语句OUTPUT将所赋的值输出到 数据集，以方便数据集的建立 . • 此处为一个双循环嵌套语句，由一个外循 环嵌套一个内循环。 • 从“DO I=1 TO 20”到第一个END为内循 环； • 从“DO C=1 TO 2”到第二个END为外循环。
– – – – The MEANS Procedure The UNIVARIATE Procedure The CHART Procedure The PLOT Procedure
Base SAS
– The GCHART Procedure – The GPLOT Procedure
Std Err 0.6843 0.5412 0.8724
没有原始观测值的样本均数与总体均数 比较的t检验
• 例3-4 某药厂制剂车间用自动装瓶机封装药 液，在装瓶机工作正常时，每瓶药液净重 500克。某日随机抽取10瓶成品，得净重量 的均数为490.3克，标准差为15.2392，问 此时的装瓶机工作是否正常?
The TTEST Procedure Statistics Lower CL Upper CL Lower CL Upper CL Variable C N Mean Mean Mean Std Dev Std Dev Std Dev X 1 20 0.6328 2.065 3.4972 2.3272 3.0601 4.4696 X 2 20 1.4922 2.625 3.7578 1.8408 2.4205 3.5353 X Diff (1-2) -2.326 -0.56 1.2062 2.2547 2.7589 3.5556 (第四部分) T-Tests Variable Method Variances DF t Value Pr > |t| X Pooled Equal 38 -0.64 0.5248 X Satterthwaite Unequal 36.1 -0.64 0.5250 X Cochran Unequal 19 -0.64 0.5286 （第五部分） Equality of Variances Variable Method Num DF Den DF F Value Pr > F X Folded F 19 19 1.60 0.315验(t test)、u(z)检验、秩和检验方法 回顾与适用资料特点 • UNIVARIATE 过程与t 检验 • TTest 过程
• NPAR1WAY 过程
t检验与u检验
• 正态性(Normality)：正态性检验 • 方差相等(Homogenity, Equal variance)：方差齐性 检验 • 独立性(Indenpendenct)： 如果不满足前两个条件，非参数检验。若两样本比 较，仅方差不齐，可考虑t’检验 大样本时可选u(Z)检验近似，但软件不能直接求得， 需写程序。
表 3-2 试验组和对照组空腹血糖下降值(mmol/L) 试验组 (n1=20) 对照组 (n2=20) -0.70 2.50 3.70 6.00 -5.60 -1.60 6.50 3.80 2.00 1.70 5.00 2.00 2.80 3.00 5.20 1.60 0.70 0.40 0.80 2.00 3.50 4.50 0.20 2.20 4.00 4.60 0.60 1.20 5.80 2.50 3.40 3.10 7.10 6.00 6.60 1.70 -0.50 -1.40 -1.10 -2.00
• • • • • • • • • • •
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;
P358数据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 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 ;
配对资料的t检验
表 3-1 铁剂和饮食两种方法治疗后患者血红蛋白值（g/L） 患者对子号 1 2 3 4 5 6 7 8 铁剂治疗组 113 120 138 120 100 118 138 123 饮食治疗组 138 116 125 136 110 132 130 110
• 程序 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;
x1-x2 a-b x1-b1 x1-b2 x2-b1 x2-b2 a-b a-c b-c
The TTEST Procedure
Statistics Lower CL Difference x1 - x2 N 8 Mean Upper CL Lower CL Mean T-Tests Difference x1 - x2 DF 7 t Value -0.65 Pr > |t| 0.5357 Mean 8.8789 9.6911 Upper CL Std Dev Std Err 29.832 5.1822 14.657 Std Dev Std Dev
数据处理技术
2015-2016-1
三 t检验 u检验 秩和检验
东南大学公共卫生学院 孙金芳 Email: sas2015seu@ Passwords: sas2015
review
• a POWERFUL learning tool and makes it easier to find relevant guidance：
U(Z)检验
• 例3-5 青海省20名男子153人，体重均数为 57.38kg。标准差为5.78kg。内蒙古自治区 同年龄的男子206人，体重均数为55.92kg， 标准差为5.17kg，问内蒙古与青海省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;
Univariate 过程与t 检验
• 单样本形式(one sample)资料检验：
– 已知样本均数与总体均数的比较 – 配对资料的t检验
例3-1正常成年男子的血液中血小板计数的均数为 225×109/L，今测得20名男性油漆作业工人的血小板 计数值（单位109/L）：220，188，162，230，145， 160，238，188，247，113 126，245，164，231， 256，183，190，158，224，175。问油漆作业工人的 血小板计数与正常成年男子有无显著差异？
导入excel数据的t检验
• 临床新药数据 • 入组年龄基线水平是否一致？
SAS Online
• SAS/STAT - SAS/STAT User's Guide
– The TTEST Procedure
• Syntax • Example
Data a; input x1 x2@@; Cards;
113 138 120 116 138 125 120 136 100 110 118 132 138 130 123 110
; proc ttest; Paired x1*x2; run;
另:
Paired x1*x2 Paired x1*x2 a*b Paired (x1-x2)*(b1-b2) Paired (a b)(b c)
• • • • • • • • •
程序： DATA EX3_1; INPUT X@@; CARDS; 220 188 162 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;
Data exp4_2; Input x@@; Cards; 504 498 496 487 509 476 482 510 469 472 ; proc ttest h0=500;/*指定总体均数*/ run;
x 0 t s/ n
The TTEST Procedure
Statistics Lower CL Upper CL Lower CL Upper CL N Mean Mean Mean Std Dev Std Dev Std Dev Std Err 10 479.4 490.3 T-Tests Variable x DF t Value Pr > |t| 9 -2.01 0.0750 501.2 10.482 15.239 27.821 4.8191
DATA EX3_3; DO C=1 TO 2; DO I=1 TO 20; INPUT X@@; OUTPUT; END; END; CARDS; -0.70 -5.60 2.00 2.80 0.70 3.50 2.50 -1.60 1.70 3.00 0.40 4.50 3.70 6.50 5.00 5.20 0.80 0.20 6.00 3.80 2.00 1.60 2.00 2.20 ; PROC UNIVARIATE NORMAL; VAR X; BY C; PROC TTEST COCHRAN; VAR X; CLASS C; RUN;
-15.63 -3.375
3 成组设计(completely random design)两 样本均数比较的t检验
例3-3 为研究国产四类新药阿卡波糖胶囊的降血糖效果，某医 院用40名II型糖尿病病人进行同期随机对照试验。试验者将 这些病人随机等分到试验组（阿卡波糖胶囊组）和对照组 （拜唐苹胶囊组），分别测得试验开始前和8周后的空腹血糖， 算得空腹血糖下降值如表3-2，能否认为该国产四类新药阿卡 波糖胶囊与拜唐苹胶囊对空腹血糖的降糖效果不同？
Variable x
2 配对设计资料的t检验
表 3-1 铁剂和饮食两种方法治疗后患者血红蛋白值（g/L） 患者对子号 1 2 3 4 5 6 7 8 铁剂治疗组 113 120 138 120 100 118 138 123 饮食治疗组 138 116 125 136 110 132 130 110
TTest过程：所有t检验资料
格式：*proc ttest options ( H0=mu/*单样本t-test*/) ; class variable;/*两样本t检验不可少*/ Paired variables;/*配对t检验*/
var variables; By variables; freq variable; Weight variable; run;
4.00 4.60 0.60 1.20
5.80 2.50 3.40 3.10
7.10 6.00 6.60 1.70
-0.50 -1.40 -1.10 -2.00
• 循环语句“DO-END”对变量进行赋值，并 通过输出语句OUTPUT将所赋的值输出到 数据集，以方便数据集的建立 . • 此处为一个双循环嵌套语句，由一个外循 环嵌套一个内循环。 • 从“DO I=1 TO 20”到第一个END为内循 环； • 从“DO C=1 TO 2”到第二个END为外循环。