R Markdown

set.seed(123)
brix <- sort.int(rnorm(54, 21, 2), 30)
mo <- sort.int(rnorm(54,2,0.25), 30)
mo_cut <- cut(mo, quantile(mo, c(0, 1/3, 2/3, 1)),include.lowest = T, labels = c('bloq1','bloq2','bloq3')) 
mo_cut <- sort(mo_cut)
variedad <- gl(3, 6, 54, c('v1','v2','v3'))
table(mo_cut, variedad)
##        variedad
## mo_cut  v1 v2 v3
##   bloq1  6  6  6
##   bloq2  6  6  6
##   bloq3  6  6  6
mod <- aov(brix ~ mo_cut * variedad) 
summary(mod)
##                 Df Sum Sq Mean Sq F value   Pr(>F)    
## mo_cut           2 101.19   50.59  71.136 1.17e-14 ***
## variedad         2  25.04   12.52  17.601 2.25e-06 ***
## mo_cut:variedad  4  16.99    4.25   5.974 0.000606 ***
## Residuals       45  32.00    0.71                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
df <- data.frame(variedad, mo_cut, brix)

df %>%
  ggpubr::ggline(
    x            = "mo_cut", 
    y            = "brix", 
    color        = "variedad",
    add          = c("mean_se"),
    palette      = c("#00AFBB", "#E7B800", "red"),
    xlab         = "Grupos",
    ylab         = "brix",
    legend.title = "graf")
## Warning: Computation failed in `stat_summary()`:
## el objeto 'mean_se_' de modo 'function' no fue encontrado

interaction.plot(mo_cut,variedad,brix)

En el bloque 1 (niveles bajos de MO) todas las variedades se comportan iguales. En niveles intermedios la variedad 1 se ve afectada y la que mejor se comporta es la v3. En niveles altos la variedad 3 es la mejor, pero no parece haber diferencias entre la v1 y v2.

ggplot(df)+
  aes(mo_cut, brix, color = variedad)+
  geom_jitter()