如何调整t-test变量顺序，计算6缸与4缸Min.Price均值差的90%置信区间？

问题背景 #

我正在处理一项任务，已从MASS包的Cars93数据框中筛选得到仅包含4缸和6缸车型的数据框df.3。我尝试使用t-test估算这两类车型Min.Price列均值差的90%置信区间，但默认t-test计算的是4缸减6缸的差值，请问如何改为计算6缸减4缸的差值？已执行的代码如下：

var.test(Min.Price ~ Cylinders, df.3)

t.test(Min.Price ~ Cylinders, df.3, conf=0.9)
t.test(df.3$Min.Price ~ df.3$Cylinders,conf=0.9)$conf.int

三种解决方法 #

• 方法一：修改因子水平顺序
  t-test默认按因子水平的先后顺序计算差值（第一组减第二组），把Cylinders的水平调整为"6"在前、"4"在后即可：

  # 重新设定因子水平顺序
  df.3$Cylinders <- factor(df.3$Cylinders, levels = c("6", "4"))
  # 执行t-test获取6缸减4缸的置信区间
  t.test(Min.Price ~ Cylinders, df.3, conf = 0.9)$conf.int
  

  这个方法会修改原数据的因子属性，后续所有基于该因子的分析都会以6缸为基准组，适合需要重复分析的场景。

• 方法二：直接指定两组数据顺序
  不修改原数据，手动把6缸和4缸的Min.Price分别作为t-test的第一、第二个参数：

  t.test(
    x = df.3$Min.Price[df.3$Cylinders == "6"],
    y = df.3$Min.Price[df.3$Cylinders == "4"],
    conf = 0.9
  )$conf.int
  

  这种方式更灵活，仅在当前t-test中调整顺序，不影响原数据。

• 方法三：对默认结果取反
  因为6缸减4缸的差值，正好是默认结果（4缸减6缸）的相反数，直接对已得到的置信区间取反即可：

  # 获取默认的4缸减6缸置信区间
  default_ci <- t.test(Min.Price ~ Cylinders, df.3, conf = 0.9)$conf.int
  # 转换为6缸减4缸的置信区间
  rev_ci <- default_ci * -1
  rev_ci
  

  如果已经算出了默认结果，这个方法最快。

内容的提问来源于stack exchange，提问作者your null hypothesis