卡方线性趋势检验_双向有序分类资料线性趋势检验的几个问题商榷

中国卫生统计 2011年 4月第 28卷第 2期 双 向有序分类资料线性趋势检验的几个 问题商榷 唐 启义 陈惠君 董 玉恒 国内相关文献¨ 在介绍双向有序列联表资料线 性趋势检验时，均将列联表 Pearson卡方检验 值分 解成线性回归分量 归与偏离线性回归分量 ，即 2总 = 2回归 + 2偏 (1) 据此计算 2偏 = 2 总 一 2 同归 (2) 并认为若 回2归和 均有统计学意义时 ，说 明两分 类变量存在的相关关系不是简单的直线关系 ；只有在 归 有统计学意义，而 无统计 学意义时 ，两分类变 量 仅是线 性关 系 。 文献[4]注意到 ， 归计算 在文献 [1，2]和文 献 [3]中给出的公式不同，且依据文献 [3]得到的 归有 时比 大，亦即 有时小于 0，并分析了它们之 问的 差别与联系。但结论仍认为式(1)成立 ，只是适 用条 件需 要调 整 。 1．关于 归及回归系数 b标准误计算 对双向有序分组资料 ，先对行变量 X、列变量 Y分 别赋予分值 (1，2，3，⋯)，再计算其 回归系数 b和 b的 标准误 se(b)。在一般的线性 回归中 6 = Ixy (3) P( )： (y_f)2／ ： ／ 2／ (4) 一 n一 ～|H 、 然后根据 b及其标准误 ，计算统计量 2 = (b／se(b)) (5) 这和文献 [3]介绍 的方法一致 ，应 用该方法计算 得到的统计量服从 F分布而不是 分布。只有 当样 本量 很大时它才服从 分布㈩。 2．非零相关统计量 文献[5]中提出用非零相关 (nonzero correlation) 统计量来检验双向有序列联表的线性趋势 ，其渐近卡 方统计量 M2=(，z一1)r2 (6) 这里 自由度等于 1；，．为相关系数 ： l浙江大学农、Ik与生物技术学院(310029) 2．杭州市上城 区中医院 3．吉林 大学 白求恩 医学部 (130021) rz= ／( lyy) (7) 在文献 [1，2]中提出线性趋势检验 ，其卡方统计 量为 ： ／ ／ 、2t 归 易。／ txy／txx] ： ． ： ． 牟 ： ， 812 ——‘ ，z。 nr ( ) lxx ‘YY ‘ YY 因此文献 [1，2]中的卡方统计量实 际上为文 献 [5]中的非零 相关统计量 ，只是文献 [1，2]中 r2是乘 以 ，z而文献 [5]中乘 以 ，z一1；故文献 [1，2]中的 S 并 不是回归系数标准误。 3．关于偏离线性回归分量 如上所述 ，目前 国内很多文献，如文献[1，2，3，4] 均认为等式(1)成立 ，亦即认为列联表 Pearson卡方值 ( )可分解 成线性 回归分量 归和偏离线性 回归分 里 2 偏 两个部分 ，从而得 到 的估计值。分析之前 ， 我们先给出 和 归的统计学含义： 列联表中 Pearson卡方检验公式是 r， ， 、2 X。=∑ (9) j e 式 (9)中 是两维列联表的每一个格子的观察频 次 ， 是理论频次，其概率为其行 (x)与列 (y)边际概 率之乘积。这样，我们就通过频数分布列联表得到了相 对频数分布的列联表(或称频率分布列联表)。因此，列 联表 Pearson卡方 )所研究分类变量之间的关联性， 实际上是通过相对频数条件分布的比较进行的。 而线性回归卡方 归，由线性 回归的知识 ，回归系 数 b统计检验 ： =0对 ：卢≠0。因样本 回归系数 b是总体 回归系数 的无偏估计 ，这 时原假设 也可写 为 Ho：E(6)：0。在 成立下 ，t=b／se(b)一t(n一 2)，