先放图
因为图要放在文章里,所以先放一部分图展示
![](https://img.haomeiwen.com/i20853302/6f3dd5a76d5d99a8.png)
主要用到下面的代码,以及ggplot2
cowplot
软件包。后期用AI处理一下
setwd("D:/迅雷下载/师兄作图/G值/")
data<-read.table("all2.txt",header = T,sep="\t")
library(doBy)
library(plyr)
library(ggplot2)
library(cowplot)
g1<-subset(data,data$chr==1)
g2<-subset(data,data$chr==2)
g3<-subset(data,data$chr==3)
g4<-subset(data,data$chr==4)
g5<-subset(data,data$chr==5)
g6<-subset(data,data$chr==6)
g7<-subset(data,data$chr==7)
g8<-subset(data,data$chr==8)
g9<-subset(data,data$chr==9)
g10<-subset(data,data$chr==10)
g11<-subset(data,data$chr==11)
g12<-subset(data,data$chr==12)
g13<-subset(data,data$chr==13)
g14<-subset(data,data$chr==14)
g15<-subset(data,data$chr==15)
g16<-subset(data,data$chr==16)
g17<-subset(data,data$chr==17)
g18<-subset(data,data$chr==18)
g19<-subset(data,data$chr==19)
g20<-subset(data,data$chr==20)
g21<-subset(data,data$chr==21)
g22<-subset(data,data$chr==22)
g23<-subset(data,data$chr==23)
g24<-subset(data,data$chr==24)
g25<-subset(data,data$chr==25)
g26<-subset(data,data$chr==26)
g27<-subset(data,data$chr==27)
g28<-subset(data,data$chr==28)
g29<-subset(data,data$chr==29)
g30<-subset(data,data$chr==30)
g31<-subset(data,data$chr==31)
G1<-ggplot()+
geom_point(aes(pos,G), data = g1,color="#ff7b42",pch=20,show.legend = FALSE, alpha=0.8,size=1) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g1,span = 1,col="#000000",size=3,se=FALSE)+ theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("1")
G2<-ggplot()+
geom_point(aes(pos,G), data = g2,color="#7fff00",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g2,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("2")
G3<-ggplot()+
geom_point(aes(pos,G), data = g3,color="#6495ed",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g3,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("3")
G4<-ggplot()+
geom_point(aes(pos,G), data = g4,color="#f000f0",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g4,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("4")
G5<-ggplot()+
geom_point(aes(pos,G), data = g5,color="#ee0065",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g5,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("5")
G6<-ggplot()+
geom_point(aes(pos,G), data = g6,color="#c2ff4c",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g6,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("6")
G7<-ggplot()+
geom_point(aes(pos,G), data = g7,color="#24c3ff",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g7,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("7")
G8<-ggplot()+
geom_point(aes(pos,G), data = g8,color="#ff7b42",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g8,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("8")
G9<-ggplot()+
geom_point(aes(pos,G), data = g9,color="#7fff00",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g9,span = 1,col="#000000",size=3,se=FALSE)+ theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("9")
G10<-ggplot()+
geom_point(aes(pos,G), data = g10,color="#6495ed",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g10,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+ theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("10")
G11<-ggplot()+
geom_point(aes(pos,G), data = g11,color="#f000f0",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g11,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("11")
G12<-ggplot()+
geom_point(aes(pos,G), data = g12,color="#ee0065",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g12,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("12")
G13<-ggplot()+
geom_point(aes(pos,G), data = g13,color="#c2ff4c",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g13,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("13")
G14<-ggplot()+
geom_point(aes(pos,G), data = g14,color="#24c3ff",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g14,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("14")
G15<-ggplot()+
geom_point(aes(pos,G), data = g15,color="#ff7b42",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g15,span = 1,col="#000000",size=3,se=FALSE)+ theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+ theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("15")
G16<-ggplot()+
geom_point(aes(pos,G), data = g16,color="#7fff00",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g16,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("16")
G17<-ggplot()+
geom_point(aes(pos,G), data = g17,color="#6495ed",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g17,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("17")
G18<-ggplot()+
geom_point(aes(pos,G), data = g18,color="#f000f0",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g18,span = 1,col="#000000",size=3,se=FALSE)+ theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("18")
G19<-ggplot()+
geom_point(aes(pos,G), data = g19,color="#ee0065",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g19,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("19")
G20<-ggplot()+
geom_point(aes(pos,G), data = g20,color="#c2ff4c",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g20,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("20")
G21<-ggplot()+
geom_point(aes(pos,G), data = g21,color="#24c3ff",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g21,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("21")
G22<-ggplot()+
geom_point(aes(pos,G), data = g22,color="#ff7b42",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g22,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("22")
G23<-ggplot()+
geom_point(aes(pos,G), data = g23,color="#7fff00",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g23,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("23")
G24<-ggplot()+
geom_point(aes(pos,G), data = g24,color="#6495ed",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g24,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("24")
G25<-ggplot()+
geom_point(aes(pos,G), data = g25,color="#f000f0",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g25,span = 1,col="#000000",size=3,se=FALSE)+ theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("25")
G26<-ggplot()+
geom_point(aes(pos,G), data = g26,color="#ee0065",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g26,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("26")
G27<-ggplot()+
geom_point(aes(pos,G), data = g27,color="#c2ff4c",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g27,span = 1,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("27")
G28<-ggplot()+
geom_point(aes(pos,G), data = g28,color="#24c3ff",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g28,span = 0.2,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("28")
G29<-ggplot()+
geom_point(aes(pos,G), data = g29,color="#ff7b42",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g29,span = 0.2,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("33")
G30<-ggplot()+
geom_point(aes(pos,G), data = g30,color="#7fff00",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g30,span = 0.2,col="#000000",size=3,se=FALSE)+ theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+xlab("LGE64")
G31<-ggplot()+
geom_point(aes(pos,G), data = g31,color="#6495ed",pch=20,show.legend = FALSE, alpha=0.8,size=5) +
geom_smooth(mapping=aes(x=pos,y=GG),data=g31,span =0.2,col="#000000",size=3,se=FALSE) + theme_bw()+
scale_y_continuous(breaks =seq(0,80,10),limits=c(0,80))+
theme(axis.title.x=element_text(size=39),axis.text.x=element_blank(),axis.ticks.x=element_blank(),axis.title.y=element_blank(),axis.text.y=element_blank(),axis.ticks.y=element_blank())+ xlab("Z")
p<-plot_grid(G1, G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G14,G15,G16,G17,G18,G19,G20,G21,G22,G23,G24,G25,G26,G27,G28,G29,G31,G30,align = "hv",ncol=31)
ggsave(p,filename = "chr33zt.png",width = 78,height =24,limitsize = FALSE )
网友评论