.:
read.table("Fig_chr1_wild_africa.bin")->data1
read.table("Fig_chr2_wild_africa.bin")->data2
read.table("Fig_chr3_wild_africa.bin")->data3
read.table("Fig_chr4_wild_africa.bin")->data4
read.table("Fig_chr5_wild_africa.bin")->data5
read.table("Fig_chr6_wild_africa.bin")->data6
read.table("Fig_chr7_wild_africa.bin")->data7
read.table("Fig_chr8_wild_africa.bin")->data8
read.table("Fig_chr9_wild_africa.bin")->data9
read.table("Fig_chr10_wild_africa.bin")->data10
pdf("Fig_wild_africa.pdf")
plot(data1[,1]/1000,data1[,2],type="l",col="blue",main="LD decay",xlab="Distance(Kb)",
xlim=c(0,300),ylim=c(0,1),ylab=expression(r^{2}),bty="n",lwd=2)
lines(data2[,1]/1000,data2[,2],type="l",col="lightblue")
lines(data3[,1]/1000,data3[,2],type="l",col="aquamarine4")
lines(data4[,1]/1000,data4[,2],type="l",col="gray0")
lines(data5[,1]/1000,data5[,2],type="l",col="blueviolet")
lines(data6[,1]/1000,data6[,2],type="l",col="brown2")
lines(data7[,1]/1000,data7[,2],type="l",col="chartreuse2")
lines(data8[,1]/1000,data8[,2],type="l",col="chocolate1")
lines(data9[,1]/1000,data9[,2],type="l",col="darkgoldenrod")
lines(data10[,1]/1000,data10[,2],type="l",col="darkred")
legend(230,1,c("chr1","chr2","chr3","chr4","chr5","chr6","chr7","chr8","chr9","chr10"), lwd=c(5,5),
col=c("blue","lightblue","aquamarine4","gray0","blueviolet","brown2","chartreuse2","chocolate1","darkgoldenrod","darkred"), y.intersp=1)
dev.off()
.:
这是我自己画的命令 发给你