comparison makeblobplot.R @ 6:b45fc20c7acd

resize legend and removed facet

5:dfa4f03b5baa

for (t in levels(mfilt[,taxlevel])) {
  g <- g + geom_point(data=mfilt[mfilt[,taxlevel]==t,],aes_string(x="gc", y="cov", colour=taxlevel), size=2, alpha=I(1/3))
}
#y_axis_breaks = c(1,2,5,10,20,50,100,200,500,1000);
g<-g +
  facet_wrap(~read_set, ncol=numcols) + 
  #scale_y_log10() + 
scale_x_continuous(limits=c(1e-307, 1),breaks = seq(0,1,.1)) +
  labs(x="GC content", y="Read coverage") + 
  guides(colour = guide_legend(nrow=3, override.aes = list(alpha = 1,size=10))) + 
  theme (
    strip.text.x = element_text(colour = "black", size = 25, vjust = 0.5),
    axis.text.x  = element_text(colour = "black", size = 25, vjust = 1),
    axis.text.y  = element_text(colour = "black", size = 25, vjust = 0.5),
    axis.title.x = element_text(colour = "black", size = 25, vjust = 0),
    axis.title.y = element_text(colour = "black", size = 25, hjust = 0.5, vjust = 0.5, angle=90),
    legend.text  = element_text(colour = "black", size = 25, vjust = 0),
    legend.title = element_text(colour = "black", size = 25, vjust = 0, hjust = 0, lineheight=1),
    legend.justification=c(1,1), legend.position="bottom", legend.direction="horizontal"
  )
g
dev.off()

6:b45fc20c7acd

for (t in levels(mfilt[,taxlevel])) {
  g <- g + geom_point(data=mfilt[mfilt[,taxlevel]==t,],aes_string(x="gc", y="cov", colour=taxlevel), size=2, alpha=I(1/3))
}
#y_axis_breaks = c(1,2,5,10,20,50,100,200,500,1000);
g<-g +
  #facet_wrap(~read_set, ncol=numcols) + 
  scale_y_log10() + 
scale_x_continuous(limits=c(1e-307, 1),breaks = seq(0,1,.1)) +
  labs(x="GC content", y="Read coverage") + 
  guides(colour = guide_legend(nrow=3, override.aes = list(alpha = 1,size=10))) + 
  theme (
    strip.text.x = element_text(colour = "black", size = 25, vjust = 0.5),
    axis.text.x  = element_text(colour = "black", size = 25, vjust = 1),
    axis.text.y  = element_text(colour = "black", size = 25, vjust = 0.5),
    axis.title.x = element_text(colour = "black", size = 25, vjust = 0),
    axis.title.y = element_text(colour = "black", size = 25, hjust = 0.5, vjust = 0.5, angle=90),
    legend.text  = element_text(colour = "black", size = 15, vjust = 0),
    legend.title = element_text(colour = "black", size = 15, vjust = 0, hjust = 0, lineheight=1),
    legend.justification=c(1,1), legend.position="bottom", legend.direction="horizontal"
  )
g
dev.off()