第十二章有序分类资料的统计分析的Stata实现.doc

《第十二章有序分类资料的统计分析的Stata实现.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、第十二章有序分类资料的统计分析的Stata实现本章使用的STATA命令：列变量有序时的分类资料CMH卡方分析opartchi行变量[weight],by(列变量)（见Stata7附加程序）双向有序时的Spearman相关spearman变量1变量2例12-2某研究欲观察人参的镇静作用，选取32只同批次的小白鼠，将其中20只随机分配到人参组：以5%人参浸液对其做腹腔注射，12只分配到对照组：以等量蒸馏水对其做同样注射。实验结果如表12-2所示。能否说明人参有镇静作用？表12-2人参镇静作用的实验结果镇

2、静等级人参组对照组-411±10+21++10+++1201．建立检验假设，确定检验水准。：人参没有镇静作用（样本来自两个相同总体）：人参有镇静作用（样本来自两个不同总体）Stata数据为：abx11412113214115122111220231240250Stata命令为：opartchib[weight=x],by(a)结果为：Chi-squaretestsdfChi-squareP-valueIndependence416.640.0023--------------------------

3、-----------------------------ComponentsofindependencetestLocation115.290.0001Dispersion1.34960.5543在的水平上，拒绝，接受H1，认为两总体之间的差别有统计学意义，可以认为人参组和对照组镇静等级的差别有统计学意义，人参有镇静作用。例12-3试根据表12-4的资料，检验针刺不同穴位的镇痛效果有无差别？表12-4针刺不同穴位的镇痛效果穴位镇痛效果＋＋＋＋＋＋＋＋＋＋合谷38441224足三里53292816扶

4、突472319331．建立检验假设，确定检验水准。：三个穴位的镇痛效果相同：三个穴位的镇痛效果不全相同Stata数据为：groupeffectw113812441312142421532229232824163147322333193433Stata命令为：opartchieffect[weight=w],by(group)结果为：Chi-squaretestsdfChi-squareP-valueIndependence622.070.0012---------------------------

5、----------------------------ComponentsofindependencetestLocation23.1180.2103Dispersion25.9140.0520有。在的水平上尚不能拒绝，即根据本例资料尚不能认为针刺不同穴位的镇痛效果差别有统计学意义。例12-4两名放射科医师对13张肺部X线片各自做出评定，评定方法是将X线片按病情严重程度给出等级，结果如表12-6所示。问他们的评定结果是否相关。表12-6两名放射科医师对13张肺部X片的评定结果X片编号1234567

6、8910111213甲医师X+++－±－+++++++++++－+++乙医师Y±++++－+++++++++++++±++++1．建立检验假设，确定检验水准。：＝0（两名医师的等级评定结果不相关）：（两名医师的等级评定结果相关）Stata数据为：ixy1212333024125006237348439341044110112331323Stata命令为：spearmanyx结果为：Numberofobs=13Spearman'srho=0.8082TestofHo:yandxareindepende

7、ntProb>

8、t

9、=0.0008P<0.05。在水平上，拒绝H0，接受H1，认为两名射科医师对肺部X片的等级评定结果有正相关关系。例12-5300名抑郁症患者按其抑郁程度和自杀意向的轻重程度的分类数据如表12-8所示。问自杀意向的轻重程度和抑郁程度之间是否存在线性变化趋势？表12-8300名抑郁症患者的分类数据自杀意向（X）抑郁程度（Y）合计轻度中等严重(Y=1)（Y=2)（Y=3）无（X=1）1357314222想要自杀（X=2）591024曾自杀过（X=3）8232354合计148105473

10、001、建立检验假设，确定检验水准。：自杀意向与抑郁症之间不存在线性变化趋势：自杀意向与抑郁症之间存在线性变化趋势Stata数据为：yxw1113512513821732292323311432103323Stata命令和结果为*：.corryx[fw=w](obs=300)

11、yx-------------+------------------y

12、1.0000x

13、0.46421.0000.returnlistscalars:r(N)=300r(rho)=.464184