文本内容:
第五章多组计量资料比较的非参数检验命令与输出结果说明本节STATA命令摘要ranksum观察变量,by（分组变量）kwallis观察变量,by（分组变量）秩和检验（Mann,Whitney andWilcoxon非参数检验）对于计量资料不满足正态分布要求或方差不齐性，但样本资料之间是独立抽取的，则可以应用秩和检验方法进行比较两组资料的中位数是否有差异STATA命令为例研究不同饲料对雌鼠体重增加的关系表中用x表示雌鼠体重增加（克），用group=l表示高蛋白饲料组以及用group=2表示低蛋白饲料组X11911312913414610112416110783group1111111111X97123701181011071329485group222222211无效假设Ho:两组增加体重的中位数相同・ranksum x,bygroupTwo-sample Wi1coxon rank-sum Mann-Whitney test
12140.
5120749.570obs ranksum expectedcombinedunadjusted variance
140.00adjustment forties-
0.12adjustedvariance
139.88Ho:xgroup==l=xgroup==2z=
1.733Prob|z|=
0.0830ranksum x,by groupTest:Equality ofmedians Two-Sample WilcoxonRank-SumSum ofRanks:
49.5group==2Expected Sum:70
②Prob|z|
0.0832z-statistic T.73
①为第二组（低饲料组）的秩的和；
②若效假设成立，则第二组的秩的和期望值为70；
③秩和统计检验量z；
④对于无效假设Ho对应的p值在本例中，虽然第二组的秩和为
49.5而期望值估计为70,但p值为
0.0832,所以根据该资料和统计结果一般不能认为用高蛋白饲料喂养能明显增加雌鼠的体重多组资料中位数比较（完全随机化设计资料的检验）对于完全随机化设计资料的比较，若各组资料不全服从正态分布（即至少有一组的资料均不服从正态分布）或各组的资料方差不齐性，则可以用Kruskal andWallis方法进行检验（Ho:各组的中位数相同）STATA命令为kwallis观察变量,by（分组变量）例50只小鼠随机分配到5个不同饲料组，每组10只小鼠在喂养一定时间后，测定鼠肝中的铁的含量（mg/g）如表所示试比较各组鼠肝中铁的含量是否有显著性差别用x表示鼠肝中铁的含量以及用group=1,2,3,4,5分别表示对应的5个组
2.
231.
142.
631.
001.
352.
011.
641.
131.
011.70X1111111111group
5.
590.
966.
961.
231.
612.
941.
963.
681.
542.59Xgroup
22222222224.
53.
9210.
338.
232.
074.
96.
846.
423.726Xgroup
33333333331.
351.
060.
740.
961.
162.
080.
690.
680.
841.34Xgroup
44444444441.
41.
512.
491.
741.
591.
3634.
815.
215.12Xgroup5555555555,kwal1is x,bygroupKruskal-Wallis equality-of-populations ranktestgroup obsRank Sum
10188.
5010280.
5010420.
001095.
0010291.00chi-squared=
27.856with4d.f.probability=.
0.0001chi-squared withties=
27.858with4d.f.probability=.
0.0001kwallis x,by groupprobability=
0.0001Test:Equality ofpopulations Kruskal-Wallis Test
①group Obs_RankSum
110188.
50210280.
50310420.
0041095.
00510291.00chi-squared=
27.856with4d.f.
②①为各组的秩和值；
②为该统计量的c2检验值；
③为无效假设检验所对应的值本例结果表明5组的中位数有显著的差异即5个不同饲料组的小鼠肝脏中铁的含量有显著差异，说明小鼠肝脏中铁的含量与喂养的饲料有关。
个人认证
优秀文档
获得点赞 0
