Função que faz gráfico a partir de uma fórmula para análise ANOVA

gistfile1.txt

grafico_serie_anova<- function(formula,  dados_origem){ 
    
    library(agricolae)
    library(ggplot2)
    
    res.aov <- aov( formula, data =dados_origem )
    
    out<-HSD.test(res.aov,res.aov[["terms"]][[3]])
    
    grupos <- out$groups[2]
    grupos$ano <- rownames(out$groups)
    
    medias<- dados_origem%>%
    inner_join(grupos) 
    
    names(medias)[1] <- "data_hora"
    names(medias) [2] <- "serie"
    
    medias %>% ggplot()+
               geom_boxplot(aes(x=data_hora, y=serie, color= groups))+
               stat_summary(aes(x=data_hora, y=serie, color= groups),fun.y=mean,        geom="point", shape=45, size=8) +
    theme(legend.title = element_text(size = 10),
          legend.text = element_text(size = 12),
          plot.title = element_text(size=16),
          axis.title=element_text(size=14,face="bold"),
          axis.text.x = element_text(size=14,face="bold"),
          axis.title.x = element_blank(),
          axis.ticks.x = element_blank(),
          axis.title.y = element_blank(),
          panel.grid = element_blank(),
          panel.background = element_blank(),
          panel.border = element_blank())
  
  
    
  }
  
##Exemplo de como chamar a função
##Não executar
formula <- componente_beneficios ~ ano

dados_origem <- df_trabalho_anova_cluster_1%>%
                select(4,15) #o df de dados de origem deve ter duas colunas, a primeira é datetime e a segunda contém os valores que se quer analisar através da ANOVA

grafico_serie_anova(formula,dados_origem)