#' @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()
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.