R Mr.How知其所以然之主题R代码模块
R Mr.How知其所以然之主题R代码模块
原创文章，版权所有，搬运文章，转发请注明polyaD，原文链接https://polyad.github.io/tags
如有任何问题或疑惑，请在后面留言或者发邮箱留言polyaluthor@gmail.com，加微信polyad或者qq数学算法开发10群 282642152进行探讨,请备注：Mr.How。
—
目录
逻辑大纲
使用R进行绘图
清空缓存，引入包，设置当前目录，清空背景，读入数据
rm(list=ls())
library(ggplot2)
library(RColorBrewer)
setwd("D:\\TempC\\")
blank_theme <- theme_minimal()+
theme(
axis.title.x = element_blank(),
axis.title.y = element_blank(),
panel.border = element_blank(),
panel.grid=element_blank(),
axis.ticks = element_blank(),
plot.title=element_text()
)
data <- read.table("lengthVSnumber_ra1.txt", sep="\t",header=T)
tf<-"Length of reads distribution"
mydata<-data.frame(data)
mydata
fn<-"lengthVSnumber.png"
fd<-"lengthVSnumber.pdf"
pdf(fd, 8, 6)
#colours=c("blue", "purple","red","violet", "yellow", "orange","lightblue")
#png(file=fn, bg="transparent",width=800,height=600)
p <- ggplot(mydata,aes(x=mydata$length,y=mydata$reads/100000,fill=mydata$group))+geom_bar(position="dodge",stat="identity")
p+blank_theme
p+scale_x_continuous(limits=c(15, 26),breaks=seq(15, 26, 1))+
scale_y_continuous(breaks=seq(0, 200, 20))+
xlab("Length (nt)") +
ylab("Number of reads(*100000) ") +
labs( title = tf,fill="types")+
theme(axis.text = element_text( size = 8, hjust = 0, vjust = 0, face = 'bold'))+
theme(plot.title = element_text(hjust = 0.5,vjust = 0.6,size=14, face="bold"))+
theme(legend.text = element_text( size = 10, hjust = 3, vjust = 3, face = 'bold'))
dev.off()
png(file=fn, bg="transparent",width=800,height=600)
p <- ggplot(mydata,aes(x=mydata$length,y=mydata$reads/100000,fill=mydata$group))+geom_bar(position="dodge",stat="identity")
p+blank_theme
p+scale_x_continuous(limits=c(15, 26),breaks=seq(15, 26, 1))+
scale_y_continuous(breaks=seq(0, 200, 20))+
xlab("Length (nt)") +
ylab("Number of reads(*100000) ") +
labs( title = tf,fill="types")+
theme(axis.text = element_text( size = 8, hjust = 0, vjust = 0, face = 'bold'))+
theme(plot.title = element_text(hjust = 0.5,vjust = 0.6,size=14, face="bold"))+
theme(legend.text = element_text( size = 10, hjust = 3, vjust = 3, face = 'bold'))
dev.off()
#参考资料
—–
一级资料文献与书籍及重要作者
文献：
书籍：
博客：
论坛：
视频：
二级资料：他人加工且有观点及大众资料
博客：
论坛：
视频：
返回顶部
评论：
技术文章推送
知其所以然主题模块分享
微信公众号：How先生polyad