154 中 国 医 药 生 物 技 术 2012 年 4 月 第 7 卷 第 2 期 Chin Med Biotechnol, April 2012, Vol. 7, No. 2 DOI:10.3969/cmba.j.issn.1673-713X.2012.02.015 讲 座 如 何 用 SAS 软 件 正 确 分 析 生 物 医 学 科 研 资 料 XVI. 用 SAS 软 件 2 2 列 联 表 资 料 的 统 计 分 析 关 雪, 胡 良 平, 王 琪 编 者 按 生 物 统 计 学 是 生 物 学 领 域 科 学 研 究 和 实 际 工 作 中 必 不 可 少 的 工 具, 在 分 子 生 物 学 迅 速 发 展 的 今 天, 生 物 统 计 学 更 显 示 出 了 它 的 重 要 性 实 验 设 计 与 数 据 统 计 分 析 是 现 代 生 物 学 的 基 石, 是 生 物 学 研 究 者 检 验 假 说 寻 找 模 式 建 立 生 物 学 理 论 的 有 利 工 具, 也 是 生 物 学 研 究 者 探 索 微 观 和 宏 观 生 物 世 界 的 必 备 基 础 知 识 对 于 每 天 甚 至 是 每 时 每 刻 涌 现 的 大 量 的 以 天 文 数 字 计 量 的 分 子 遗 传 数 据, 必 须 借 助 统 计 学 知 识 加 以 分 析 处 理, 才 能 从 中 获 得 有 意 义 的 信 息 生 物 多 样 性 数 据 分 析 是 开 展 生 物 多 样 性 研 究 的 一 个 重 要 方 面, 数 据 分 析 能 力 的 高 低 极 大 地 影 响 着 我 们 对 各 种 生 态 学 现 象 认 识 的 深 度 和 广 度 现 在, 电 子 计 算 机 的 普 及 使 得 生 物 统 计 分 析 过 程 大 大 简 化, 生 物 统 计 分 析 软 件 包 的 普 及 将 生 物 统 计 学 从 统 计 学 家 的 书 本 里 解 放 了 出 来, 简 化 了 生 物 统 计 分 析 过 程, 使 之 成 为 生 物 学 研 究 者 的 常 用 工 具 本 刊 特 邀 军 事 医 学 科 学 院 生 物 医 学 统 计 学 咨 询 中 心 主 任 胡 良 平 教 授, 以 如 何 用 SAS 软 件 正 确 分 析 生 物 医 学 科 研 资 料 为 题, 撰 写 系 列 统 计 学 讲 座, 希 望 该 系 列 讲 座 能 对 生 物 医 学 科 研 工 作 者 有 所 帮 助 分 析 定 性 资 料 时, 首 先 应 正 确 判 断 资 料 所 对 应 的 列 联 表 类 型 其 次 根 据 不 同 的 分 析 目 的, 并 结 合 统 计 分 析 方 法 的 应 用 条 件, 选 择 合 适 的 分 析 方 法 通 常, 列 联 表 可 分 为 :2 2 表 (4 类 ) R C 表 (5 类 ) 和 高 维 列 联 表 (3 类 ) 以 及 具 有 重 复 测 量 因 素 的 高 维 列 联 表 [1] 本 期 主 要 介 绍 2 2 列 联 表 资 料 统 计 分 析 的 SAS 1 横 断 面 研 究 设 计 2 2 表 资 料 统 计 分 析 的 SAS 在 某 个 时 间 断 面 ( 时 点 或 很 短 时 间 内 ) 进 行 的 调 查 或 试 验 研 究 中, 对 一 组 受 试 对 象 同 时 按 两 个 定 性 变 量 来 划 分, 每 个 定 性 变 量 都 只 有 2 个 水 平, 其 目 的 是 了 解 某 个 时 间 点 的 现 状, 由 此 获 得 的 2 2 表 资 料 称 为 横 断 面 研 究 设 计 的 2 2 列 联 表 资 料 [2] 医 学 中, 通 常 用 于 对 疾 病 或 健 康 状 况 分 布 特 征 的 描 述 对 防 治 措 施 效 果 的 评 价 对 疾 病 与 某 些 因 素 相 互 联 系 的 探 讨 等 例 1 为 探 讨 术 前 3 d 给 予 大 剂 量 麝 香 保 心 丸 对 于 经 皮 冠 状 动 脉 介 入 疗 法 病 人 心 肌 损 伤 的 影 响, 随 机 选 择 133 例 2008 2009 年 拟 行 冠 状 动 脉 造 影 的 不 稳 定 型 心 绞 痛 患 者, 分 组 表 1 两 组 患 者 血 脂 测 定 结 果 ( 人 ) 血 脂 正 常 血 脂 异 常 合 计 对 照 组 39 28 67 试 验 组 34 32 66 合 计 73 60 133 分 为 对 照 组 和 试 验 组, 试 验 组 在 常 规 治 疗 的 基 础 上 加 服 麝 香 保 心 丸 试 验 前 对 两 组 患 者 的 血 脂 情 况 进 行 比 较, 数 据 见 表 1, 试 分 析 两 组 患 者 血 脂 异 常 率 之 间 是 否 存 在 差 别 [1] 分 析 与 SAS : 该 试 验 目 的 是 测 定 试 验 前 两 组 患 者 的 血 脂 异 常 率 是 否 存 在 差 别, 可 视 为 横 断 面 研 究 设 计 2 2 表 资 料 因 总 频 数 133 大 于 40, 且 理 论 频 数 均 大 于 5, 故 可 用 一 般 χ 2 检 验 具 体 的 SAS 程 序 如 下 : data tjfx1 do A=1 to 2 do B=1 to 2 input f@@ output tables A*B/chisq end end ods html close 39 28 34 32 程 序 说 明 : 数 据 步, 建 立 名 为 tjfx1 的 数 据 集, 建 立 数 值 型 变 量 A B f, 分 别 读 入 原 因 变 量 的 水 平 结 果 变 量 的 水 平 和 对 应 单 元 格 中 的 实 际 频 数 过 程 步, 调 用 FREQ 过 程, 指 定 频 数 变 量 为 f, 用 tables A*B 语 句 表 示 生 成 二 维 列 联 表 资 料, 斜 杠 / 后 的 选 项 chisq 要 求 进 行 一 般 χ 2 检 验 SAS 输 出 结 果 与 结 果 解 释 如 下 : 作 者 单 位 :100850 北 京, 军 事 医 学 科 学 院 生 物 医 学 统 计 学 咨 询 中 心 通 讯 作 者 : 胡 良 平,Email:lphu812@sina.com
中 国 医 药 生 物 技 术 2012 年 4 月 第 7 卷 第 2 期 Chin Med Biotechnol, April 2012, Vol. 7, No. 2 155 A*B 表 的 统 计 量 卡 方 1 0.6016 0.4380 似 然 比 卡 方 1 0.6021 0.4378 连 续 校 正 卡 方 1 0.3617 0.5476 Mantel-Haenszel 卡 方 1 0.5971 0.4397 Phi 系 数 0.0673 列 联 系 数 0.0671 Cramer 的 V 0.0673 单 元 格 (1,1) 频 数 (F) 39 左 侧 Pr <= F 0.8289 右 侧 Pr >= F 0.2739 计 算 相 对 危 险 度 RR, 并 进 一 步 对 总 体 RR 值 是 否 等 于 1 进 行 假 设 检 验 具 体 的 SAS 程 序 如 下 : data tjfx2 do a=1 to 2 do b=1 to 2 input f @@ tables a*b /chisq cmh output end end ods html close 3 9 10 1 表 概 率 (P) 0.1028 双 侧 Pr <= P 0.4879 样 本 大 小 = 133 统 计 和 专 业 结 论 : 由 于 χ 2 = 0.6016,P = 0.4380 > 0.05, 表 明 两 组 患 者 的 血 脂 异 常 率 之 间 的 差 别 无 统 计 学 意 义, 即 两 组 患 者 的 血 脂 异 常 情 况 相 同, 差 别 仅 是 抽 样 误 差 所 致 对 于 横 断 面 研 究 设 计 的 2 2 表 资 料 而 言, 当 总 n > 40 且 无 单 元 格 的 理 论 频 数 T < 5 时, 可 选 用 一 般 χ 2 检 验 当 总 n > 40 但 至 少 有 一 个 单 元 格 的 理 论 频 数 1 < T 5 时, 可 选 用 校 正 的 χ 2 检 验 当 总 频 数 n < 40 或 至 少 有 一 个 单 元 格 的 理 论 频 数 T < 1 时, 必 须 选 用 Fisher 精 确 检 验 2 队 列 研 究 设 计 2 2 表 资 料 统 计 分 析 的 SAS 所 谓 队 列 研 究 设 计 的 2 2 表 资 料, 就 是 依 据 专 业 知 识, 通 过 对 不 同 暴 露 水 平 的 对 象 进 行 追 踪 观 察, 确 定 其 疾 病 发 生 情 况, 从 而 分 析 暴 露 因 素 与 疾 病 发 生 之 间 的 因 果 关 系 它 是 由 因 溯 果 的 前 瞻 性 研 究 设 计, 是 论 证 疾 病 因 果 关 系 的 一 个 重 要 研 究 方 法 例 2 为 探 讨 两 种 药 物 治 疗 胃 炎 的 预 后, 将 胃 炎 患 者 随 机 分 为 两 组, 一 组 患 者 服 用 A 药, 另 一 组 患 者 服 用 B 药, 治 愈 后 追 踪 观 察 一 段 时 间, 观 察 胃 炎 复 发 的 情 况, 资 料 见 表 2 试 比 较 两 种 药 物 复 发 率 之 间 的 差 异 有 无 统 计 学 意 义 药 物 种 类 表 2 不 同 药 物 治 疗 胃 炎 的 复 发 情 况 复 发 未 复 发 合 计 A 3 9 12 B 10 1 11 合 计 13 10 23 分 析 与 SAS : 表 2 资 料 为 来 自 队 列 研 究 设 计 的 2 2 列 联 表 资 料, 可 先 按 横 断 面 研 究 设 计 2 2 列 联 表 资 料 的 检 验 方 法 考 察 两 组 频 数 分 布 是 否 相 同, 如 果 不 相 同, 再 程 序 说 明 : 建 立 名 为 tjfx2 的 数 据 集, 建 立 数 值 型 变 量 a b f, 分 别 读 入 行 号 列 号 每 格 实 际 频 数 过 程 步, 调 用 FREQ 过 程, 指 定 频 数 变 量 为 f, 用 tables a*b 语 句 表 示 生 成 二 维 列 联 表 资 料, 斜 杠 / 后 的 选 项 chisq 表 示 要 求 进 行 一 般 χ 2 检 验, 选 项 cmh 要 求 输 出 Cochran- Mantel-Haenszel 统 计 量, 当 列 联 表 为 2 2 表 时, 它 可 以 输 出 OR( 优 势 比 ) 和 RR( 相 对 危 险 度 ) 估 计 值 及 置 信 区 间 SAS 输 出 结 果 与 结 果 解 释 : a*b 表 的 统 计 量 卡 方 1 10.1449 0.0014 似 然 比 卡 方 1 11.2943 0.0008 连 续 校 正 卡 方 1 7.6402 0.0057 Mantel-Haenszel 卡 方 1 9.7038 0.0018 Phi 系 数 0.6641 列 联 系 数 0.5532 Cramer 的 V 0.6641 WARNING: 25% 的 单 元 格 的 期 望 计 数 比 5 小 卡 方 可 能 不 是 有 效 检 验 单 元 格 (1,1) 频 数 (F) 3 左 侧 Pr <= F 0.0022 右 侧 Pr >= F 0.9999 表 概 率 (P) 0.0021 双 侧 Pr <= P 0.0028 以 上 给 出 了 若 干 χ 2 检 验 和 统 计 量 a*b 的 汇 总 统 计 量 Cochran-Mantel-Haenszel 统 计 量 ( 基 于 表 得 分 ) 统 计 量 对 立 假 设 自 由 度 值 概 率 1 非 零 相 关 1 9.7038 0.0018 2 行 均 值 得 分 差 值 1 9.7038 0.0018 3 一 般 关 联 1 9.7038 0.0018
156 中 国 医 药 生 物 技 术 2012 年 4 月 第 7 卷 第 2 期 Chin Med Biotechnol, April 2012, Vol. 7, No. 2 普 通 相 对 风 险 的 估 计 值 ( 行 1/ 行 2) 研 究 类 型 方 法 值 95% 置 信 限 案 例 对 照 Mantel-Haenszel 0.0333 0.0029 0.3808 ( 优 比 ) Logit 0.0333 0.0029 0.3808 Cohort Mantel-Haenszel 0.2750 0.1014 0.7458 ( 第 1 列 风 险 ) Logit 0.2750 0.1014 0.7458 Cohort Mantel-Haenszel 8.2500 1.2375 54.9982 ( 第 2 列 风 险 ) Logit 8.2500 1.2375 54.9982 以 上 第 一 部 分 是 Cochran-Mantel-Haenszel 统 计 量 的 计 2 算 结 果, 相 当 于 χ MH ( 即 Mantel-Haenszel χ 2 ), χ 2 MH = 9.7038,P = 0.0018, 说 明 总 体 RR 与 1 之 间 的 差 别 有 统 计 学 意 义 第 二 部 分 是 普 通 相 对 风 险 度 的 估 计 值, 本 资 料 属 于 队 列 研 究 资 料, 因 此 需 要 看 Cohort 对 应 的 RR 值, 第 1 列 ( 复 发 ) 的 相 对 危 险 度 RR = 0.2750, 其 95% 置 信 区 间 为 (0.1014,0.7458) 统 计 和 专 业 结 论 : 因 本 例 资 料 不 满 足 χ 2 检 验 的 前 提 条 件, 故 根 据 一 般 结 果 可 知,P = 0.0028 < 0.01, 按 照 0.05 的 检 验 水 准, 拒 绝 原 假 设 ( 两 药 物 治 疗 对 应 的 胃 炎 复 发 率 相 等 ), 接 受 备 择 假 设 ( 两 药 物 治 疗 对 应 的 胃 炎 复 发 率 不 相 等 ), 可 以 认 为 不 同 药 物 治 疗 胃 炎 的 复 发 率 不 同 根 据 RR 值 及 其 置 信 区 间, 可 以 认 为, 一 段 时 间 后, 使 用 A 药 治 疗 胃 炎 复 发 的 风 险 是 B 药 的 0.2750 倍, 换 言 之, 相 对 于 A 药 而 言,B 药 治 疗 胃 炎 更 易 于 复 发 值 得 注 意 的 是 : 输 出 结 果 中 的 ( 优 比 ) 部 分 的 结 果 仅 当 四 格 表 资 料 来 自 病 例 - 对 照 研 究 设 计 时 才 用 另 外, 与 队 列 研 究 设 计 四 格 表 资 料 分 析 结 果 对 应 的 有 两 部 分, 第 一 部 分 标 有 ( 第 一 列 风 险 ), 意 思 是 与 第 1 列 对 应 的 相 对 危 险 度 计 算 结 果 第 二 部 分 标 有 ( 第 二 列 风 险 ), 意 思 是 与 第 2 列 对 应 的 相 对 危 险 度 计 算 结 果 若 将 资 料 输 入 程 序 时, 第 1 列 上 代 表 研 究 者 关 心 的 结 果, 其 RR 的 计 算 结 果 就 应 采 用 第 一 部 分 的 内 容, 反 之, 就 应 采 用 第 二 部 分 的 内 容 分 析 与 SAS : 上 表 资 料 为 来 自 病 例 - 对 照 研 究 设 计 的 2 2 列 联 表 资 料, 因 总 频 数 n = 1418 > 40, 且 无 小 于 5 的 理 论 频 数, 可 先 按 横 断 面 研 究 设 计 2 2 列 联 表 资 料 的 检 验 方 法 考 察 两 组 频 数 分 布 是 否 相 同, 如 果 不 相 同, 则 进 一 步 对 总 体 OR( 优 势 比 ) 是 否 等 于 1 进 行 检 验 具 体 的 SAS 程 序 如 下 : data tjfx3 do a=1 to 2 do b=1 to 2 input f @@ output endend tables a*b/chisq cmh ods html close 688 650 21 59 程 序 说 明 : 建 立 名 为 tjfx3 的 数 据 集, 建 立 数 值 变 量 a b f, 分 别 读 入 行 号 列 号 每 格 实 际 频 数 过 程 步, 调 用 FREQ 过 程, 指 定 频 数 变 量 为 f, 用 tables a*b 语 句 表 示 生 成 二 维 列 联 表 资 料, 斜 杠 / 后 的 chisq 和 cmh 分 别 要 求 进 行 一 般 卡 方 检 验 和 cmh χ 2 检 验 SAS 输 出 结 果 与 结 果 解 释 : 卡 方 1 19.1292 <.0001 似 然 比 卡 方 1 19.8780 <.0001 连 续 校 正 卡 方 1 18.1357 <.0001 Mantel-Haenszel 卡 方 1 19.1157 <.0001 Phi 系 数 0.1161 列 联 系 数 0.1154 Cramer 的 V 0.1161 3 病 例 - 对 照 研 究 设 计 2 2 表 资 料 统 计 分 析 的 SAS 为 什 么 有 些 人 会 患 病, 而 另 一 些 人 却 没 有 患 病? 研 究 者 依 据 专 业 知 识, 提 出 产 生 该 病 的 可 疑 危 险 因 素, 寻 找 与 患 者 条 件 相 近 的 另 一 组 人 作 为 正 常 对 照 者, 对 他 们 进 行 回 顾 性 调 查, 并 按 此 危 险 因 素 的 暴 露 与 否 将 调 查 者 细 分 为 两 类, 记 录 两 组 中 暴 露 与 非 暴 露 的 人 数, 这 样 收 集 的 资 料 称 为 来 自 病 例 - 对 照 研 究 设 计 的 2 2 列 联 表 资 料 例 3 一 项 关 于 吸 烟 与 肺 癌 关 系 的 病 例 对 照 研 究, 资 料 见 表 3, 试 探 讨 吸 烟 史 对 肺 癌 发 生 的 影 响 [3] 吸 烟 史 表 3 吸 烟 史 在 肺 癌 及 非 肺 癌 组 中 的 情 况 病 例 对 照 合 计 有 688 650 1338 无 21 59 80 合 计 709 709 1418 单 元 格 (1,1) 频 数 (F) 688 左 侧 Pr <= F 1.0000 右 侧 Pr >= F 7.382E-06 表 概 率 (P) 5.016E-06 双 侧 Pr <= P 1.476E-05 以 上 是 卡 方 检 验 及 法 对 两 组 频 数 分 布 是 否 相 同 的 检 验 结 果,χ 2 = 19.1292,P < 0.05,Fisher 精 确 检 验 法 得 P < 0.05( 双 侧 ), 说 明 两 组 频 数 分 布 不 同 Cochran-Mantel-Haenszel 统 计 量 ( 基 于 表 得 分 ) 统 计 量 对 立 假 设 自 由 度 值 概 率 1 非 零 相 关 1 19.1157 <.0001 2 行 均 值 得 分 差 值 1 19.1157 <.0001 3 一 般 关 联 1 19.1157 <.0001
中 国 医 药 生 物 技 术 2012 年 4 月 第 7 卷 第 2 期 Chin Med Biotechnol, April 2012, Vol. 7, No. 2 157 以 上 三 种 检 验 方 法 是 对 总 体 OR 是 否 等 于 1 的 假 设 检 验,P < 0.05, 说 明 总 体 OR 1 普 通 相 对 风 险 的 估 计 值 ( 行 1/ 行 2) 研 究 类 型 方 法 值 95% 置 信 限 案 例 对 照 Mantel-Haenszel 2.9738 1.7867 4.9494 ( 优 比 ) Logit 2.9738 1.7867 4.9494 Cohort Mantel-Haenszel 1.9589 1.3517 2.8387 ( 第 1 列 风 险 ) Logit 1.9589 1.3517 2.8387 Cohort Mantel-Haenszel 0.6587 0.5716 0.7591 ( 第 2 列 风 险 ) Logit 0.6587 0.5716 0.7591 本 例 属 于 病 例 - 对 照 研 究 设 计, 故 只 需 看 前 两 行 的 结 果 OR = 2.9738, 其 95% 置 信 区 间 为 (1.7867,4.9494), 不 包 含 1, 说 明 总 体 OR 与 1 之 间 的 差 别 有 统 计 学 意 义 统 计 和 专 业 结 论 : 这 是 病 例 - 对 照 研 究 设 计 的 2 2 表 资 料, 总 体 OR 与 1 之 间 的 差 别 有 统 计 学 意 义, 计 算 得 OR = 2.9738, 其 95% 置 信 区 间 为 (1.7867,4.9494), 可 认 为 吸 烟 史 是 导 致 肺 癌 的 危 险 因 素, 并 且 有 吸 烟 史 的 人 患 肺 癌 的 概 率 对 于 无 吸 烟 史 的 人 患 肺 癌 的 概 率 优 势 比 OR 值 为 2.9738, 也 就 是 说 有 吸 烟 史 的 人 患 肺 癌 的 危 险 度 是 无 吸 烟 史 人 的 2.9738 倍 4 配 对 研 究 设 计 2 2 表 资 料 统 计 分 析 的 SAS 所 谓 配 对 研 究 设 计 2 2 表 资 料 是 指 按 照 配 对 原 则, 受 试 者 或 样 品 分 别 接 受 两 种 不 同 的 处 理 方 法, 每 种 处 理 方 法 的 结 果 都 可 分 为 阳 性 和 阴 性 两 种, 数 出 两 种 处 理 方 法 同 时 判 定 为 阳 性 阴 性 的 频 数 以 及 它 们 结 果 不 一 致 的 频 数, 将 结 果 列 成 配 对 设 计 的 2 2 表 形 式 其 中 一 种 处 理 方 法 区 分 为 阳 性 阴 性 的 标 准 即 为 金 标 准, 医 学 中 常 见 的 金 标 准 有 活 组 织 检 查 尸 体 解 剖 手 术 探 查 和 跟 踪 随 访 结 果 等, 另 一 种 处 理 方 法 区 分 阳 性 阴 性 的 结 果 与 金 标 准 区 分 结 果 进 行 比 较, 检 验 其 区 分 能 力 例 4 为 了 探 讨 噬 菌 体 生 物 扩 增 法 (PhaB 法 ) 直 接 检 测 痰 结 核 分 支 杆 菌 (MTB) 对 异 烟 肼 的 耐 药 性 的 效 果, 以 建 立 一 种 特 异 性 强, 敏 感 性 高 的 MTB 药 敏 试 验 方 法 现 有 临 床 上 初 步 诊 断 为 耐 药 肺 结 核 的 40 例 患 者 痰 标 本, 用 PhaB 法 和 绝 对 浓 度 法 分 别 检 测 痰 MTB 对 异 烟 肼 的 耐 药 性, 结 果 见 表 4 以 绝 对 浓 度 法 检 测 结 果 为 判 断 标 准 试 比 较 两 种 方 法 检 测 结 果 之 间 不 一 致 部 分 的 差 别 是 否 具 有 统 计 学 意 义, 并 检 测 其 一 致 性 是 否 具 有 统 计 学 意 义 [4] 表 4 PhaB 法 与 绝 对 浓 度 法 检 测 结 果 比 较 PhaB 法 绝 对 浓 度 法 耐 药 敏 感 合 计 耐 药 22 1 23 敏 感 6 11 17 合 计 28 12 40 分 析 与 SAS : 绝 对 浓 度 法 是 检 测 痰 结 核 分 支 杆 菌 对 异 烟 肼 耐 药 性 的 金 标 准, 故 该 资 料 为 特 设 金 标 准 的 配 对 设 计 四 格 表 资 料 值 得 注 意 的 是 : 若 绝 对 浓 度 法 不 是 金 标 准, 则 仅 当 检 测 结 果 中 隐 含 金 标 准 ( 即 两 种 检 测 方 法 得 到 的 阳 性 结 果 都 是 真 阳 性, 不 存 在 假 阳 性 结 果 ) 时, 可 以 按 如 下 的 方 法 进 行 统 计 分 析 若 既 没 有 金 标 准 又 不 属 于 隐 含 金 标 准 的 试 验 研 究 问 题, 则 没 有 分 析 的 价 值 根 据 不 同 的 分 析 目 的, 有 两 种 统 计 分 析 方 法 ⑴ 若 原 文 的 目 的 是 为 了 比 较 两 种 检 测 方 法 不 一 致 部 分 的 频 数 之 间 的 差 别 是 否 具 有 统 计 学 意 义, 应 选 用 配 对 设 计 定 性 资 料 的 χ 2 检 验, 即 McNemar χ 2 检 验 McNemar χ 2 检 验 分 析 的 SAS 程 序 如 下 : data tjfx4 do a=1 to 2 do b=1 to 2 input f @@ tables a*b/agree output endend ods html close 22 1 6 11 程 序 说 明 : 数 据 步, 建 立 名 为 tjfx4 的 数 据 集, 建 立 数 值 变 量 a b f, 分 别 读 入 行 号 列 号 每 格 实 际 频 数 过 程 步, 调 用 FREQ 过 程, 指 定 频 数 变 量 为 f, 用 tables a*b 语 句 表 示 二 维 列 联 表 资 料, 斜 杠 / 后 的 选 项 agree 表 示 要 求 进 行 McNemar χ 2 检 验, 同 时 计 算 简 单 的 Kappa 系 数 及 其 95% 置 信 区 间 SAS 输 出 结 果 与 解 释 : McNemar 检 验 统 计 量 (S) 3.5714 自 由 度 1 Pr > S 0.0588 以 上 是 McNemar χ 2 检 验 的 结 果, 检 验 统 计 量 为 3.5714,P = 0.0588, 可 认 为 两 种 检 测 方 法 不 一 致 部 分 的 差 别 无 统 计 学 意 义 简 单 Kappa 系 数 Kappa 0.6277 渐 近 标 准 误 差 0.1229 95% 置 信 下 限 0.3868 95% 置 信 上 限 0.8685 以 上 是 Kappa 检 验 ( 即 一 致 性 检 验 ) 结 果,Kappa 值 为 0.6277, 虽 然 未 对 其 进 行 假 设 检 验, 但 这 里 给 出 了 Kappa
158 中 国 医 药 生 物 技 术 2012 年 4 月 第 7 卷 第 2 期 Chin Med Biotechnol, April 2012, Vol. 7, No. 2 值 的 95% 置 信 区 间 (0.3868,0.8685), 不 包 含 0, 说 明 Kappa 值 与 0 之 间 的 差 异 有 统 计 学 意 义 ⑵ 若 原 文 的 目 的 是 为 了 检 验 PhaB 法 绝 对 浓 度 法 两 种 方 法 的 检 测 结 果 是 否 一 致, 就 应 采 用 一 致 性 检 验, 即 Kappa 检 验 Kappa 检 验 分 析 的 SAS 程 序 如 下 : data tjfx5 do a=1 to 2 do b=1 to 2 input f @@ tables a*b output test kappa endend ods html close 22 1 6 11 程 序 说 明 : 数 据 步 与 McNemar χ 2 检 验 部 分 是 一 致 的 过 程 步 的 不 同 之 处 在 于, 去 掉 了 选 项 agree, 增 加 test Kappa 语 句 要 求 进 行 Kappa 检 验, 在 进 行 Kappa 检 验 的 同 时, 也 进 行 了 McNemar χ 2 检 验 SAS 输 出 结 果 与 结 果 解 释 : McNemar 检 验 统 计 量 (S) 3.5714 自 由 度 1 Pr > S 0.0588 这 是 用 McNemar χ 2 检 验 的 结 果, 结 论 同 前, 从 略 简 单 Kappa 系 数 Kappa 0.6277 渐 近 标 准 误 差 0.1229 95% 置 信 下 限 0.3868 95% 置 信 上 限 0.8685 以 上 是 Kappa 检 验 结 果, 结 论 同 前, 从 略 这 是 对 H 0 : 总 体 Kappa = 0 的 假 设 检 验 结 果,Z = 4.1180, 双 侧 P < 0.0001, 说 明 两 种 方 法 检 测 结 果 之 间 的 一 致 率 P 0 = (23+11)/40 = 85.0% 与 期 望 的 一 致 率 ( 相 当 于 随 机 产 生 的 )P e = (23 28/40 + 17 12/40)/40 = 53.0% 之 间 的 差 别 有 统 计 学 意 义, 即 实 际 一 致 率 比 期 望 的 要 高, 两 种 方 法 的 检 测 结 果 具 有 较 高 的 一 致 性 专 业 结 论 : 经 McNemar χ 2 检 验 可 知, 两 种 检 测 方 法 之 间 不 一 致 部 分 的 差 别 无 统 计 学 意 义 经 一 致 性 检 验 可 知, PhaB 法 与 绝 对 浓 度 法 在 检 测 痰 MTB 对 异 烟 肼 的 耐 药 性 方 面 具 有 较 高 的 一 致 性 通 常 最 好 的 结 果 是,McNemar χ 2 检 验 结 果 P > 0.05 且 两 个 不 一 致 的 频 数 都 接 近 于 0, 而 一 致 性 检 验 结 果 是 P < 0.05 且 观 察 的 一 致 率 符 合 专 业 要 求 ( 例 如, 专 业 上 要 求 观 察 的 一 致 率 大 于 85%) 此 时, 若 在 专 业 上, 新 的 方 法 较 金 标 准 经 济 便 捷, 方 可 认 为 : 可 用 新 方 法 取 代 金 标 准 参 考 文 献 [1] Hu LP, Wang Q. Scientific research design and statistical analysis of Chinese traditional medicine. Beijing: People s Medical Publishing House, 2011:199-221. (in Chinese) 胡 良 平, 王 琪. 中 医 药 科 研 设 计 与 统 计 分 析. 北 京 : 人 民 卫 生 出 版 社, 2011:199-221. [2] Hu LP, Wang Q. Statistical analysis using SAS software package for data of 2 2 contingency table. J Chin Integr Med, 2009, 7(7):678-682. (in Chinese) 胡 良 平, 王 琪. 2 2 列 联 表 资 料 的 统 计 分 析 与 SAS 软 件. 中 西 医 结 合 学 报, 2009, 7(7):678-682. [3] Li YH. Achievement of statistical analysis of case-control study contingency table data using SPSS. J Pub Health Prev Med, 2007, 18(6):114-115. (in Chinese) 李 永 红. 病 例 对 照 研 究 列 联 表 资 料 的 SPSS 分 析. 公 共 卫 生 与 预 防 医 学, 2007, 18(6):114-115. [4] Bian ZY, Yang YH, He C, et al. Study on the isoniazid susceptibility of mycobacterium tuberculosis in sputum detected by phage amplified biologically assay. Exp Lib Med, 2009, 29(3):219-220, 334. (in Chinese) 边 泽 源, 杨 芸 红, 何 川, 等. 噬 菌 体 生 物 扩 增 法 直 接 检 测 痰 标 本 中 结 核 分 枝 杆 菌 异 烟 肼 耐 药 性 的 研 究. 实 验 与 检 验 医 学, 2009, 29(3): 219-220, 334. H 0 检 验 :Kappa = 0 H 0 下 的 渐 近 标 准 误 差 0.1524 Z 4.1180 单 侧 Pr > Z <.0001 双 侧 Pr > Z <.0001