2017-10-17 48 views
0

我正尝试在辅助轴上创建一个点阵图。然而,一旦我绘制第二个点阵图(带有par(new = T)),我无法弄清楚如何不在轴侧= 1中显示轴滴答。这是我的代码用模拟数据:带辅助轴的点阵图

y1_i <- c(2,8,2,14,2) 
y2_i <- c(15,17,28,22,30) 
y1_f <- c(4,9,11,16,7) 
y2_f <- c(13,11,16,11,21) 

y=c(y1_i,y2_i,y1_f,y2_f) 

x <- c("AAEG","AALO","AGAM","ACHR","AALB") 

y1=c(y1_i,y1_f) 
y2=c(y2_i,y2_f) 

dotchart(y1_i,labels=x,xlab="N50 length",xlim = c(0,max(y1))) 
par(new=T) 
dotchart(y2_i,labels=x,xlim = c(0,max(y2))) 
axis(side=3) 

enter image description here

此外,如果可能的话,我想补充这将被稍微推压所述第一数据集的垂直上方(到不重叠的话)的第二数据集,但仍然对应于相同的y轴类别。

谢谢你的任何建议:)

回答

0

发现,通过使用dotchart2从Hmisc包

library(Hmisc) 

y1_i <- c(2,8,2,14,2) 
y2_i <- c(15,17,28,22,30) 
y1_f <- c(4,9,11,16,7) 
y2_f <- c(13,11,16,11,21) 

y=c(y1_i,y2_i,y1_f,y2_f) 

x <- c("AAEG","AALO","AGAM","ACHR","AALB") 

y1=c(y1_i,y1_f) 
y2=c(y2_i,y2_f) 
y1_i <- c(2,8,2,14,2) 
y2_i <- c(15,17,28,22,30) 
y1_f <- c(4,9,11,16,7) 
y2_f <- c(13,11,16,11,21) 

y=c(y1_i,y2_i,y1_f,y2_f) 

x <- c("AAEG","AALO","AGAM","ACHR","AALB") 

y1=c(y1_i,y1_f) 
y2=c(y2_i,y2_f) 

dotchart2(y1_i,labels=x,xlab="N50 length",xlim = c(0,max(y1))) 
par(new=T) 
dotchart2(y2_i,labels=x,xlim = c(0,max(y2)),xlab="Scaffold number",lines=F,xaxis=F) 
axis(side=3,xlab="Scaffold number")