Stata第五章 多组计量资料比较的非参数检验命令与输出结果说明

2012-04-17 生物谷 生物谷

本节STATA 命令摘要 ranksum观察变量,by(分组变量) kwallis 观察变量,by(分组变量) ·       秩 和检验(Mann,WhitneyandWilcoxon非参数检验)  对于计量资料不满足正态分布要求或方差不齐性,但样本资料之间是独立

本节STATA 命令摘要

ranksum观察变量,by(分组变量)

kwallis 观察变量,by(分组变量)

·        和检验(Mann,WhitneyandWilcoxon非参数检验)

 对于计量资料不满足正态分布要求或方差不齐性,但样本资料之间是独立抽取的,则可以应用秩和检验方法进行比较两组资料的中位数是否有差异。STATA命令为:

ranksum  观察变量,by(分组变量)

 例:研究不同饲料对雌鼠体重增加的关系(摘自医学统计方法,金丕焕主编,p218)。表中用x表示雌鼠体重增加(克),用group=1表示高蛋白饲料组以及用group=2表示低蛋白饲料组。

x

134

146

104

119

124

161

107

83

113

129

group

1

1

1

1

1

1

1

1

1

1

x

97

123

70

118

101

85

107

132

94

 

group

1

1

2

2

2

2

2

2

2

无效假设Ho:两组增加体重的中位数相同。

ranksumx,by(group)

Test:Equalityofmedians(Two-SampleWilcoxonRank-Sum)

                                                                

SumofRanks:49.5(group==2)

ExpectedSum:70                                                                                                  

           

z-statistic-1.73

Prob>|z|0.0832                                                                                                

①为第二组(低饲料组)的秩的和;②若无效假设成立,则第二组的秩的和期望值为70;③秩和统计检验量z;④对于无效假设Ho对应的p值。

 在本例中,虽然第二组的秩和为49.5而期望值估计为70,但p值为0.0832,所以根据该资料和统计结果一般不能认为 用高蛋白饲料喂养能明显增加雌鼠的体重。

·       多组资料中位数比较(完全随机化设计资料的检验)

 对于完全随机化设计资料的比较,若各组资料不全服从正态分布(即:至少有一组的资料均不服从正态分布)或各组的资料方差不齐性,则可以用KruskalandWallis方法进行检验(Ho:各组的中位数相同)。STATA命令为:

kwallis 观察变量,by(分组变量)

例:505组,10鼠。后,(mg/g)示:(法,编,p220)xgroup=123455组。

x:

2.23

1.14

2.63

1

1.35

2.01

1.64

1.13

1.01

1.70

group:

1

1

1

1

1

1

1

1

1

1

x:

5.59

0.96

6.96

1.23

1.61

2.94

1.96

3.68

1.54

2.59

group:

2

2

2

2

2

2

2

2

2

2

x:

4.5

3.92

10.33

8.23

2.07

4.9

6.84

6.42

3.72

6

group:

3

3

3

3

3

3

3

3

3

3

x:

1.35

1.06

0.74

0.96

1.16

2.08

0.69

0.68

0.84

1.34

group:

4

4

4

4

4

4

4

4

4

4

x:

1.4

1.51

2.49

1.74

1.59

1.36

3

4.81

5.21

5.12

group:

5

5

5

5

5

5

5

5

5

5

kwallis x,by(group)

Test:Equalityofpopulations(Kruskal-WallisTest)

                                            

group_Obs_RankSum

110188.50

210280.50

310420.00

41095.00

510291.00

chi-squared=27.856with4d.f.                                    

probability=0.0001                              

①为各组的秩和值;②为该统计量的c2检验值;③为无效假设检验所对应的p值。

 本例结果表明:5组的中位数有显著的差异。即:5 个不同饲料组的小鼠肝脏中铁的含量有显著差异,说明小鼠肝脏中铁的含量与喂养的饲料有关。

相关资讯

Stata第四章 t检验和单因素方差分析(下)

·单因素方差分析  单因素方差分析又称为OnewayANOVA,用于比较多组样本的均数是否相同,并假定:每组的数据

第四章 t检验和单因素方差分析(上)

本 节STATA  命 令 摘 要 ttest  varname = #                        &

第三章 正态检验与作图命令

本 节STATA  命 令 摘 要 swilk  变 量 名1  变 量 名2  …    变 量 名m graph 变量名  [, bin(#) ] graph  变量名1  变量名2[, 连接设置  曲线上符号设置] ·   &nb

第二章 描述性统计命令与输出结果说明

本 节STATA  命 令 摘 要 by   分组变量:]summarize 变 量 名1  变 量 名2  …    变 量 名m[,detail] ci 变 量 名1  变 量 名2  …    变 量 名m  [, level(#)

Stata 第一章数据输入,存盘和调用文件命令以及数据管理命令

 本 节  STATA 命 令 摘 要[ [1] ]  输 入 数 据 命 令: 1) input    变量名1 变量名2 … 变量名m    ( 各 变 量 数 据 之 间 用 空 格 隔 开, 每 行 一 个 记 录, 最 后 以  end  表 示 结 束) 2)in

第一章 Stata 概貌

  第一章 Stata 概貌 §1.1 Stata的功能、特点和背景 Stata是一个用于分析和管理数据的功能强大又小巧玲珑的实用统计分析软件,由美国计算机资源中心(Computer Resource Center)研制。从1985至1998的十四年时间里,已连续推出1.1,1.2,1.3,1.4,1.5,……及2.0,2.1,3.0,3.1,4.0,5.0,6.0等多个版本,通过