|||
script如下:
data <- read.table("chrome_bs.csv",header=T,sep=",")
#csv文件有标题行,至少含两列数据染色体长度(size)和染色体遗传力(h2)
library("ggplot2")
tiff("chrome_bs.tiff",height=1600,width=2400,res=300,compression="lzw")
ggplot(data,aes(gene,bs))+geom_point(cex=10,pch=21,bg="green",col="red")+geom_smooth(method=lm,colour="blue",size=1.5)+geom_text(aes(label=chr))+xlab("Chromosome genes")+ylab("Chromosomewise heritability")+theme(axis.text=element_text(size=12,face="bold"),axis.line = element_line(size =1.2,colour="black"),axis.title=element_text(face="bold",size=16))+scale_x_continuous(breaks=seq(0,3000,by=400))+ggtitle("Bivariate Stream")
dev.off()
结果类似下图:
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-19 21:33
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社