手动计算标准差 (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。
    注意

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