ggplot2 tips

General usages, click here to see list of opts
  1. theme_get() will show you the “hidden” options that you can use in theme()
  2. ?opts or ?themes see lists of opetions
  3. args(function_name) will list the arguments for this function
  4. The alpha channel is normally used as an opacity channel. If a pixel has a value of 0% in its alpha channel, it is fully transparent (and, thus, invisible), whereas a value of 100% in the alpha channel gives a fully opaque pixel (traditional digital images).
Common variable
  • Two dots variable

The two dots are a visual indicator highlighting that variable is not present in the original data, but has been computed by the statistic.

  • ..density..: bin_count / sum(count)
  • ..scaled..: bin_count / max(count)
  • ..ndensity..: bin_density / max(abs(density))
  • ..ncount..: bin_count / max(abs(count))
  • count: number of points in bin
  • density: density of points in bin, scaled to integrate to 1
  • ncount: count, scaled to maximum of 1
  • ndensity: density, scaled to maximum of 1
#chang axis size
theme(axis.text.x=element_text(size=X)) + theme(axis.text.y=element_text(size=X))
theme(axis.text.x=element_text(angle=90,hjust=1))

#Erase labels and ticks on x,y axis  
p + theme(axis.ticks=element\_blank(), axis.text.x=element\_blank(), \  
axis.text.y=element\_blank(),axis.ticks.x = element\_blank())  
Add Title, xlab, ylab
#wrap the title, you can use "\n" to move the remaining text to a new line:
#ggplot2 doesn't have "subtitle" functionality. 
#But you can use the \n term in any of the labels to drop down a line.
theme(title="text \n more text")
xlab(NULL) + ylab(NULL)
ggtitle("strings")
# A scatterplot with regular (linear) axis scaling
sp <- ggplot(dat, aes(xval, yval)) + geom_point()
sp

# log2 scaling of the y axis (with visually-equal spacing)
library(scales) # Need the scales package
sp + scale_y_continuous(trans=log2_trans())

# log2 coordinate transformation (with visually-diminishing spacing)
sp + coord_trans(y="log2")

#scale_y_log2() will do the transformation first 
#and then calculate the geoms
#coord_trans() will do the opposite: calculate the geoms first, 
#and then transform the axis.
#So sometimes you need coord_trans(ytrans = "log2") instead of #
#scale_y_log2() to avoid data loss.
Options for legends
theme(legend.key.width=unit(1, "in"),
legend.text = theme_text(size=30),
legend.title=element_blank(), #no legend title
legend.key=element_blank(), #no border for legend
legend.position="none"  #"right","left","top"
#0.08 means right away from y-axis, 
#0.8 means above from x-axis, 
#relative to the size of picture
legend.position=c(0.08,0.8) 
legend.direction = "vertical",
legend.justification = "center"
)
Set the order of legends or other variables
#Default, ggplot2 uses alphabetical order. #
#One can change it by given vectors
foomelt$COG <- factor(foomelt$COG, levels = c("first","second",...,"last"), ordered=T)
foomelt$COG <- factor(foomelt$COG, levels = c("first","second",...,"last"))
Use facets to set the layout of pictures, also check here
#ncol=6 means horizontally, six pics one row.
#nrow=6 means vertically, six pics one column.
#scale='free' means each pic can have different axis ranges.
#each pic can also have different x-axis by 'free_x' or y-axis by 'free-y'
facet_wrap(~Size,  ncol=6,  scale='free')   
facet_grid(. ~level, nrow=6, scale="fixed")
facet_grid(vertical_level ~ horizontal_level)
# Calculate correlation coefficient
with(mtcars,cor(wt, mpg, use = "everything", method = "pearson"))
[1] -0.8676594
#annotate the plot
+ geom_abline(intercept = 37, slope = -5) + 
geom_text(data = data.frame(), aes(4.5, 30, label = "Pearson-R = -.87"))
Remove grid line and use white background
theme_bw()
theme(panel.grid.major = element_blank(), 
panel.grid.minor = element_blank())
#theme_blank for old version
[ggplot2 layout] (https://ggplot2-dev.googlegroups.com/attach/5aa16afece3d5bc6/theme0.html?gda=9XgBzEYAAABCncUW0npTUN_veVgl3inYi0oNsf4Sjxsz8g3AimkTHy2Q5nwgitdzQrQMmMK7aytx40jamwa1UURqDcgHarKEE-Ea7GxYMt0t6nY0uV5FIQ&view=1&part=4)
geom_boxplot, also check oneline boxplot
1.Hidden outliers
geom_boxplot(outlier.colour='NA')
2.Adjust ylim
stats <- boxplot.stats(value)$stats
ylim_zoomin <- c(stats[1]/2,stats[5]*2)
p + coord_cartesian(ylim=ylim_zoomin)
Manually set line type and line color, ref
ggplot(mort3, aes(x = year, y = BCmort, col = State, linetype = State)) +
  geom_line(lwd = 1) +
  scale_linetype_manual(values = c(rep("solid", 10), rep("dashed", 6))) +
  scale_color_manual(values = c(brewer.pal(10, "Set3"), brewer.pal(6, "Set3"))) +
  opts(title = "BC mortality") +
  theme_bw()

scale_color_manual(values = c("red",'green','blue')
scale_color_manual(values = c(rgb(255/255,0/255,0/255),
  rgb(0/255,255/255,0/255),rgb(0/255,0/255,255/255))
Manually set ytics and xtics
scale_x_continuous(breaks=round(seq(min(dat$x), mx(dat$x), by=0.5),1))

sclale_y_continuous(breaks=c(8,16,100,128,512,1000))  #any number
Color usage

[color bars] (http://www.colbyimaging.com/wiki/statistics/color-bars)

color plate

CHENTONG
版权声明:本文为博主原创文章,转载请注明出处。
alipay.png WeChatPay.png

CHENTONG

CHENTONG
积微,月不胜日,时不胜月,岁不胜时。凡人好敖慢小事,大事至,然后兴之务之。如是,则常不胜夫敦比于小事者矣!何也?小事之至也数,其悬日也博,其为积也大。大事之至也希,其悬日也浅,其为积也小。故善日者王,善时者霸,补漏者危,大荒者亡!故,王者敬日,霸者敬时,仅存之国危而后戚之。亡国至亡而后知亡,至死而后知死,亡国之祸败,不可胜悔也。霸者之善著也,可以时托也。王者之功名,不可胜日志也。财物货宝以大为重,政教功名者反是,能积微者速成。诗曰:德如毛,民鲜能克举之。此之谓也。

生信宝典文章集锦

### 程序学习心得* [生物信息之程序学习](http://mp.weixin.qq.com/s?__biz=MzI5MTcwNjA4NQ==&mid=2247483927&idx=1&sn=23adf2b9d13400f2081f790e674e...… Continue reading

生信宝典Linux学习系列文章整理

Published on October 18, 2017

在线画图

Published on October 17, 2017