医学STATA软件及其应用PPT课件

.,1,主要内容：定量资料的假设检验,t检验单样本t检验配对t检验成组t检验及方差齐性检验方差分析单因素方差分析及两两比较两因素以上的方差分析,.,2,样本均数与总体均数比较的t检验,有原始数据时ttest变量名=#val仅有数据之统计量时ttesti#obs#mean#sd#val,.,3,例1(Ex_1),样本：某医生随机抽查10名某病患者的血红蛋白，如下：问题：该病患者的平均Hb含量是否与正常人的平均Hb含量相同(正常人的平均Hb含量为14.02(g/dl)。,.,4,例,样本：某医生随机抽查10名某病患者的血红蛋白，求得其均数为12.59(g/dl),标准差为1.632619(g/dl)。问题：该病患者的平均Hb含量是否与正常人的平均Hb含量相同(正常人的平均Hb含量为14.02(g/dl)。,.,5,操作,ttest变量名=#valttestx=14.02ttesti#obs#mean#sd#valttesti1012.591.63261914.02,.,6,One-samplettest-|ObsMeanStd.Err.Std.Dev.95%Conf.Interval-+-x|1012.59.51627951.63261911.4220913.75791-Degreesoffreedom:9Ho:mean(x)=14.02Ha:mean14.02t=-2.7698t=-2.7698t=-2.7698P|t|=0.0218Pt=0.9891,P(|t|2.7698)=0.0218,假设检验与可信区间的分析具有相同的量化判定，其结论是等价的,ttesti1012.51.63261914.02,.,7,配对设计检验,ttest变量1=变量2注意：变量1和变量2必须是成对输入的配对样本昆明治疗羊角风,.,8,例3(配对设计)(Ex_3),例3现用两种测量肺活量的仪器对12名妇女测得最大呼气率(PEER)(L/min)，资料如下，问两种方法的检测结果有无差别？（P60）,.,9,ttestx2=x1,Pairedttest-Variable|ObsMeanStd.Err.Std.Dev.95%Conf.Interval-+-x2|12419.333326.0085590.09624362.0889476.5778x1|12402.166727.4159494.97161341.8246462.5087-+-diff|1217.1666711.6416440.32782-8.4564142.78974-Ho:mean(x2-x1)=mean(diff)=0Ha:mean(diff)0t=1.4746t=1.4746t=1.4746P|t|=0.1684Pt=0.0842,ttestx1=x2,是否还有其他方法?,.,10,独立性正态性方差齐性方差齐：t检验方差不齐:t检验,t检验应用条件,.,11,两组资料间的方差齐性检验,sdtest变量名1=变量名2/*两组数据各用一个变量sdtest变量名，by（分组变量）/*两组数据用同一个变量输入，用分组变量加以区分sdtesti#obs1#mean1#sd1#obs2#mean2#sd2/*mean1,mean2可被替换为点号。,.,12,成组设计检验,ttest变量1=变量2,unpairedunequalwelchttest变量,by（分组变量）unequalwelchttesti#obs1#mean1#sd1#obs2#mean2#sd2,unequalwelchunpaired表示非配对的，如不选就作配对检验unequal表示假设两组方差不齐，如不选表示假设两组方差达到齐性welch当方差不齐时，用welch法比较注：welch必须和unequal同时选才有效,.,13,例（成组设计）(Ex_5),某医生研究转铁蛋白对病毒性肝炎诊断的临床意义，测得12名正常人和15名病毒性肝炎患者血清转铁蛋白含量(g/dl)，结果如下，问患者和正常人转铁蛋白含量是否有差异?（同例）,.,14,Varianceratiotest-Group|ObsMeanStd.Err.Std.Dev.95%Conf.Interval-+-1|12271.89172.99558910.37702265.2984278.48492|15235.20673.71566214.3907227.2374243.176-+-combined|27251.51114.31252922.40856242.6466260.3756-Ho:sd(1)=sd(2)F(11,14)observed=F_obs=0.520F(11,14)lowertail=F_L=F_obs=0.520F(11,14)uppertail=F_U=1/F_obs=1.923Ha:sd(1)sd(2)PF_U=0.2649PF_obs=0.8596,Sdtestx,by(group),.,15,ttestx,by(group)Two-samplettestwithequalvariances-Group|ObsMeanStd.Err.Std.Dev.95%Conf.Interval-+-1|12271.89172.99558910.37702265.2984278.48492|15235.20673.71566214.3907227.2374243.176-+-combined|27251.51114.31252922.40856242.6466260.3756-+-diff|36.6854.9500326.4902246.87978-Degreesoffreedom:25Ho:mean(1)-mean(2)=diff=0Ha:diff0t=7.4111t=7.4111t=7.4111P|t|=0.0000Pt=0.0000.,P（|t|7.4111）=0.0000,.,16,思考,t检验两个均数比较t检验那么，三个、四个均数如何比较？,方差齐,方差不齐,？,方差分析,.,17,单因素方差分析及方差齐性检验,oneway响应变量分组变量，选择项noanova/*不打印方差分析表missing/*将缺省值作为单独的一组tabulate/*打印各组的基本统计量表scheffe/*Scheffe法bonferroni/*Bonferroni法sidak/*Sidak法注意：方差分析时数据必须按特定形式输入，即新定义一分组变量g，而将不同处理组的观察值以某变量x表示。,.,18,例6.4三组小鼠的FDP酶活力,.,19,资料输入形式,XgABC2.7912.793.835.412.6912.6.8323.152.,Stackabc,into(x),.,20,onewayxg,tsch,|Summaryofxg|MeanStd.Dev.Freq.-+-1|2.7025.5001356982|3.66125.9850806983|3.90125011.01644258-+-Total|3.4216667.9827320724AnalysisofVarianceSourceSSdfMSFProbF-Betweengroups6.4368089723.218404494.280.0275Withingroups15.775724621.751224983-Total22.212533623.965762331Bartlettstestforequalvariances:chi2(2)=3.4559Probchi2=0.178,.,21,两两比较结果,Comparisonofxbyg(Scheffe)RowMean-|ColMean|12-+-2|.95875|0.111|3|1.19875.24|0