R/bp_function.R

In AgronomiaR/IS:

#' boxplot
#' @export
#' @example
#' resp1=c(10,11,12,7,6,5,4,4,5,12,7,16,5,4,4,5)
#' resp2=c(10,11,12,7,6,5,4,5,5,12,7,6,15,4,4,5)
#' resp3=c(10,11,12,7,6,5,4,4,5,12,7,6,10,4,7,5)
#' resp4=c(6,1,2,7,6,3,4,3,5,12,7,6,5,14,4,15)
#' trat=rep(paste("Genotipo",1:16))
#' dados=data.frame(trat,resp1,resp2,resp3,resp4)
#' library(IS)
#' a=mulamba(dados)
#' b=plot_mulamba2(a,dados)
#' bp(b,a,dados)

bp=function(b,a,dados){
  sel=b[[1]]$data$sel
  var=b[[1]]$data$vari
  resp=unlist(c(dados[-1]))
  data=data.frame(var,sel,resp)
  library(ggplot2)
  graph=list()
  for(i in 1:length(unique(var))){
  gra=ggplot(data[var==unique(var)[i],],aes(y=resp,x=sel,fill=sel))+
    geom_boxplot(show.legend = F)+
    stat_boxplot(geom ='errorbar', width = 0.4) +
    ylab(unique(var)[i])+
    theme_bw()+xlab("")+
    theme(axis.text = element_text(size=12,color="black"),
          axis.title = element_text(size=13,color="black"))+
    scale_fill_manual(values=c("lightgreen","darkgreen"))
  graph[[i]]=gra
  }
  n <- length(graph)
  nCol <- floor(sqrt(n))
  do.call("grid.arrange", c(graph, ncol=nCol))
  graficos=graph}