《SAS置信区间t检验》PPT课件.ppt

《《SAS置信区间t检验》PPT课件.ppt》由会员分享，可在线阅读，更多相关《《SAS置信区间t检验》PPT课件.ppt（20页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 置信区间和置信区间和t检验检验总体均数的可信区间总体均数的可信区间t t 分布法分布法 SASSAS函数：函数：TINVTINV t=t=TINV(p,dfTINV(p,df)可求可求t t分位数的函数，分位数的函数，p p是从是从-到当前到当前t t分位数位置的面积。分位数位置的面积。dfdf=n-1=n-1（自由度）自由度）例：随机抽取例：随机抽取15名学生，记录他们的性别名学生，记录他们的性别（sex）、）、年龄（年龄（age）、体重（）、体重（w）和身高（和身高（h），），求：学生身高求：学生身高的的95%置信区间。置信区间。f 15 46 156 f 14 41 149 f 18

2、65 165m 15 50 160 m 13 48 155 m 18 70 180m 14 38 150 m 16 55 165 m 17 68 176m 16 60 170 f 17 50 160 f 17 58 160f 16 60 165 m 17 65 175 f 18 61 162data a;input sex$age w h;cards;f 15 46 156 f 14 41 149 f 18 65 165m 15 50 160 m 13 48 155 m 18 70 180m 14 38 150 m 16 55 165 m 17 68 176m 16 60 170 f 17 5

3、0 160 f 17 58 160F 16 60 165 m 17 65 175 f 18 61 162；proc means;var h;Output out=b mean=mean std=s n=n;data c;set b;t=tinv(0.975,n-1);/*也可使用也可使用t=tinv(0.025,n-1);但此时计算出来的但此时计算出来的t界界值为负数。值为负数。*/yl=mean-t*s/sqrt(n);/*可信区间的下限可信区间的下限*/y2=mean+t*s/sqrt(n);/*可信区间的上限可信区间的上限*/proc print;run;单样本资料的 t检验 目的目的:

4、根据样本均数根据样本均数 推断其总体均数推断其总体均数 是否与已知是否与已知总体均数总体均数 0相等相等 与与 0的比较。的比较。应用条件应用条件:独立性独立性 正态性正态性例1 某镇痛新药在药厂投入量产，已知该药某镇痛新药在药厂投入量产，已知该药的生产技术标准为：每片药中平均有效的生产技术标准为：每片药中平均有效药物含量为药物含量为40.4mg40.4mg。为评价生产质量是为评价生产质量是否达标，现抽测否达标，现抽测1010片药物，结果为：片药物，结果为：41.341.3，40.840.8，38.738.7，42.742.7，43.043.0，38.638.6，42.742.7，39.539

5、.5，42.142.1，42.142.1，问该厂是，问该厂是否达到生产要求？否达到生产要求？data aa;input a;x=a-40.4;cards;41.3 40.8 38.7 42.7 43.0 38.6 42.7 39.5 42.1 42.1;数据步数据步 proc means mean std stderr t prt;var x;run;proc univariate normal;var x;run;proc ttest H0=40.4;var a;run;方法方法1方法方法2方法方法3方法1 方法2方法3结果分析 本例资料服从正态分布本例资料服从正态分布（W=0.876128

6、，P=0.1178），故选），故选t检验。本例检验。本例t=1.41，双侧检验，双侧检验p=0.1917，按，按a=0.05水准，不拒绝水准，不拒绝H0，差别无统计学意义，尚不认为该厂生产，差别无统计学意义，尚不认为该厂生产质量不达标。质量不达标。配对资料t检验资料类型：资料类型：两个同质对象接受不同处理；两个同质对象接受不同处理；同一受试对象分别接受不同的处理，同一同一受试对象分别接受不同的处理，同一受试对象处理前后。受试对象处理前后。条件：条件：差值差值d服从正态分布服从正态分布例例2数据步数据步data aa;input x1 x2;d=x1-x2;cards;0.940.92 1.02

7、1.01 1.141.111.231.22 1.311.32 1.411.421.531.51 1.611.61 1.721.721.811.82 1.931.93 2.022.04;proc means t prt；var d；run；proc univariate normal；var d；run；proc ttest;paired x1*x2;run;方法1方法2方法3方法一 方法二方法三结果解释 本本例例t0.771，P=0.4569，在，在=0.05水平水平上不能拒绝上不能拒绝H0，差别没有统计学意义，差别没有统计学意义，尚不能认为两法测定结果不同。尚不能认为两法测定结果不同。作业作业作业