卡方检验

出自 MBA智库百科(http://wiki.mbalib.com/)

卡方检验(Chi-square test/Chi-Square Goodness-of-Fit Test)

目录

什么是卡方检验

  卡方检验是一种用途很广的计数资料的假设检验方法。它属于非参数检验的范畴,主要是比较两个及两个以上样本率( 构成比)以及两个分类变量的关联性分析。其根本思想就是在于比较理论频数和实际频数的吻合程度或拟合优度问题。

  它在分类资料统计推断中的应用,包括:两个率或两个构成比比较的卡方检验;多个率或多个构成比比较的卡方检验以及分类资料的相关分析等。

卡方检验的基本原理[1]

卡方检验的基本思想

  卡方检验是以χ2分布为基础的一种常用假设检验方法,它的无效假设H0是:观察频数与期望频数没有差别。

  该检验的基本思想是:首先假设H0成立,基于此前提计算出χ2值,它表示观察值与理论值之间的偏离程度。根据χ2分布及自由度可以确定在H0假设成立的情况下获得当前统计量及更极端情况的概率P。如果P值很小,说明观察值与理论值偏离程度太大,应当拒绝无效假设,表示比较资料之间有显著差异;否则就不能拒绝无效假设,尚不能认为样本所代表的实际情况和理论假设有差别。

卡方值的计算与意义

  χ2值表示观察值与理论值之问的偏离程度。计算这种偏离程度的基本思路如下。

  (1)设A代表某个类别的观察频数,E代表基于H0计算出的期望频数,A与E之差称为残差。

  (2)显然,残差可以表示某一个类别观察值和理论值的偏离程度,但如果将残差简单相加以表示各类别观察频数与期望频数的差别,则有一定的不足之处。因为残差有正有负,相加后会彼此抵消,总和仍然为0,为此可以将残差平方后求和。

  (3)另一方面,残差大小是一个相对的概念,相对于期望频数为10时,期望频数为20的残差非常大,但相对于期望频数为1 000时20的残差就很小了。考虑到这一点,人们又将残差平方除以期望频数再求和,以估计观察频数与期望频数的差别。

  进行上述操作之后,就得到了常用的χ2统计量,由于它最初是由英国统计学家Karl Pearson在1900年首次提出的,因此也称之为Pearson χ2,其计算公式为
\chi^2=\sum \frac{(A-E)^2}{E}=\sum_{i=1}^k \frac{(A_i-E_i)^2}{E_i}=\sum_{i=1}^k \frac{(A_i-np_i)^2}{np_i}  (i=1,2,3,…,k)

  其中,Ai为i水平的观察频数,Ei为i水平的期望频数,n为总频数,pi为i水平的期望频率。i水平的期望频数Ti等于总频数n×i水平的期望概率pi,k为单元格数。当n比较大时,χ2统计量近似服从k-1(计算Ei时用到的参数个数)个自由度的卡方分布

  作为学术界的领袖,Pearson先生当初发表在《哲学杂志》上的χ2论文题目为:On the criterion that a given system of deviations from the probable in the case of a correlated system of variables is such that it can be reasonably supposed to have arisen from random sampling.

  由卡方的计算公式可知,当观察频数与期望频数完全一致时,χ2值为0;观察频数与期望频数越接近,两者之间的差异越小,χ2值越小;反之,观察频数与期望频数差别越大,两者之间的差异越大,χ2值越大。换言之,大的χ2值表明观察频数远离期望频数,即表明远离假设。小的χ2值表明观察频数接近期望频数,接近假设。因此,χ2是观察频数与期望频数之间距离的一种度量指标,也是假设成立与否的度量指标。如果χ2值“小”,研究者就倾向于不拒绝H0;如果χ2值大,就倾向于拒绝H0。至于χ2在每个具体研究中究竟要大到什么程度才能拒绝H0,则要借助于卡方分布求出所对应的P值来确定。

卡方检验的样本量要求

  卡方分布本身是连续型分布,但是在分类资料的统计分析中,显然频数只能以整数形式出现,因此计算出的统计量是非连续的。只有当样本量比较充足时,才可以忽略两者问的差异,否则将可能导致较大的偏差具体而言,一般认为对于卡方检验中的每一个单元格,要求其最小期望频数均大于1,且至少有4/5的单元格期望频数大于5,此时使用卡方分布计算出的概率值才是准确的。如果数据不符合要求,可以采用确切概率法进行概率的计算。

卡方检验的类型[2]

  1、四格表资料的卡方检验

  四格表资料的卡方检验用于进行两个率或两个构成比的比较。

  1)专用公式:

  若四格表资料四个格子的频数分别为a,b,c,d,则四格表资料卡方检验的卡方值=(ad-bc)2 \times n/(a+b)(c+d)(a+c)(b+d),自由度v=(行数-1)(列数-1)

  2)应用条件:

  要求样本含量应大于40且每个格子中的理论频数不应小于5。当样本含量大于40但理论频数有小于5的情况时卡方值需要校正,当样本含量小于40时只能用确切概率法计算概率。

  2、行×列表资料的卡方检验

  行×列表资料的卡方检验用于多个率或多个构成比的比较。

  1)专用公式:

  r行c列表资料卡方检验的卡方值=n[(A_{11}/n_1n_1+A_{12}/n_1n_2+\cdots+A_{rc}/n_rn_c)-1]

  2)应用条件:

  要求每个格子中的理论频数T均大于5或1<T<5的格子数不超过总格子数的1/5。当有T<1或1<T<5的格子较多时,可采用并行并列、删行删列、增大样本含量的办法使其符合行×列表资料卡方检验的应用条件。而多个率的两两比较可采用行×列表分割的办法。

  3、列联表资料的卡方检验

  同一组对象,观察每一个个体对两种分类方法的表现,结果构成双向交叉排列的统计表就是列联表。

  1)R×C 列联表的卡方检验:

  R×C 列联表的卡方检验用于R×C列联表的相关分析,卡方值的计算和检验过程与行×列表资料的卡方检验相同。

  2)2×2列联表的卡方检验:

  2×2列联表的卡方检验又称配对记数资料或配对四格表资料的卡方检验,根据卡方值计算公式的不同,可以达到不同的目的。当用一般四格表的卡方检验计算时,卡方值=(ad-bc)2n/(a+b)(c+d)(a+c)(b+d),此时用于进行配对四格表的相关分析,如考察两种检验方法的结果有无关系;当卡方值=( | bc | − 1)2/(b+c)时,此时卡方检验用来进行四格表的差异检验,如考察两种检验方法的检出率有无差别。

  列联表卡方检验应用中的注意事项同R×C表的卡方检验相同。

卡方检验的用途[1]

  卡方检验最常见的用途就是考察某无序分类变量各水平在两组或多组间的分布是否一致实际上,除了这个用途之外.膏方检验还有更广泛的应用。具体而言,其用途主要包括以下几个方面:

  (1)检验某个连续变量的分布是否与某种理论分布相一致。如是否符合正态分布、是否服从均匀分布、是否服从Poisson分布等。

  (2)检验某个分类变量各类的出现概率是否等于指定概率。如在36选7的彩票抽奖中,每个数字出现的概率是否各为1/36;掷硬币时,正反两面出现的概率是否均为0.5。

  (3)检验某两个分类变量是否相互独立。如吸烟(二分类变量:是、否)是否与呼吸道疾病(二分类变量:是、否)有关;产品原料种类(多分类变量)是否与产品合格(二分类变量)有关。

  (4)检验控制某种或某几种分类因素的作用以后,另两个分类变量是否相互独立。如在上例中,控制性别、年龄因素影响以后,吸烟是否和呼吸道疾病有关;控制产品加工工艺的影响后,产品原料类别是否与产品合格有关。

  (5)检验某两种方法的结果是否一致。如采用两种诊断方法对同一批人进行诊断,其诊断结果是否一致;采用两种方法对客户进行价值类别预测,预测结果是否一致。

卡方检验的应用条件[3]

  适用于四格表应用条件:

  1)随机样本数据。两个独立样本比较可以分以下3种情况:

  (1)所有的理论数T≥5并且总样本量n≥40,用Pearson卡方进行检验。

  (2)如果理论数T<5但T≥1,并且1≥40,用连续性校正的卡方进行检验。

  (3)如果有理论数T<1或n<40,则用Fisher’s检验。

  2)卡方检验的理论频数不能太小。

  R×C表卡方检验应用条件:

  (1)R×C表中理论数小于5的格子不能超过1/5;

  (2)不能有小于1的理论数。如果实验中有不符合R×C表的卡方检验,可以通过增加样本数、列合并来实现。

卡方检验应用实例[4]

  1.应用实例——适合度检验

  实际执行多项式试验而得到的观察次数,与虚无假设的期望次数相比较,称为卡方适度检验,即在于检验二者接近的程度,利用样本数据以检验总体分布是否为某一特定分布的统计方法。这里以掷骰子为例介绍适度检验的方法。

  【例1】

  (1)假设掷一骰子120次,各点数共出现次数为a,b为各点数出现的期望值120×1/6=20,建立工作表文件,如图1所示。

  (2)设置零假设H0:观察分布等于期望分布。

  (3)计算卡方检验统计量,如图2所示。

  D2=(B2-C2)^2/C2

  D8=SUM(D2:D7)

Image:适合度检验.jpg

  (4)确定自由度,(6-1)×(2-1)=5;选择显著水平α=0.05。

  (5)利用Excel提供的CHIINV函数求临界值,在D9单元格中键入“=CHIINV(0.05,5)”按回车键,得临界值11.07。

  (6)比较临界值和统计量,11.07>2.3,即临界值大于统计量,故差异不显著,接受H0

  2.应用实例2——独立性检验

  卡方独立性检验是用来检验两个属性间是否独立。一个变量作为行,另一个变量作为列。下面一例便是介绍卡方独立性检验的方法。

  【例2】某机构欲了解现在性别与收入是否有关,他们随机抽样500人,询问对此的看法,结果分为“有关、无关、不好说,,三种答案,图3中县调查得到的数据。 Image:图3 调查数据表.jpg

  下面是利用Excel解决此问题的步骤。

  (1)零假设H0:性别与收入无关。

  (2)确定自由度为(3-1)×(2-1)=2,选择显著水平α=0.05。

  (3)求解男女对收入与性别相关不同看法的期望次数,这里采用所在行列的合计值的乘机除以总计值来计算每一个期望值,如图4所示,在单元格B9中键入“=B5*E3/E5”,同理(第一个等于号理解为在单元格中键入):

  • B10=“=B5*E4/E5”,
  • C9=“=C5*E3/E5”,
  • C10=“=C5*E4/E5”,
  • D9=“=D5*E3/E5”,
  • D10=“=D5*E4/E5”。

  (4)利用卡方统计量计算公式计算统计量,在单元格B15中键入“=(B3-B9)^2/B9”,其余单元格依次类推,结果如图5所示。 Image:独立性检验.jpg

  \Chi^2=\sum \frac{(f_0 - f_e)^2}{f_e}=21.4675

  (5)利用Excel提供的CHIINV函数计算显著水平为0.05,自由度为2卡方分布的临界值,在Excel单元格中键入“=CHIINV(0.05,2)”按回车键,得临界值为5.9915。

  (6)比较统计量度和临界值,统计量21.4675大于临界值5.9915,故拒绝零假设。

  3.应用实例3——统一性检验

  检验两个或两个以上总体的某一特性分布,也就是各“类别”的比例是否统一或相近,一般称为卡方统一性检验或者卡方同质性检验。下面一例便是利用卡方统一性检验的例子。

  【倒3】某咨询公司想了解南京和北京的市民对最低生活保障的满意程度是否相同。他们从南京抽出600居民,北京抽取600居民,每个居民对满意程度(非常满意、满意、不满意、非常不满意)任选一种,且只能选一种。将统计结果键入Excel工作表中,如图6所示。 Image:图6 抽查统计结果.jpg

  下面是利用Excel解决此问题的步骤。

  (1)零假设H0:南京和北京居民对最低生活保障满意程度的比例相同。

  (2)确定自由度为(4-1)×(2-1)=3,选择显著水平α=0.05。

  (3)求解卡方检验的l临界值,在Excel单元格中键入“=CHIINV(0.05,3)”,按回车键得临界值为7.81。

  (4)计算北京和南京不同满意程度的期望值,在单元格B11和C11中分别键入“=$B$7*D3/$D$7”和“=$C$7*D3/$D$7”,选中B11:C11,按住C11右下角填充控制点,填充至C14。

  (5)计算卡方统计量,在单元格B19中键入“=(B3-B11)^2/B11”,其余单元格依次类推,结果如图7所示。 Image:图7 期望值和统计量计算.jpg

  \Chi^2=\sum \frac{(f_0 - f_e)^2}{f_e} = 1.3875

  (6)比较统计量和临界值,统计量1.3875小于临界值7.81,故接受零假设。

参考文献

  1. 1.0 1.1 张文彤,邝春伟编著.第16章 无序分类变量的统计推断——卡方检验 SPSS统计分析基础教程 第2版.高等教育出版社,2011.11.
  2. 刘铁民主编.第十章 安全生产统计分析 安全生产管理知识 2006版.中国大百科全书出版社,2006年05月第1版.
  3. 卡方检验的应用条件[J].中国医药科学,2013,(第2期).
  4. 孙志刚,杨聪.第八章 假设检验 Excel在经济与数理统计中的应用.中国电力出版社,2004年01月第1版.
本条目对我有帮助151
分享到:
  如果您认为本条目还有待完善,需要补充新内容或修改错误内容,请编辑条目

评论(共19条)

提示:评论内容为网友针对条目"卡方检验"展开的讨论,与本站观点立场无关。
202.114.136.* 在 2008年8月28日 21:02 发表

很好,谢谢。但请问一下,什么是pearson chi-squsre

回复评论
Hnoju (Talk | 贡献) 在 2008年8月29日 13:48 发表

202.114.136.* 在 2008年8月28日 21:02 发表

很好,谢谢。但请问一下,什么是pearson chi-squsre

您指的pearson chi-squsre,不知是否为Pearson Chi-Square,二者有一单词不一样,Pearson Chi-Square(皮尔森卡方检验)即常用的卡方检验

回复评论
211.100.210.* 在 2009年1月13日 14:20 发表

谢谢

回复评论
118.250.2.* 在 2009年10月21日 23:05 发表

请问 卡方检验用于2行c列表的计算时 若不用其专用公式计算卡方值 ,想用理论频数计算时 怎么算

回复评论
58.254.92.* 在 2010年6月2日 13:17 发表

那个四格表不知怎么看的

回复评论
117.27.27.* 在 2010年9月8日 15:21 发表

谢谢

回复评论
116.253.106.* 在 2010年11月7日 23:57 发表

四个表这样看不出来,要经过计算差值、平方值等最后计算卡方值得出结果!

回复评论
114.83.229.* 在 2011年6月25日 14:43 发表

孤烽到似家将热热健康了gre

回复评论
114.83.229.* 在 2011年6月25日 14:45 发表

表1:暗示疗法对癌痛患者止痛剂效果 组别 有效 显效 合计 显效率(%) 实验组 5 39 44 88.64 对照组 13 32 45 71.11 合计 18 71 89 79.78 请帮我计算卡方值,P值有没有意义.谢谢

回复评论
114.83.229.* 在 2011年6月25日 14:47 发表

组别 有效 显效 合计 显效率(%) 实验组 5 39 44 88.64 对照组 13 32 45 71.11 合计 18 71 89 79.78 请帮我计算卡方值,P只是否有意义。万分感谢

回复评论
61.178.18.* 在 2011年10月19日 20:41 发表

谢谢,想问一下,用excel可以卡方检验么?

回复评论
Dan (Talk | 贡献) 在 2011年10月20日 12:16 发表

61.178.18.* 在 2011年10月19日 20:41 发表

谢谢,想问一下,用excel可以卡方检验么?

添加了excel的应用实例,希望对你有帮助~

回复评论
220.248.45.* 在 2012年3月5日 14:10 发表

114.83.229.* 在 2011年6月25日 14:47 发表

组别 有效 显效 合计 显效率(%) 实验组 5 39 44 88.64 对照组 13 32 45 71.11 合计 18 71 89 79.78 请帮我计算卡方值,P只是否有意义。万分感谢

Chi-Square Value: 4.235, Critical Value: 3.841, Accept Ha, There is significant difference between Tested and Control Group. It shows that your treatment is valid.

回复评论
123.182.67.* 在 2012年5月13日 16:29 发表

非常好,看起来简单易懂!谢谢!

回复评论
119.98.146.* 在 2012年5月29日 11:45 发表

第二个例子漏洞百出

回复评论
泡芙小姐 (Talk | 贡献) 在 2012年5月30日 10:02 发表

119.98.146.* 在 2012年5月29日 11:45 发表

第二个例子漏洞百出

已核对参考文献并进行了修改,希望对您有帮助!

  MBA智库百科是可以自由参与的百科,如有发现错误和不足,您也可以参与修改编辑,只要通过网页右上的创建新帐号,创建用户名后即可参与,期待您的加入!~

回复评论
58.246.161.* 在 2012年12月26日 14:29 发表

恩,很好,但是看不懂。

回复评论
180.160.6.* 在 2013年7月21日 20:55 发表

58.246.161.* 在 2012年12月26日 14:29 发表

恩,很好,但是看不懂。

确实看不懂哇,怎么办?还是不知道卡方检验是干嘛的,什么时候用,怎么看!

回复评论
218.4.189.* 在 2014年5月20日 21:00 发表

分类变量的关联系检验。不过卡方检验只能说明统计学上有意义,并且变数略大,还要结合列联系数吧应该是,来对此,sig值一般要小于0.05,才能拒绝原假设。

回复评论

发表评论请文明上网,理性发言并遵守有关规定。

返回顶部