library(ggplot2) library(cowplot) library(plyr) library(scales) library(gridExtra) library(grid) library(ggplotify) library(ggforce) cbPalette <- c('Sn'='#ffb366', 'Pr'='#2d8659') dat <- read.table('output/statsFiles/all.noFilt.tmerge.min2reads.vs.SIRVs.stats.tsv', header=T, as.is=T, sep='\t') dat <- subset(dat, seqTech=='ont-Cshl-CapTrap') dat <- subset(dat, capDesign=='HpreCap') dat <- subset(dat, sizeFrac=='0+') dat <- subset(dat, sampleRep=='Brain01Rep1') dat$seqTech <- gsub('-', '\n', dat$seqTech) dat$sampleRep <- gsub('HpreCap_', '', dat$sampleRep) horizCats <- length(unique(dat$capDesign)) * length(unique(dat$sampleRep)) vertCats <- length(unique(dat$seqTech)) wXyPlot = (horizCats * 0.9) +1.7 hXyPlot = (vertCats * 0.6) + 1.7 geom_textSize=1.4 themeSize = (14/5) * geom_textSize # https://stackoverflow.com/questions/25061822/ggplot-geom-text-font-size-control/25062509 lineSize=geom_textSize/8 minorLineSize=lineSize/2 plotBase <- "p <- ggplot(dat, aes(x=level, y=value)) + #geom_mark_rect(aes(filter = level == 'Transcript' & metric == 'Pr'), size=2, expand = unit(7, 'mm'), radius = unit(4, 'mm'), color='#4d4d4d', fill='#ffff00') + geom_point(aes(color=metric), shape=18, size=lineSize*15, alpha=0.7) + scale_colour_manual (values=cbPalette, name='Metric', breaks=c('Sn', 'Pr'))+ ylab('Sn | Pr (%)') + xlab('Evaluation level') + scale_x_discrete(limits=c('Base', 'Locus', 'Exon', 'Intron', 'Intronchain', 'Transcript')) + scale_y_continuous() + expand_limits(y=c(0,100))+ theme(axis.text.x = element_text(angle = 45, hjust = 1)) + theme(axis.text= element_text(size=themeSize*1.8), axis.ticks = element_line(size=lineSize), axis.line = element_line(colour = '#595959', size=lineSize), axis.title=element_text(size = themeSize*2), panel.grid.major = element_line(colour='#d9d9d9', size=lineSize),panel.grid.minor = element_line(colour='#e6e6e6', size=minorLineSize),panel.border = element_blank(),panel.background = element_blank(), strip.background = element_rect(colour='#737373',fill='white'), legend.key.size=unit(0.5,'line'), legend.title=element_text(size=themeSize*1.2), legend.text=element_text(size=themeSize), strip.text = element_text(size = themeSize)) + " plotFacetXy <- parse(text =paste(plotBase, "facet_grid( seqTech ~ capDesign + sampleRep, scales='free_y')")) plotFacetYx <- parse(text=paste(plotBase, "facet_grid( capDesign + sampleRep ~ seqTech, scales='free_y')")) plotFacetXy <- parse(text =paste(plotFacetXy, " + theme(strip.text.x = element_blank(), strip.text.y = element_blank())")) plotFacetYx <- parse(text=paste(plotFacetYx, " + theme(strip.text.x = element_blank(), strip.text.y = element_blank())")) pXy <- eval(plotFacetXy) pYx <- eval(plotFacetYx) legend <- get_legend(pXy) pXyNoLegend <- pXy + theme(legend.position='none') pYxNoLegend <- pYx + theme(legend.position='none') legendOnly <- grid.arrange(legend) pXyGrob <- as.grob(pXy) pYxGrob <- as.grob(pYx) pXyNoLegendGrob <- as.grob(pXyNoLegend) pYxNoLegendGrob <- as.grob(pYxNoLegend) hLegendOnly <- convertUnit(sum(legend$heights), 'in', valueOnly=TRUE) wLegendOnly <- convertUnit(sum(legend$widths), 'in', valueOnly=TRUE) hYxPlot <- wXyPlot wYxPlot <- hXyPlot hXyNoLegendPlot<- hXyPlot wXyNoLegendPlot<- wXyPlot - wLegendOnly hYxNoLegendPlot<- hYxPlot wYxNoLegendPlot<- wYxPlot - wLegendOnly save_plot('output/plots/tmerge.vs.SIRVs.stats/ont-Cshl-CapTrap/HpreCap/ont-Cshl-CapTrap_HpreCap_0+_Brain01Rep1.noFilt.tmerge.min2reads.vs.SIRVs.stats.legendOnly.png', legendOnly, base_width=wLegendOnly, base_height=hLegendOnly) save_plot('output/plots/tmerge.vs.SIRVs.stats/ont-Cshl-CapTrap/HpreCap/ont-Cshl-CapTrap_HpreCap_0+_Brain01Rep1.noFilt.tmerge.min2reads.vs.SIRVs.stats.xy.wLegend.png', pXy, base_width=wXyPlot, base_height=hXyPlot) save_plot('output/plots/tmerge.vs.SIRVs.stats/ont-Cshl-CapTrap/HpreCap/ont-Cshl-CapTrap_HpreCap_0+_Brain01Rep1.noFilt.tmerge.min2reads.vs.SIRVs.stats.xy.woLegend.png', pXyNoLegend, base_width=wXyNoLegendPlot, base_height=hXyNoLegendPlot) save_plot('output/plots/tmerge.vs.SIRVs.stats/ont-Cshl-CapTrap/HpreCap/ont-Cshl-CapTrap_HpreCap_0+_Brain01Rep1.noFilt.tmerge.min2reads.vs.SIRVs.stats.yx.wLegend.png', pYx, base_width=wYxPlot, base_height=hYxPlot) save_plot('output/plots/tmerge.vs.SIRVs.stats/ont-Cshl-CapTrap/HpreCap/ont-Cshl-CapTrap_HpreCap_0+_Brain01Rep1.noFilt.tmerge.min2reads.vs.SIRVs.stats.yx.woLegend.png', pYxNoLegend, base_width=wYxNoLegendPlot, base_height=hYxNoLegendPlot)