手动计算标准差 (sigmas) 的 Bonferroni 置信区间

完成以下步骤来手动计算因子水平的标准差(西格玛)的 Bonferroni 置信区间,而不是使用统计 > 方差分析 > 等方差检验 当检验基于正态分布时。

假设您使用 Minitab 样本数据集。医院比较.MWX 并分析响应 评级 和因子 医院,其具有 2 个水平。要使用整体置信水平 95% (0.95) 计算医院 A 的置信区间,请执行以下操作:

  1. 计算 K 并将值存储在一个名为 K1 的常量中。
    1. 打开 Minitab 样本数据集 医院比较.MWX.
    2. 选择计算 > 计算器
    3. 将结果存储在变量中中,输入 K1
    4. 表达式中,输入 0.05 / (2 * 2)
    5. 单击确定
  2. 计算方法的每个水平的方差和样本大小,并将结果存储在工作表中。
    1. 选择统计 > 基本统计 > 显示描述性统计量
    2. 变量中,输入评级
    3. 分组变量(可选的)中,输入医院
    4. 单击统计量。选择单方差N 非缺失
    5. 单击每个对话框中的确定
    对于 医院 A,方差为 66.9579,样本大小为 20。
  3. 计算 U,并将结果存储在名为 U1 的工作表常量中。
    1. 选择计算 > 概率分布 > 卡方 カイ二乗分布の詳細については、 カイ二乗分布を参照してください。
    2. 选择逆累积概率
    3. 自由度中,键入 19
    4. 选择输入常量,并输入 K1
    5. 可选存储中,输入 U1
    6. 单击确定
  4. 计算置信区间的上限。
    1. 选择计算 > 计算器
    2. 将结果存储在变量中中,输入 UpperL1
    3. 表达式中,输入 ((19 * 66.9579) / U1)**0.5
    4. 单击确定
    医院 A 的 95% Bonferroni 置信区间的上限是 12.6776。下限的计算方法相同,使用 L 代替 U。
  5. 计算 L 并将结果存储在一个名为 L1 的工作表常量中。
    1. 选择计算 > 计算器
    2. 将结果存储在变量中中,输入 K3
    3. 表达式中,输入 1 - K1
    4. 单击确定
    5. 选择计算 > 概率分布 > 卡方
    6. 选择逆累积概率
    7. 自由度中,键入 19
    8. 选择输入列并输入 K3
    9. 可选存储中,输入 L1
    10. 单击确定
  6. 计算置信区间的下限。
    1. 选择计算 > 计算器
    2. 将结果存储在变量中中,输入 LowerL1
    3. 表达式中,输入 ((19 * 66.9579) / L1)**0.5
    4. 单击确定
    医院 A 的 95% Bonferroni 置信区间的下限是 5.99494。
    注意

    当具有多个因子时,需将每个不同的因子水平组合视为一个单独的因子水平。