2017-10-17 50 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")