Парный фактор группировки 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 с теми же данными (скопированные значения, метки переведены с местного языка на английский sh):

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-значения для двусторонней гипотезы одинаковы.