R/plotExperiments.R

In dosorio/ccCINTROP: Calculate and plot the CC50

#' @export plotExperiments
plotExperiments <- function(file, xlab=NULL, ylab=NULL){
  data <- read.table(file,header = TRUE,row.names = 1)
  data_id <- rownames(data)
  data_mean <- round(rowMeans(data),3)
  data_sd <- round(apply(data,1,sd),3)
  data_se <- round(data_sd/sqrt(ncol(data)),3)
  summaryData <- as.data.frame(cbind(data_mean,data_sd,data_se))
  anova <- analyzeExperiments(file)
  ggplot(summaryData, aes(x=data_id, y=data_mean, fill=data_id)) +
    scale_fill_grey(start = .9, end = 0) + 
    geom_bar(position=position_dodge() , stat="identity", color="gray50") +
    geom_errorbar(aes(ymin=data_mean-data_sd, ymax=data_mean+data_sd), width=.2, position=position_dodge(.9)) +
    xlab(xlab) +
    ylab(ylab) +
    guides(fill=FALSE) + 
    theme_classic()
}