第 七 章 相 关 分 析
作者：赵耐青 授权刊登：医学统计之星
本 节STATA 命 令 摘 要
|
correlate 变量名1 变量名2 … 变量名m spearman 变量1 变量2 |
· 线 性 关 系 的 相 关 分 析
若 计 量 资 料 变 量1，…， 变 量m 服 从 正 态 分 布， 对 于 它 们 之 间 是 否 存 在 线 性 相 关 关 系 可 以 通 过 相 关 分 析 方 法， 相 应 的 STATA 命 令 为：
pwcorr 变量名1 变量名2 … 变量名m, sig
例： 上 海 医 科 大 学 儿 科 医 院 研 究 某 种 代 乳 粉 的 营 养 价 值 时， 用 大 白 鼠 作 试 验， 得 大 白 鼠 进 食 量( 克) 和 增 加 体 重( 克) 之 间 的 关 系 的 原 始 数 据 如 下 表， 用 x 表 示 大 白 鼠 进 食 量 和 用 y 表 示 大 白 鼠 增 加 体 重，试 作 相 关 分 析( 资 料 摘 自 医 学 统 计 分 析， 金 丕 焕 主 编，p101)。
|
进食量(克) |
820 |
780 |
720 |
867 |
690 |
787 |
934 |
679 |
639 |
820 |
|
增加体重(克) |
165 |
158 |
130 |
180 |
134 |
167 |
186 |
145 |
120 |
158 |
作 关 于 x 和 y 相 关 分 析
pwcorr x y,sig
|
| x y ---------------+----------------------------- x | 1.0000 | | y | ① 0.9395 1.0000 | ② 0.0001 | |
① 为 x 和 y 的 相 关 系 数 r。 ② 为 相 关 系 数 显 著 性 检 验(Ho:r=0) 所 对 应 的 p 值。 因 此 根 据 本 例 资 料 和 上 述 结 果 可 以 认 为 大 白 鼠 进 食 量 与 大 白 鼠 体 重 增 加 呈 线 性 正 相 关( 若 数 据 x 和 y 服 从 正 态 分 布)。
· 秩 和 相 关 分 析
由 于 以 上 的 相 关 分 析 中， 要 求 数 据 x 和 y 均 服 从 正 态 分 布，但 是 对 于 那 些 不 能 满 足 正 态 分 布 的 数 据 之 间 的 相 关 分 析 应 采 用 非 参 数 相 关 分 析， 通 常 采 用 秩 和 相 关 分 析， 即： Spearman 相 关 分 析 方 法。STATA 命 令 为：
spearman 变 量1 变 量2
以 上 例 资 料 为 例， 作 秩 和 相 关 分 析：
spearman x y
|
Number of obs = 10 Spearman's rho = 0.8994 ① Test of Ho: x and y independent ② ③ Pr > |t| = 0.0004 |
① 为Spearman 相 关 系 数；② 为 无 效 假 设Ho:x 与 y 独 立；③ 无 效 假 设 的Spearman 相 关 显 著 性 检 验 所 对 应 的 p 值。
由 上 述 结 果 表 明：Spearman 相 关 系 数 为 0.8994， 相 应 的 p 值 为 0.0004， 因 此 由 本 例 资 料 和 检 验 结 果 可 以 认 为 大 白 鼠 进 食 量 与 大 白 鼠 体 重 增 加 呈 正 相 关。