STATA软件应用三定量资料的统计分析

医学统计学实习课,STATA软件的应用(3)定量资料的统计分析,定量资料的统计分析,样本均数与总体均数比较的t检验配对设计t检验成组设计t检验单因素方差分析两因素方差分析方差齐性检验正态性检验,样本均数与总体均数比较的t检验,仅有数据之统计量时ttesti#obs#mean#sd#val仅有原始数据时ttest变量名=#val注：#val/*总体均数命令可以加一些if和in条件限制,样本均数与总体均数比较的t检验,例4.4命令：ttesti255.10.884.6,样本均数与总体均数比较的t检验,例：某医生随机抽查10名某病患者的血红蛋白，如下：问：该病患者的平均Hb含量是否与正常人的平均Hb含量相同(正常人的平均Hb含量为14.02(g/dl)。data:ttest_1,问题：,正常人患者,？,统计量与参数不同的两种可能,其一：抽样误差(偶然的、随机的、较小的)其二：本质上的差别(必然的、大于随机误差),样本均数与总体均数比较的t检验,ttestx=14.02One-samplettest-Variable|ObsMeanStd.Err.Std.Dev.95%Conf.Interval-+-x|1012.59.51627941.63261911.4220913.75791-mean=mean(x)t=-2.7698Ho:mean=14.02degreesoffreedom=9Ha:mean14.02Pr(T|t|)=0.0218Pr(Tt)=0.9891,配对设计t检验,ttest变量1=变量2变量1与变量2必须成对输入，样本含量必须相等，如有缺省值，则用小数点表示，但与之对应的记录在计算时被忽略,配对设计t检验,例：10例矽肺患者经某药治疗，其血红蛋白(g/dl)如下：data:ttest_2,配对设计t检验,.ttestx1=x2Pairedttest-Variable|ObsMeanStd.Err.Std.Dev.95%Conf.Interval-+-x1|1012.59.51627941.63261911.4220913.75791x2|1013.27.34158131.08017512.4972914.04271-+-diff|10-.6799999.52042721.645735-1.857288.4972881-mean(diff)=mean(x1-x2)t=-1.3066Ho:mean(diff)=0degreesoffreedom=9Ha:mean(diff)0Pr(T|t|)=0.2237Pr(Tt)=0.8881,配对设计t检验,另一种思路：gend=x0-x1ttestd=0结果一致,成组设计t检验,ttest变量1=变量2,unpairedunequalttest变量,by（分组变量）unequalttesti#obs1#mean1#sd1#obs2#mean2#sd2,unequalunpaired表示非配对的，如不选就作配对检验unequal表示假设两组方差不齐，如不选表示假设两组方差达到齐性,成组设计t检验,例：分别测得14例老年人煤饼病人及11例正常人的尿中17酮类固醇排出量(mg/dl)如下，试比较两组的均数有无差别病人：2.905.415.484.604.035.104.974.244.362.722.372.097.105.92健康人：5.188.793.146.463.726.645.604.577.714.994.01data:ttest_3.dtattest_4.dta,成组设计t检验,ttestx1=x2,unpairedttestx,by(g),独立性正态性方差齐性方差齐：t检验方差不齐:t检验,均数的假设检验应用条件,两组资料间的方差齐性检验,sdtest变量名1=变量名2sdtest变量,by（分组变量）sdtesti#obs1#mean1#sd1#obs2#mean2#sd2,sdtestx,by(g)Varianceratiotest-Group|ObsMeanStd.Err.Std.Dev.95%Conf.Interval-+-1|144.377857.38751.4498923.5407145.2152|115.528182.52324311.7354014.3623246.69404-+-combined|254.884.33064531.6532274.2015825.566418-ratio=sd(1)/sd(2)f=0.6980Ho:ratio=1degreesoffreedom=13,10Ha:ratio1Pr(Ff)=0.7325,单因素方差分析,oneway因变量分组变量，选择项noanova/*不打印方差分析表missing/*将缺省值作为单独的一组tabulate/*打印各组的基本统计量表简写：tscheffe/*Scheffe法简写：schBonferroni/*Bonferroni法简写：bonsidak/*Sidak法简写：si,单因素方差分析,三组小鼠的FDP酶活力,单因素方差分析,onewayxgrouponewayxgroup,tonewayxgroup,tschdata:oneway.dta,.onewayxg,tsch|Summaryofxg|MeanStd.Dev.Freq.-+-1|2.7025.5001356982|3.66125.9850806983|3.90125011.01644258-+-Total|3.4216667.9827320724,AnalysisofVarianceSourceSSdfMSFProbF-Betweengroups6.4368089723.218404494.280.0275Withingroups15.775724621.751224983-Total22.212533623.965762331Bartlettstestforequalvariances:chi2(2)=3.4559Probchi2=0.178,Comparisonofxbyg(Scheffe)RowMean-|ColMean|12-+-2|.95875|0.111|3|1.19875.24|0.0380.859,两因素方差分析,四种抗癌药物抑瘤效果,anova因变量分组变量1分组变量2,两因素的方差分析,egena=seq(),b(5)egenb=seq(),f(1)t(5)anovaxabdata：anova.dta,两因素的方差分析,.anovaxgroupblockNumberofobs=20R-squared=0.7058RootMSE=.134818AdjR-squared=0.5341Source|PartialSSdfMSFProbF-+-Model|.5231700047.0747385724.110.0157|group|.4108399983.1369466667.530.0043block|.1123300064.0280825011.550.2514|Residual|.21811000112.018175833-+-Total|.74128000519.039014737,正态性检验,命令：sktest变量该命令要求资料的样本含量至少为8Skewness/KurtosistestsforNormality-joint-Variable|Pr(Skewness)Pr(Kurtosis)adjchi2(2)Probchi2-+-db|0.3040.3042.530.2829,小结,样本均数与总体均数比较的t检验ttestx=14.02配对设计检验ttestx1=x2成组设计检验及方差齐性检验sdtestx1=x2ttestx1=x2,unpairedunequalttestx,by（group）unequal单因素方差分析onewayxg,tsch两因素方差分析anovaxab,练习,例4.5（P55）:配对t检验例4.7（P58）：成组t检验例5.1（P68）：单因素方差分析,思考题,某医师为评价两种药物对血沉的影响，用药前后结果如下，问两种药物对血沉影响有无差别。A药组B药组前后前后403639204335402536303930412943234030382437283620382