1.是；2.否
•

你赞成还是反对这一改革方案?
1.赞成；2.反对
4. 5.
对分类数据的描述和分析通常使用列联表 可使用检验
分类数据与列联表
一.
分类数据
二.
三.
列联表的构造
列联表的分布
列联表的构造
列联表
(contingency table)
1. 2. 3.
4.
5.
6.
由两个以上的变量交叉分类的频数分布表 行变量的类别用 r 表示， ri 表示第 i 个类别 列变量的类别用 c 表示， cj 表示第 j 个类别 每种组合的观察频数用 fij 表示 表中列出了行变量和列变量的所有可能的组合， 所以称为列联表 一个 r 行 c 列的列联表称为 r c 列联表
2 ( f e ) 2 e 3.0319
合计：3.0319
拟合优度检验
分类数据的假设检验
分类数据
比例检验
一个总体 两个以上总体
独立性检验
Z 检验
Z 检验
两个总体
检验
检验
拟合优度检验
(goodness of fit test)
1.
2.
检验多个比例是否相等 检验的步骤
63.35 13.6% 23.4% 36.7% 7.9%
21.4%
71.8% 18.8% 22.0% 28.2% 7.4%
26.2%
—
— 33.6% — — 100%
期望频数的分布
1. 2.
假定行变量和列变量是独立的 一个实际频数 fij 的期望频数 eij ，是总频 数的个数 n 乘以该实际频数 fij 落入第 i 行 和第j列的概率，即
2 2 2 2 ( a e ) ( b e ) ( c e ) ( d e ) 11 12 21 22 2 e11 e12 e21 e22
•
计算检验的统计量
2
i 1 j 1
r
c
( f ij eij ) 2 eij

进行决策

根据显著性水平和自由度(r-1)(c-1)查出临界值2 若22，拒绝H0；若2<2，接受H0
独立性检验
(例题分析)
【例】一种原料来自三个不同的地区，原料质量被分成三个 不同等级。从这批原料中随机抽取 500件进行检验，结果如 下表。检验各地区与原料之间是否存在依赖关系（ 0.05)
因素 X x1 x2 合计
因素 Y y1 y2 合计
a c a+c
b d b+d
a+b c+d n
相关系数
(原理分析)

列联表中每个单元格的期望频数分别为 (a b)( a c) (a c)(c d ) e11 e21 n n (a b)(b d ) (b d )(c d ) e12 e22 n n 将各期望频数代入 的计算公式得
拟合优度检验


H0: 1= 2= 3= 4 H1: 1234 不全相等 = 0.1 df = (2-1)(4-1)= 3 临界值(s):
(例题分析)
统计量:
2
i 1 j 1 r c
( f ij eij ) 2 eij
3.0319
决策:
=0.1 在 = 0.1的水平上不能拒绝H0
立

= 0.05 df = (3-1)(3-1)= 4 临界值(s): =0.05
2
i 1 j 1
r
c
( f ij eij ) 2 eij
19.82
决策:
在 = 0.05的水平上拒绝H0
结论:
0
9.488
19.82

地区和原料等级之间存在依 赖关系
似然比卡方统计量
•

提出假设
H0：1 = 2 = … = j；H1： 全相等
1
, 2 , … , j 不
•

计算检验的统计量 进行决策

2
i 1 j 1
r
c
( f ij eij ) 2 eij
根据显著性水平和自由度(r-1)(c-1)查出临界值2 若2>2，拒绝H0；若2<2，接受H0

1.
如果两个变量之间为完全关联，那么测度值等于+1(正相关) 或-1(负相关)。
如果两个变量无关，那么测度值等于0。
2.
3.
测度值的符号表明关联的方向。大于0的值 (正数)对应于正 关系;小于0的值 (负数)对应于负关系。
两个变量之间的关联越紧密，测度值越大。统计量的绝对值 (不考虑符号)用来衡量关联度的大小。
一分公司 二分公司 三分公司 四分公司
赞成该方案 反对该方案
合计 合计 279 141 420
68 32
100
75 75
120
57 33
90
79 31
110
列联表的分布
观察值的分布
1.
边缘分布
•

行边缘分布
行观察值的合计数的分布 例如，赞成改革方案的共有279人，反对改革方案的141人
•

一分公司 二分公司 三分公司 四分公司
赞成该 方案
实际频数 期望频数
实际频数 期望频数
68 66
32 34
75 80
75 40
57 60
33 30
79 73
31 37
反对该 方案
拟合优度检验
一. 二.
统计量
拟合优度检验

统计量
统计量
1. 2. 3.
用于检验列联表中变量间拟合优度和独立性 用于测定两个分类变量之间的相关程度 计算公式为
地区
甲地区
一级
二级
三级
合计
140
52
64
24
乙地区
丙地区
合计
60
50
162
59
65
188
52
74
150
171
189 500
独立性检验
(例题分析)
1.
提出假设
• • H0：地区与原料等级之间独立 H1：地区与原料等级之间不独立
2.
计算检验的统计量

2 i 1 j 1 r c
( f ij eij ) 2 eij
第二讲 分类数据统计分析
1. 2.
分类数据与列联表 卡方检验
3.
相关的度量
数据的类型与列联分析
数据
数值型数据
离散数据 连续数据
分类数据(顺序数据)
列联分析
分类数据
1.
分类变量的结果表现为类别
• 例如：性别 (男, 女)
2. 3.
各类别用符号或数字代码来测度 使用分类或顺序尺度
•

你吸烟吗?
r1 c1 n n
由于观察频数的总数为n ，所以f11 的期望频数 e11 应为
r1 c1 r1c1 279 100 e11 n 66.43 66 420 n n n
期望频数的分布
(例题分析)
1. 2. 3.
测度22列联表中数据相关程度 对于22 列联表， 系数的值在0～1之间 相关系数计算公式为
n
2
2 r c （f ij eij ） 2 式中： i 1 j 1 e ij n为实际频数的总个数，即样本容量
相关系数
(原理分析)

一个简化的 22 列联表
列联表的结构
(2 2 列联表)
列(cj) 行 (ri) j =1 列( cj ) 合计 j =2
i =1 i =2
合计
f11 f21
f11+ f21
f12 f22
f12+ f22
f11+ f12 f21+ f22
n
列联表的结构
(r c 列联表的一般表示)
列(cj) 行(ri)
列(cj)
j =1
结论:
可以认为四个分公司对改革方 案的赞成比例是一致的
0
3.0319 6.215
2
四格表资料使用条件
（1）n≥40且T≥5 用Pearson χ2检验
（2）n≥40且1≤T<5连续性校正χ2检验
（3）n<40或T<1 用Fisher确切概率法
29
行×列表资料使用条件

要求每个格子中的理论频数T均大于5或1<T<5 的格子数不超过总格子数的1/5；
(例题分析)
实际频数 期望频数
(fij)
68 75 57 79 32 45 33 31
(eij)
66 80 60 73 34 40 30 37
fij - eij
2 -5 -3 6 -2 5 3 -6
(fij- eij)2
eij
0.0606 0.3125 0.1500 0.4932 0.1176 0.6250 0.3000 0.9730
19.82
3. 根据显著性水平 ＝ 0.05 和自由度 (3-1)(3-1)=4 查 出 相 应 的 临 界 值 2=9.488 。 由 于 2=19.82>2＝9.448，拒绝H0
独立性检验
(例题分析)

H0:地区与原料等级之间独立 统计量: H1:地区与原料等级之间不独
2
i 1 j 1 r c
( f ij eij ) 2 eij
其自由度为(r 1)(c 1) 式中：f ij