Знаете ли вы, почему размеры [продукт 0] не соответствуют длине объекта [1] при попытке парыw.anova

У меня есть эти данные:

a <- c(12, 14, 16, 12, 14, 16, 12, 14, 16)
b <- c(107.13, 83.01,  164.01, 113.01, 85.48,  86.08,  112.96, 82.18,  51.92)

a — фактор с тремя уровнями 12, 14 и 16, а b — моя переменная ответа

Когда я пытаюсь бежать и строить планы pairs.anova

pairw.anova(y=a, x=b,  conf.level = 0.95, method = "tukey", MSE = NULL, df.err = NULL)

Я получаю эту ошибку:

размеры [продукт 0] не соответствуют длине объекта [1]

Знаете почему?

Спасибо за вашу помощь.


r statistics anova tukey data-comparison
person javier J    schedule 12.03.2020    source источник
comment
Вы должны добавить, из какого пакета находится функция pairw.anova. Поскольку вопрос стоит сейчас, вам трудно помочь, так как вам трудно воспроизвести вашу ошибку.   -  person henrik_ibsen    schedule 12.03.2020
comment
Привет, я использую библиотеку (asbio)   -  person javier J    schedule 27.03.2020

Ответы (1)


arrow_upward
0
arrow_downward

Возможно, вы могли бы просто поменять местами значения ваших аргументов, поскольку переменная ответа (количественная), вероятно, равна b, а ваш категориальный вектор, содержащий группы, вероятно, равен a.

И вы должны убедиться, что a является фактором:

library(asbio)

pairw.anova(y=b, x=factor(a),  conf.level = 0.95, 
            method = "tukey", MSE = NULL, df.err = NULL)

95% Tukey-Kramer confidence intervals 

               Diff      Lower    Upper Decision Adj. p-value
mu12-mu14  27.47667  -55.80207 110.7554   FTR H0     0.596586
mu12-mu16  10.36333   -72.9154 93.64207   FTR H0     0.923774
mu14-mu16 -17.11333 -100.39207  66.1654   FTR H0     0.809506
person Edward    schedule 12.03.2020
comment
Привет, Эдвард, да, ты был прав, моя ошибка, когда я печатал, я сделал то, что ты показываешь, но у меня все та же проблема. - person javier J; 13.03.2020
comment
Да, я забыл сказать вам изменить a на фактор. Извини за это! - person Edward; 13.03.2020
comment
Привет, надеюсь, что все идет хорошо. Я сделал это, но у меня та же проблема, какие библиотеки вы используете? Спасибо - person javier J; 27.03.2020
comment
У меня все хорошо спасибо! Эмм, какие библиотеки? Гоша, как давно это было. Мне нужно прочитать это снова... Кажется, я использовал пакет asbio. - person Edward; 27.03.2020
comment
Я использую библиотеку (tcltk) библиотеку (mvtnorm) библиотеку (asbio) библиотеку (multcompView) - person javier J; 27.03.2020
comment
Итак, я думаю, что использую функцию pairw.anova из пакета asbio, потому что вы не упомянули, какой пакет вы использовали в вопросе (так что я догадался). Попробуйте использовать этот пакет. - person Edward; 27.03.2020