Lesson 7 Box Plot in R Multiple Box plot

>> YOUR LINK HERE: ___ http://youtube.com/watch?v=o9CCoCgC8BU

Box Plot in R • Multiple Box plot • setwd( F:/R ) • library(tidyr) • library(dplyr, warn.conflicts = FALSE) • library(purrr) • library(agricolae) • library(ggplot2) • library(devtools) • library(ggbiplot) • library(ggfortify) • library(corrplot) • library(RColorBrewer) • library(ggpubr) • library(ggplot2) • library(magrittr) • library(gridExtra);library(gtable);library(grid) • library(reshape) • library(dplyr); library(plyr) • library(ggpubr) • library(ggplot2) • library(gridExtra);library(gtable);library(grid) • library(reshape) • library(plyr); library(dplyr) • library(lme4) • library(lmerTest) • library(agricolae) • library(scatterplot3d) • library( ScapeMCMC ) • library( car ) • library( MASS ) • library( mlmRev ) • library( agridat ) • library( MCMCglmm ) • library(devtools) • library(ggbiplot) • library(corrplot) • library(RColorBrewer) • library( tidyverse ) • library(data.table) • library(heatmap.plus) • library(RColorBrewer) • library(gplots) • library(pheatmap) • library(devtools) • library(ggbiplot) • d1=read.csv( dataA.csv ) • p1=boxplot(RHL~Genotype, data=d1, • col=(c( #882255 , #44AA99 , #66AA11 , Orange , pink )), • main= Nitrogen , xlab= Genotypes ) • p1 • • • d2=read.csv( dataB.csv ) • A=ggplot(data = d2,aes(x=Genotypes, y=Gene, fill=Nitrogen)) • B=A+stat_boxplot(geom ='errorbar', width = 0.3,size=1.5,position=position_dodge(width=0.6),stat= identity ) • C=B+geom_boxplot(width = 0.6,outlier.size=4,outlier.shape =1,lwd=1.5) • F=C+scale_y_continuous(sec.axis = dup_axis(label = NULL, • name = NULL), • expand = expansion(mult = c(0, 0)), • breaks = pretty(c(0,50), n = 5), • limits = c(0,50)) • G=F+xlab( Genotypes ) • H=G+ylab( Relative expression ) • I=H+scale_fill_discrete(name = Nitrogen level ) • J=I+theme(panel.background = element_rect(fill = 'white', colour = 'white'),plot.margin = margin(1, -1, 1, 1, cm ),text = element_text(size=20,face= bold ),axis.line = element_line(color= black , size =1.5),axis.text.x = element_text(colour= black ,size=36,angle=0,hjust=.1,vjust=.1,face= bold ),axis.text.y = element_text(colour= black ,size=36,angle=0,hjust=.5,vjust=.5,face= bold )) • K=J+theme(axis.text=element_text(size=36),axis.title=element_text(size=36,face= bold )) • L=K+ theme(legend.position= top )+ theme( • axis.ticks.length.y = unit(.3, cm ), • axis.ticks.length.x = unit(.3, cm )) • O=L+ theme(axis.ticks = element_line(size = 1.5)) • P=O+scale_fill_brewer(palette= Set2 ) • P

#############################

New on site