バープロット上のX軸の両側にプロットするR

Question

私は現在、バープロットのX軸の両側にプロットしようとしているので、さまざまなY点で2つの異なるベクトルを簡単に比較できます

私はそれの大部分を理解しました。私はちょうど第2のベクトルを追加できません。私が望むものをより明確にする最も簡単な方法は、以下のイメージを示すことです。

私が望むのは、cleanvector2をX軸から「反転」してプロットすることです.0から始まり、下に向かって拡大します。好ましくはY軸上の正の値でもよい。現在使用されて

やコード：

naam=" " 
    filenaam=paste(t, ".png", sep="") 
    bitmap(filenaam, type="png16m", res=300) 
    cleanvector1=c(1.29 1.93 1.71 1.35 1.54 0.73 1.41 1.52 1.43 1.09 0.96 0.87 2.29 1.50 0.67) 
    cleanvector2=c(1.29 1.93 1.71 1.35 1.54 0.73 1.41 1.52 1.43 1.09 0.96 0.87 2.29 1.50 0.67) 
    cleanvector2=cleanvector2*-1 
    barplot(rollapply(cleanvector1, 2,mean, by=1, fill=NA), main=naam, ylim=c(max(cleanvector1)*-1,max(cleanvector1)), xlab="200 basepair bins") 
# barplot(rollapply(cleanvector2, 2,mean, by=1, fill=NA)) 
    abline(h=1, col="darkred") 
    abline(h=-1, col="darkred") 
    axis(side=1) 
    dev.off() 

Answer 1

ここでは、最小限の作業例です。注意すべき2つの点：コードにcleanvector*にカンマがないため、人々が結果を再現するのが難しくなります。たとえば、次のような出力をコピーして貼り付けてください。 dput(cleanvector1)。第二に、rollapplyがzooパッケージから来ているので、この例ではパッケージをロードしているか、zoo::rollapplyを使用しています。

barplot(..., add = TRUE)を使用し、barplotに第二のベクターを追加するには：

cleanvector1 <- c(1.29, 1.93, 1.71, 1.35, 1.54, 0.73, 1.41, 1.52, 1.43, 1.09, 
0.96, 0.87, 2.29, 1.5, 0.67) 
cleanvector2 <- cleanvector1 * -1 
barplot(zoo::rollapply(cleanvector1, 2, mean, fill = NA), 
     ylim=c(max(cleanvector1) * -1, max(cleanvector1)), 
     xlab="200 basepair bins") 
barplot(zoo::rollapply(cleanvector2, 2, mean, fill = NA), 
     add = TRUE) 
abline(h = 1, col = "darkred") 
abline(h = -1, col = "darkred")