第9讲 定量资料的统计推断（讲）.ppt

样本均数的比较 检验分类 检验 基本命令： Ttest过程： Proc ttest options; Class variable;/*分组语句，必不可少*/ Var variable; By variable; Run; Proc ttest options; 选择项为data=sasdatasets; Class语句： Class variable; 该变量为字符型（16字节）或者数值型均可，但只能有两个分类。 Var 语句； 列出要比较分析的变量，如果不列出，则用在其他语句中未出现的数值型变量分析。 By语句； By variable; 按by 语句的变量分组进行统计推断；by语句使用前要进行排序，或者用notsorted 与descending T检验分类 单样本均数t检验 两独立样本t检验 配对样本t检验 几何均数t检验（略） 等效性检验 可信区间的计算 单样本均数t检验(one-sample test) Data one_sample test; N=36; S_m=130.83; Std=25.74; P_m=140; Df=n-1; T-(s_m-p_m)/(std/sqrt(n)); P=(1-probt(abs(t),df))*2; Proc print; Var t p; Run; 两独立样本t检验示例：stscores程序 Data stscores; input sex$ score@@; Cards; F 75 f 76 f 80 f 77 f 80 f 77 f 73 m 82 m 80 m 85 m 85 m 78 m 87 m 82 ; proc ttest; class sex; var score; run; 配对样本t检验 先求各对子差值，然后对差值是否为零，运用proc univariate完成。 示例：程序paired t test（univariate） data work1; input f m @@; dif=m-f; datalines; 169 170 102 96 67 74 206 178 109 86 97 92 156 179 73 92 83 87 ; proc univariate data=work1; var dif; title '成对差值T检验'; run; Paired ttest (means) data a; input id pretest posttest@@; dif=posttest-pretest; cards; 1 80 82 2 73 71 3 70 95 4 60 69 5 88 100 6 84 71 7 65 75 8 37 80 9 91 95 10 98 99 11 52 65 12 78 83 13 40 60 14 79 86 15 59 62 ; proc means mean std t prt; var dif; title 'paired conparisons t test'; run; 等效性检验 data euqal_test; n1=102;n2=100; m1=0.57;m2=0.48; s1=0.89;s2=0.82; delta=0.52; sc2=(s1**2*(n1-1)+ s2**2*(n2-1))/(n1+n2-2); 单一总体均数的可信区间计算 Data onesample_cl; N=10; Mean=166.95; Std=3.64; T=tinv(0.975,n-1); In=t*std/sqrt(n); Lclm=mean-in; Uclm=mean+in; Run; Proc print; Var lclm ulclm; Run; 两个总体均数差的可信区间 方差分析 可用于方差分析的方法包括：GLM ANOVA NESTED VARCOMP GLM：进行方差分析、回归、协方差分析、多元方差分析. ANOVA：对平衡设计进行方差分析 NESTED：纯嵌套随机模型的方差分析 VARCOMP：估计方差成分 ANOVA 用于平衡设计 对实验的每种组合安排相同的实验单位，称为平衡设计。 研究设计类型 完全随机设计 随机区组设计 拉丁方设计 析因设计 正交设计 系统设计 裂区设计 命令： Proc anova options; Class variables; Model dependent=effects/options; Means effects/options; Absorb variables; Freq variables;