计算单样本 Wilcoxon 检验的中位数估计值和置信区间

中位数估计值

假设下面的数据位于 C1 中,且 C2 和 C3 为空:

24 12 24 19 12 21 23 11 17 19 23 14 23 15 14 6

  1. 计算 W(i) = (Xi + Xj)/2 = i 的所有配对平均数 j.
    1. 选择统计 > 非参数 > 配对平均数
    2. 变量中,输入 C1
    3. 将平均数存储在中,键入 C2。单击确定
  2. 计算中位数估计值
    1. 选择统计 > 基本统计 > 存储描述性统计量
    2. 变量中,输入 C2
    3. 单击统计量,并且仅选中中间距离
    4. 在每个对话框中单击确定
    所得到的存储值 (17.5) 是单样本 Wilcoxon 检验的中位数估计值。

置信区间

假设下面的数据位于 C1 中,且 C2 和 C3 为空:

24 12 24 19 12 21 23 11 17 19 23 14 23 15 14 6

  1. 计算 W(i) = (Xi + Xj)/2 = i 的所有配对平均数 j.
    1. 选择统计 > 非参数 > 配对平均数
    2. 变量中,输入 C1
    3. 将平均数存储在中,键入 C2。单击确定
  2. 按从小到大的数字顺序排列配对平均数。
    1. 选择数据 > 排序
    2. 要作为排序条件的列中的下面,输入 C2
    3. 要排序的列中,选择指定的列
    4. 中,输入 C2
    5. 已排序列的存储位置下面,选择在原始列中。单击确定
    排序的配对平均数位于 C2 中。
  3. 要获得 (1-α)*100% 置信区间的端点,请首先为 Z(1-α/2) 求解。对于 95% 的置信区间:
    1. 选择计算 > 概率分布 > 正态
    2. 选择逆累积概率
    3. 选择输入常量并输入 0.975。单击确定
  4. 接下来,计算 d,它约为
    1. 选择计算 > 计算器
    2. 将结果存储在变量中中,输入 C3
    3. 表达式中,输入 16*17/4-.05-1.96*sqrt(16*17*33/24)。单击确定
  5. 使用上面的表示法,置信区间的下端点为 W(d+1),上端点为 W(nw-d),其中 nw 为配对平均值的数量。

    W(d+1) = W(31)。C2 中的第 31 个配对平均值为 14.5。

    W(nw-d) = W(136-30) = W(106)。C2 中的第 106 个配对平均值为 21。

    Wilcoxon 符号秩置信区间: C1

    方法 η: C1 的中位数
    描述性统计量 η 的置 取得的 样本 N 中位数 信区间 置信度 C1 16 17.5 (14.5, 21) 94.75%
    Note

    由于 d(与 Wilcoxon 检验统计量接近)为正值,因此很难实现指定的置信。该过程打印最接近的值,该值是结合使用正态近似和连续校正计算得出的。