数理统计复习笔记七——列联表的独立性检验

这篇博客介绍了二维列联表的结构及其在统计学中的应用,特别是在检验两个分类指标A和B是否独立时。通过建立列联表,利用χ2拟合优度检验,当原假设H0:A与B独立时,计算统计量χ2,并根据渐进分布进行显著性检验。在H0成立且样本量足够大时,χ2统计量遵循χ2分布,从而确定拒绝域,完成独立性检验。

1. 一般的二维列联表

B1B_1B1B2B_2B2⋯\cdotsBsB_sBs合计
A1A_1A1n11n_{11}n11n12n_{12}n12⋯\cdotsn1sn_{1s}n1sn1⋅n_{1\cdot}n1
A2A_2A2n21n_{21}n21n22n_{22}n22⋯\cdotsn2sn_{2s}n2sn2⋅n_{2\cdot}n2
⋮\vdots⋮\vdots⋮\vdots⋯\cdots⋮\vdots⋮\vdots
ArA_rArnr1n_{r1}nr1nr2n_{r2}nr2⋯\cdotsnrsn_{rs}nrsnr⋅n_{r\cdot}nr
合计n⋅1n_{\cdot1}n1n⋅2n_{\cdot2}n2⋯\cdotsn⋯n_{\cdots}nnnn

其中,ni⋅=∑j=1snijn_{i\cdot}=\sum\limits_{j=1}^sn_{ij}ni=j=1snijn⋅j=∑i=1rnijn_{\cdot j}=\sum\limits_{i=1}^rn_{ij}nj=i=1rnij,此时指标A,BA, BA,B分别有r,sr, sr,s个水平,且以nijn_{ij}nij表示在nnn个样本中属于Ai∩BjA_i\cap B_jAiBj的样本个数。

2. 假设

考虑两个指标之间是否独立,即H0:指标A与B独立或指标A与B没有关系(1)H_0:指标A与B独立或指标A与B没有关系\tag1H0:ABAB(1)
如记pij=P{X∈Ai∩Bj}p_{ij}=P\{X\in A_i\cap B_j\}pij=P{XAiBj},则这nnn个样本可以看成来自多项分布XXX的样本。

再记pi⋅=P{X∈Ai},i=1,⋯ ,rp_{i\cdot}=P\{X\in A_i\}, i=1,\cdots,rpi=P{XAi},i=1,,rp⋅j=P{X∈Bj},j=1,⋯ ,sp_{\cdot j}=P\{X\in B_j\}, j=1,\cdots,spj=P{XBj},j=1,,s,则有pi⋅=∑j=1spijp_{i\cdot}=\sum\limits_{j=1}^sp_{ij}pi=j=1spijp⋅j=∑i=1rpijp_{\cdot j}=\sum\limits_{i=1}^rp_{ij}pj=i=1rpij,且有如下约束∑i=1rpi⋅=∑j=1sp⋅j=1(2)\sum\limits_{i=1}^rp_{i\cdot}=\sum\limits_{j=1}^sp_{\cdot j}=1\tag2i=1rpi=j=1spj=1(2)

H0H_0H0成立时,应该有pij=pi⋅p⋅jp_{ij}=p_{i\cdot}p_{\cdot j}pij=pipj,于是假设(2)(2)(2)等价于H0:pij=pi⋅p⋅j(3)H_0:p_{ij}=p_{i\cdot}p_{\cdot j}\tag3H0:pij=pipj(3)

3.检验

由于我们可以把上述列联表数据看作时多项分布的样本,故可以用χ2\chi^2χ2拟合优度检验对其独立性假设(3)(3)(3)进行显著性检验。

不过由于pi⋅p_{i\cdot}pip⋅jp_{\cdot j}pj均未知,且有约束(2)(2)(2),故当H0H_0H0成立时,共有r+s−2r+s-2r+s2个未知参数,此时,其未知参数的极大似然估计为p^i⋅=ni⋅n,p^⋅j=n⋅jn\hat p_{i\cdot}=\frac{n_{i\cdot}}{n}, \hat p_{\cdot j}=\frac{n_{\cdot j}}{n}p^i=nni,p^j=nnj

于是有统计量为χ2=n∑i=1r∑j=1s(nij−ni⋅n⋅jn)2ni⋅n⋅j(4)\chi^2=n\sum_{i=1}^r\sum_{j=1}^s\frac{(n_{ij}-\frac{n_{i\cdot}n_{\cdot j}}{n})^2}{n_{i\cdot}n_{\cdot j}}\tag4χ2=ni=1rj=1sninj(nijnninj)2(4)

且当H0H_0H0成立及n→∞n\to\inftyn时,有χ2→χ2((r−1)(s−1))\chi^2\to\chi^2((r-1)(s-1))χ2χ2((r1)(s1))
于是,拒绝域为W={χ2≥χα2((r−1)(s−1))}(5)W=\{\chi^2\ge\chi^2_\alpha((r-1)(s-1))\}\tag5W={χ2χα2((r1)(s1))}(5)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值