首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >R编程-实例

R编程-实例

作者头像
生信补给站
发布2020-08-06 10:33:09
8700
发布2020-08-06 10:33:09
举报
文章被收录于专栏:生信补给站生信补给站

data <- read.csv("11.csv",header=T)

data <- as.matrix(data)

pdf("Number_of_cock_SNPs(M).pdf")

cols =c("royalblue2","yellow2")

par(fig=c(0.1,1,0,0.8))

barplot(data,col = cols,las=2,xlim=c(0,10),sub="Number of cock_SNPs(M)",horiz =TRUE,axis.lty = 1)

j=2.1

for (i in 2:14){

text(4,j,data[1,i],cex=0.9)

text(5,j,paste("(",round(data[1,i]/(data[1,i]+data[2,i]),2),"%",")",sep=""),cex=0.9)

j=j+1.17

}

j=2.1

par(fig=c(0.1,1,0.55,1),new=TRUE)

plot(data,axes=FALSE,col="white",ann=FALSE)

legend("bottom",cex=0.8,pch=c(15,15),col=c("royalblue2","yellow2"),

legend=c("Overlapped SNPs with chicken dbSNPs","Novel SNPs"))

dev.off()

本文参与 腾讯云自媒体分享计划,分享自微信公众号。
原始发表:2016-05-17,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 生信补给站 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档