Q:如何在SPSS中实现四格表的卡方检验?A:在多数统计软件中,四格表(和行*列表)的数据格式均为行变量、列变量和频数变量。
如下面这个四格表的数据及相应格式如下:分析时首先选择菜单Data->Weight Cases,将频数变量选入Frequency格中,按OK确认。
此时系统就会以频数表的形式来读取所输入的数据,既记录数应为34+12+23+26=95例,而不是4例。
然后选择菜单Analyze->Descriptive Statistics->Cross Tables,将行、列变量分别选入相应的Row、Column格中,再按下方的Statistics钮,选中左上角的Chi-square复选框,按Continue钮,最后按OK即可。
Spss电脑实验-第三节(1) 您要打印的文件是:Spss电脑实验-第三节(1) 打印本文Spss电脑实验-第三节(1)作者:佚名转贴自:本站原创点击数:74第三节不同对象有关指标发生率(百分比)间的比较Ⅰ.两种对象率(百分比)间的比较—四格表χ2 检验χ2 检验(chi-square test),χ为希腊文字母,读作 [kai](卡);chi-square读作“卡方”;χ2检验即“卡方检验”。
它是一种用途较广的假设检验方法,是分析计率或百分比及某些等级资料常用的方法,可分析两个或两个以上率(或百分比)差别的显著性。
1. 两个样本率(百分比)比较—一般四格表的χ2检验四格表的χ2检验用于分析两组或两组以上率(或百分比)差别的显著性。
χ2 = ∑[(∣A - T∣)2 / T ].....................................(3-1)式中 A 为四格表各格子中的实际数,T 为理论数。
χ2 =(ad-bc)2 n / [(a+b)(c+d)(a+c)(b+d)].........................(3-2)例如: 某医院用两种疗法(drug)治疗慢性肾炎病人,结果(effect)西药组有效率为 79.7%(63/79);中药组有效率为87.0(47/54),见表 3-1-a。
问中药治疗的有效率是否较高?(引自何清波等主编,医学统计学及其软件包,2002年,P.183)。
表3-1-a 两种疗法治疗慢性肾炎病人的结果治疗用药治疗例数有效例数有效率(%)西药 79 63 79.7中药54 4787.0--------------------------------------------------------------------------------合计 133 110 82.7 表 3-1是学术论文或总结报告中的常用形式;而在进行χ2检验时,不能用这个表的形式,必须列出象表 3-1-b的表格形式,其中要从“治疗例数合计”中分出“无效例数”,这才是所谓著名的“四格表”(fourfold table,或2×2表,即 a =63、b = 16、c = 47、d = 7 这 4 格),而不是表 3-1-a 中的 79、63、54、47 这 4 格!表3-1-b 两种疗法治疗慢性肾炎病人的结果治疗结果(effect)治疗用药---------------------------------(drug)有效例数 (1) 无效例数 (2)治疗例数合计有效率(%)西药 (1) 63 (a) 16 (b) 79 (a+b) 79.7中药 (2) 47 (c) 7 (d) 54 (c+d)87.04-----------------------------------------------------------------------------------合计 110 (a+c) 23 (b+d) 133 (a+b+c+d=n) 82.7运算所用程序如下:对照表3-1-b按行、按列输入数据到 BEGIN DATA 与 END DATA 之间。
命令中用了一句 VALUELABELS,是用来标示疗法(drug)与治疗效果(effect),使之输出结果更易阅读,对计算结果并无影响。
所用程序文件名是Chi2-test.sps 的例 *1。
程序中 WEIGHT BY count 命令的产生方法是:在数据编辑窗 Data Editor 中录入数据→点击 Data 主菜单→点击Weight Cases 主菜单→改为勾选Weight cases by →将 count选入 Frequency Variable 框中→ Paste →出现WEIGHTBY count 命令。
程序中 CROSSTABS 命令的产生方法是:SPSS 程序编辑窗主菜单 Analyze →选 Descriptive Statistics 菜单→选 Crosstabs →将 drug 选入右边的 Rows 框中(表示“药物 drug”在“行 Row”的方向上)→将 effect(疗效)选入右边的Column(s) 框中→点击下方的 Statistics 钮,出现另一 Statistics窗→勾选Chi-Square(χ2 值)→点击下方的 Cells 钮,出现另一 Cell Display 窗→勾选 Observed(观察值)、Expected (期望值)、Row (行方向上的百分率),即两种药物的有效率→ Continue → Paste,即出现“CROSSTABS”命令。
* Some Example for Chi-Square Test; Filename:Chi2-test.sps.*--------------------------------------------------------------------------.*1.Chi-square test for fourfold table (2 BY 2 table):.DATA LIST FREE /drug effect count.BEGIN DATA.1 1 63 1 2 16 2 1 47 2 2 7ENDDATA.VALUE LABELS drug 1 '西药' 2 '中药'/effect 1 '有效' 2 '无效'.WEIGHTBY count.CROSSTABS TABLES=drug BY effect /CELLS=COUNT ROW EXPECTED/STATISTICS=CHISQ.*-------------------------------------------------------------------------------.Chi-Square TestsValue df Asymp. Sig. (2-sided) Exact Sig. (2-sided) Exact Sig.(1-sided)Pearson Chi-Square 1.192 1 .275Continuity Correction .737 1 .391Likelihood Ratio 1.226 1 .268Fisher's Exact Test .353 .196Linear-by-Linear Association 1.183 1 .277N of Valid Cases 133a Computed only for a 2x2 tableb 0 cells (.0%) have expected count less than 5. The minimumexpected count is 9.34.程序中的“CROSSTABS”一句用于计算χ2值。
其中的关键词即用于计算各格子的理论数 T。
主要运算结果:χ2 (PearsonChi-Square)= 1.192, P = 0.275。
可见两种药物的有效率差异无统计学意义,因而不能认为中药治疗的有效率较高。
2. 两种方法效果的比较—配对计数资料的χ2检验有些调研工作要进行两种方法效果的比较,可用配对计数资料的χ2检验。
目的是通过单一样本数据,推断两种处理的效果差异有无显著性。
例如:采集了205份标本,将每份标本各分别接种于甲种与乙种培养基,观察其生长情况,结果见表3-2,问哪种培养基效果效果较好?(引自杨树勤主编.中国医学百科全书—医学统计学,P.93)。
从表中可见:甲种培养基生长率为 34.1% (70/205),乙种为 17.6% (36/205)。
注意:表 3-2 与表 3-1-b 对比,形式上类似,也象“四格表”,但意义不同!因为其中a、b、c、d的每份标本都接种了两种培养基,所以就数据而言加了倍。
因而所用的χ2检验公式也不同,要用下式计算χ2值:χ2 = (∣b -c∣-1)2 / (b+c) .......................................(3-3)表3-2 甲乙两种培养基培养结果的比较乙种培养基(B)甲种培养基--------------------------- 合计(A)生长不生长生长36 34 70不生长0135135-------------------------------------------------------------合计36 169 205所用程序文件名为Chi2-test.sps 的例 *2:*---------------------------------------------------------------------------------.*2.Yang Shu-Qing: CHN Encyclop. Med. Stat. Chi-square for pairedsample:.DATA LIST FREE /a b c d.BEGIN DATA.36 34 0 135ENDPUTE Chi2=(ABS(b-c)-1)**2/(b+c).LIST VARIABLES=Chi2/CASES=FROM 1 TO 1.*So getting Chi2=PUTE Chi2=PUTE df=PUTE P = PDF.CHISQ(Chi2,df)*2.FORMATS P(F6.4).LISTVARIABLES=Chi2 P /CASES=FROM 1 TO 1.*Or using non-parameter method, McNemar:.DATA LIST FREE /a b f.BEGIN DATA.1 1 361 0 340 1 00 0 135END DATA.VALUE LABELS a 1 '+' 0 '0' /b 1 '+' 0 '0'. WEIGHT BYf.NPAR TEST /MCNEMAR= a WITH b(PAIRED).*----------------------------------------------------------------------------------.在进行配对χ2 检验时,用一句“COMPUTE Chi2 = (ABS(b-c)-1)**2/(b+c).”此命令按公式(3-3)计算χ2 值,结果χ2 = 32.03。