У меня есть датафрейм df с 2 столбцами:
df <- data.frame(runif(80,0.5,1), as.factor(rep(1:8, 10)))
names(df) = c('A','B')
Я хотел бы создать коробочные диаграммы, которые показывают результаты для специального теста Тьюки на разных уровнях факторов с пакетом multcompview (можно добавить метки на график, указывающие различные группы из теста posthoc результаты)
> multcompBoxplot(A ~ B, data = df, horizontal = FALSE,
compFn = "TukeyHSD",
sortFn = NULL,
plotList=list(
boxplot=list(fig=c(0.1, 0.98, 0, 1), las=1,ylim=c(0.5,1),
cex.axis=0.8),
multcompLetters=list(fig=c(0.1, 0.12, 0.15, 0.9))))
Дает:
Ошибка в xy.coords (x, y, xlabel, ylabel, log):
«x» - это список, но не содержит компонентов «x» и «y»
Почему x рассматривается как список, когда он записан в виде формулы (A ~ B)? Или я что-то упускаю? Также: когда я использую данные для df $ A, содержащие нули, все работает отлично. График изготовлен кстати, но без меток. Запуск только с первыми аргументами multcompBoxplot(A ~ B, data = df)
приводит к той же ошибке.
Пожалуйста, помогите,
Спасибо!