Stata基本操作和数据分析入门：第四讲 两组计量资料平均水平的统计检验.doc

《Stata基本操作和数据分析入门：第四讲 两组计量资料平均水平的统计检验.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、Stata软件基本操作和数据分析入门第四讲两组计量资料平均水平的统计检验一、配对设计的平均水平检验统计方法选择原则：如果配对的差值服从近似正态分布(小样本)或大样本，则用配对t检验小样本的情况下，配对差值呈明显偏态分布，则用配对秩符号检验(matched-pairssigned-rankstest)。例110例男性矽肺患者经克矽平治疗，其血红蛋白（g/dL）如下：表10例男性矽肺患者血红蛋白值（g/dL）病例号12345678910治疗前11.315.015.013.512.810.011.012.013.012.3治疗后14.01

2、3.814.013.513.512.014.711.413.812.0问：治疗前后的血红蛋白的平均水平有没有改变这是一个典型的前后配对设计的研究(但不提倡，因为对结果的解释可能会有问题)Stata数据输入结构X1X211.3141513.8151413.513.512.813.510121114.71211.41313.812.312操作如下：gend=x1-x2产生配对差值的变量dsktestd正态性检验正态性检验结果如下：.sktestdSkewness/KurtosistestsforNormality-------joint

3、------Variable

4、Pr(Skewness)Pr(Kurtosis)adjchi2(2)Prob>chi2-------------+-------------------------------------------------------d

5、0.2790.7741.430.4885正态性检验的无效假设为：资料正态分布相应的备选假设为：资料非正态分布a=0.05，由于正态性检验的P值=0.40189>>a，故可以认为资料近似服从正态分布。ttestd=0配对t检验:H0:md=0vsH1:md¹0，a=0.05结果如下

6、：One-samplettest------------------------------------------------------------------------------Variable

7、ObsMeanStd.Err.Std.Dev.[95%Conf.Interval]---------+--------------------------------------------------------------------d

8、10-.6799999.52042721.645735-1.857288.4972881-

9、-----------------------------------------------------------------------------Degreesoffreedom:9Ho:mean(d)=0Ha:mean<0Ha:mean~=0Ha:mean>0t=-1.3066t=-1.3066t=-1.3066P

10、t

11、=0.2237P>t=0.8881P值=0.2237>a，故认为治疗前后的血红蛋白的平均数差异没有统计学意义。即：没有足够的证据可以认为治疗前后的血红蛋白的总体平均数不同。如果已知差值

12、的样本量，样本均数和样本标准差，可以用立即命令如下(如，已知样本量为10，差值的样本均数为-0.66，差值的标准差为1.65，则输入命令如下：ttesti样本量样本均数样本标准差0本例为：ttesti10-0.661.650得到下列结果如下：.ttesti10.661.650One-samplettest------------------------------------------------------------------------------

13、ObsMeanStd.Err.Std.Dev.[95%Conf.Inter

14、val]---------+--------------------------------------------------------------------x

15、10.66.52177581.65-.52033891.840339------------------------------------------------------------------------------Degreesoffreedom:9Ho:mean(x)=0Ha:mean<0Ha:mean~=0Ha:mean>0t=1.2649t=1.264

16、9t=1.2649P

17、t

18、=0.2377P>t=0.1188结果解释与结论同上述相同。如果对于小样本的情况下，差值不满足正态分布，则用Match-Sign-ranktest，操作如下：signrank差值变量名=