配对 t.test 分组因子必须恰好有 2 个级别错误

Question

我为一个项目创建了一个调查。调查参与者使用李克特量表 0-10 进行回答。这是之前和之后的场景。数据被清理以删除不完整的响应。我还有 1135 个配对响应。数据显示数值。我使用 t.test(Q15_2~Q15_3,paired = TRUE) 来运行我的测试。

t.test.formula(Q15_2 ~ Q15_3,paired = T) 中的错误： 分组因子必须正好有 2 个级别

再次...数据已被审核了十几次。之前和之后的响应是相等的。这确实需要成对的 t 检验，所以我有目的地使用 ~.我可以让它在 excel 中正常工作，但我总是喜欢通过在 Excel 和 R 中运行它们来仔细检查我的答案，以确保它们是正确的。欢迎任何意见！

Answer 1

我认为你应该使用：

t.test(Q15_2, Q15_3, paired = TRUE) 

用逗号而不是波浪号。如果两个向量的列中都包含数据对，则应用此逗号语法。如果左侧包含两个集合的数据，而右侧向量是具有相应两个处理的因子，则应用带有~ 的公式语法。

例子：

x1 <- c(4.71, 6.66, 5.01, 4.6, 3.98)
x2 <- c(5.42, 7.1, 5.52, 6.05, 6.23)
t.test(x1, x2, paired=TRUE)

    Paired t-test

data:  x1 and x2
t = -3.1096, df = 4, p-value = 0.03589
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -2.029145 -0.114855
sample estimates:
mean of the differences 
                 -1.072 

在 Excel 中使用相同的数据（复制的值、从当地语言翻译成英文的标签）：

t-Statistics                       -3.109611652
P(T<=t) one sided                   0.017942691
critical t-value (one sided t-Test) 2.131846786
P(T<=t) two-sided                   0.035885381
critical t-value (two sided t-Test) 2.776445105

我们看到两侧假设的 p 值是相同的。