[TOC]
ggplot2 学习笔记
1.安装
install.packages('ggplot2')
2.qplot
set.seed(1410)##让样本可重复
dsmall<-diamonds[sample(nrow(diamonds),100),]#取前100行
head(dsmall,4)
qplot(carat,price,data=diamonds)
> head(diamonds,4)
carat cut color clarity depth table price x y z
1 0.23 Ideal E SI2 61.5 55 326 3.95 3.98 2.43
2 0.21 Premium E SI1 59.8 61 326 3.89 3.84 2.31
3 0.23 Good E VS1 56.9 65 327 4.05 4.07 2.31
4 0.29 Premium I VS2 62.4 58 334 4.20 4.23 2.63
Rplot01
qplot(carat,price,data=dsmall,colour=color,shape=cut)
Rplot02
qplot(log(carat),log(price),data=diamonds,alpha=I(1/100))
# I()手动设定图形属性，如color=I('red');alpha设定透明度，0为完全透明
Rplot03
#geom="point",散点图
#geom="smooth",拟合平滑曲线
#geom="boxplot",箱线图
#geom="path"和geom="line",可在数据点之间绘制连线
#geom="histogram",直方图
#geom="freqploy",多边形
#geom="density",密度曲线
#geom="bar",条形图
qplot(carat,price,data=dsmall,geom=c("point","smooth")) # 若不想绘制标准误差，se=FALSE
Rplot04
qplot(color,price/carat,data=diamonds,geom="jitter",alpha=I(1/5)) #扰动点图
Rplot05
qplot(carat,data=diamonds,geom="histogram",binwidth=0.05,xlim=c(0,4),fill=color)
Rplot06
qplot(carat,data=diamonds,geom="density",xlim=c(0,4),color=color)
Rplot07
qplot(cut,price,data=diamonds,geom="boxplot",fill=cut)
Rplot09
2.1分面
#row_var ~ col_var;row_var ~ .会创建一个单列多行的图形矩阵
qplot(carat,data=diamonds,facets=color ~ .,geom="histogram",binwidth=0.1,xlim=c(0,3))
Rplot08
2.2选项
# xlim,ylim: x轴和y轴范围，xlim=c(0,20)
# log： log='x',表示x轴取对数
# main：图形主标题
# xlab,ylab: x,y轴标签文字
# 建立图形
p<-qplot(....)
# 保存图形对象
save(p,file="plot.rdata")
# 读入图形对象
load("plot.rdata")
# 将图片保存成png格式
ggsave("plot.png",width=5,height=5)
2.3甲基化的部分图
v=''
v$a=c(10,20,40,60,70)
v$b=c(-0.01,-0.02,-0.01,-0.02,-0.01)
v$d=c('Upstream','TSS','Genebody','TTS','Downstream')
v<-data.frame(v)
sort<-df[240]
sort<-sort(sort)
sort[240]
library('ggplot2')
df<-read.table('element.final',header=TRUE)
svg('TSS.svg')
ggplot(data=df)+geom_line(aes(x=pos,y=value/sort[240],color=type))+
ylab('normalized methylation level')+geom_text(data=v,aes(x=a,y=b,label=d),size=3,color=I('red'))+
xlab('')
dev.off()
TSS
setwd('/lustre/project/og04/shichunwei/project/bisulfite_sequencing/mouse/out/SRR1040649/test')
df<-read.table("NC_000084.6.txt",header=TRUE)
library('ggplot2')
qplot(count,value,data=df,geom="line",color=color,facets=type~.,alpha=I(8/10))+theme(legend.position='none')+
xlab('')+ylab('')
chr
2.4 其他选项
library('ggplot2')
p<-qplot(cyl,wt,data=mtcars)
p+scale_x_continuous(breaks=c(5.5,6.5)) # breaks 横坐标只显示5.5，6.5
p+scale_x_continuous(limits=c(5.5,6.5)) # limits 只显示5.5~6.5范围
p<-qplot(wt,cyl,data=mtcars,color=cyl)
p+scale_color_gradient(breaks=c(5.5,6.5)) # 与上面的类似，只是作用范围不是横坐标而是图注
p+scale_color_gradient(limits=c(5.5,6.5)) #
scale_x_log10() # x 轴转换成log10
#图例
theme(legend.position=right,left,top,bottom,none) # 控制图例位置
Rplot10
pie<-ggplot(mtcars,aes(x=factor(1),fill=factor(cyl)))+geom_bar(width=1)
pie + coord_polar(theta='y')+theme_bw()+xlab('')+ylab('')+labs(fill='')+
theme(axis.text = element_blank())+theme(axis.ticks = element_blank())
# theme_bw() 更改主题背景，改为白色
Rplot11