stat_compare_means anova в выбранных группах

Я хотел бы рассчитать только anova, включая gr1-gr3, исключая gr4. Возможно ли это с помощью stat_compare_means()?

x   value
gr1 3.543
gr1 2.12
gr1 4.56
gr2 3.3
gr2 2.12
gr2 4.1
gr3 3.32
gr3 2.3
gr3 3.3
gr4 3.325
gr4 2.355
gr4 3.34

library(ggpubr)
ggplot(df(aes(x,value)) + geom_boxplot() + stat_compare_means()

r ggplot2 ggpubr
person user2300940    schedule 23.10.2018    source источник
comment
Вы всегда можете отфильтровать это. df %>% filter(x != "gr4").   -  person TheRimalaya    schedule 23.10.2018
comment
Я буду держать группу в сюжете   -  person user2300940    schedule 23.10.2018

Ответы (1)


arrow_upward
3
arrow_downward

Это может помочь. Для этого можно использовать аргумент comparisons.

ggplot(dta, aes(x, value)) + 
    geom_boxplot() + 
    stat_compare_means(comparisons = list(1:2, 2:3, c(1, 3), c(1, 2, 3)))

введите описание изображения здесь

РЕДАКТИРОВАТЬ: для дисперсионного анализа

ggplot(dta, aes(x, value)) + 
    geom_boxplot() + 
    stat_compare_means(method = "anova", label.y = 4.8, 
                       data = dta %>% filter(x != "gr4"))

введите описание изображения здесь

person TheRimalaya    schedule 23.10.2018
comment
Это не дает ANOVA - person user2300940; 23.10.2018
comment
Попробуйте использовать method="anova" и дать отфильтрованный набор данных. - person TheRimalaya; 23.10.2018