03 SAS统计分析-T检验

1t检验邓伟2012.122t检验简介均数的抽样分布规律:来自正态总体或来自偏态总体但样本量比较大的样本其均数服从正态分布标准变换后u统计量:u~N(0,1),以样本标准差s代替总体标准差σ,此时变换后的统计量称为t统计量,t统计量服从自由度为ν-1的t分布t检验的种类单样本均数的比较：样本所代表的总体均数和某已知总体均数比较配对设计两样本的比较：每对差值的总体均数是否为0完全随机设计两样本：两样本所代表的总体均数是否相等/xxxtSSn/xxxun2~(,)xxN3配对t检验原理:检验差值d的总体均数是否为0假设:公式:常见配对的情况同一对象治疗（或处理）前后的比较同一样本接受不同处理的比较配对的两个受试对象分别给予两种处理01:0:0ddHH0//ddddtsnsn4datattest1;inputx1x2;/*输入治疗前、后的血红蛋白值*/d=x1-x2;/*计算治疗前后的差值*/cards;113140150138150140135135128135100120110147120114130138123120;procmeansmeanstderrtprt;/*输出对均数是否为0的t检验及其p值*/vard;run;配对t检验5AnalysisVariable:DMeanStdErrorTProb|T|---------------------------------------------------6.80000005.2042717-1.30661890.2237--------------------------------------------------配对t检验6成组t检验原理:两样本所代表的总体均数是否相等,检验两个样本均数的差是否为0假设:公式121201:0:0xxxxHH121201::xxxxHH12121212()0xxxxxxxxtSS7成组t检验应用条件独立性正态性－－univariate正态性检验方差齐性－－齐性检验－t‘检验8datattest3;doc=1to2;/*c循环开始,循环次数为2*/inputn;/*输入n值*/doi=1ton;/*i循环开始,循环次数为n*/inputx@@;/*输入x值*/output;/*将x输出*/end;/*i循环结束*/end;/*c循环结束*/cards;52793343033381983229274310;procprint;procttest;classc;varx;run;成组t检验9循环语句的应用do………;/*循环1开始*/do………;/*循环2开始*/input……;output;end;/*循环2结束*/end;/*循环1结束*/注意：1）@@的合理应用2）output的作用：10/*循环的结果*/OBSCNIX11512792152334315330341543385155198623122972322748233310成组t检验11TTESTPROCEDUREVariable:XCNMeanStdDevStdError---------------------------------------------------------------15290.4000000056.9938593225.4884287523271.0000000040.5832477823.43074903VariancesTDFProb|T|---------------------------------------Unequal0.56035.60.5971Equal0.50996.00.6284ForH0:Variancesareequal,F'=1.97DF=(4,2)ProbF'=0.7272成组t检验12单样本ttest原理:检验某样本所来自的总体均数与已知总体均数是否相等假设:公式:例:某疾病的平均住院天数是否是80天？xxtS01::xxHH13/*程序1*/datatime;inputtime@@;cards;439084871169586999392121716698791026011210598;procttesth0=80;vartime;run;/*程序2*/Datatime;Settime;D=time-80;Procmeansstderrtprt;Run;单样本ttest14SAS程序:procttestProcttest[options];/*过程名*/Classvar;/*分组变量名,决定哪两组进行比较*/Varlist;/*反应变量列表*/PAIREDvar1*var2;/*专用于配对变量检验,class和var语句在使用该语句时无效*/Byvariables;/*指明分层变量,按该变量分别分析,需先排序*/Freqvariable;/*指明频数变量*/Run;15optionsAlpha＝定义可信区间的alpha水平，默认为0.05H0＝定义特定的总体均数值，默认为0by：根据by变量分别计算分析结果，要求事先根据by变量排序class：定义分组变量，决定进行哪两组比较。在ttest中只允许有两组SAS程序:procttest16Freq－weight非常接近Freq－－频数，只能是整数Weight－－可以是小数SAS程序:procttest17语句作用PAIREDA*BA-BPAIREDA*BC*DA-BandC-DPAIRED(AB)*(CD)A-C,A-D,B-C,andB-DPAIRED(AB)*(CB)A-C,A-B,andB-CPAIRED(A1-A2)*(B1-B2)A1-B1,A1-B2,A2-B1,andA2-B2PAIRED(A1-A2):(B1-B2)A1-B1andA2-B2Paired－－定义配对变量SAS程序:procttest18例1的分析用ttest过程解决:Procttestdata=ttest1;Pairedx1*x2;Run;LowerCLUpperCLLowerCLUpperCLDifferenceNMeanMeanMeanStdDevStdDevStdDevStdErrx1-x210-18.57-6.84.972911.3216.45730.0455.2043T-TestsDifferenceDFtValuePr|t|x1-x29-1.310.223719SAS的循环语句介绍记数循环语法do记数变量＝起始值to终止值by增减值…...end；条件循环语句Dowhile（条件）；……End；Dountil（条件）；……End；20记数循环语法dataa;amount=1000;rate=0.08;year_int=amount*rate;doi=1to12;mont_int+(mont_int+amount)*rate/12;end;run;procprint;run;21条件循环语句datab;amount=1000;rate=0.08;year_int=amount*rate;i=0;dowhile(i12);i=i+1;mont_int+(mont_int+amount)*rate/12;end;procprint;run;22条件循环语句datac;amount=1000;rate=0.08;year_int=amount*rate;i=0;dountil(i=12);i=i+1;mont_int+(mont_int+amount)*rate/12;end;procprint;run;