如何在密度圖上創建x軸斷裂？

Question

我正在嘗試繪制范圍廣泛的數據集的密度圖。

數據= c（-10，-20，-20，-18，-17,1000,10000，500，500，500，500000）

地塊（密度（數據））

正如您在圖中看到的，我們看不到太多 。

有沒有一種方法可以使x軸上的一個（或多個）軸斷裂，以便更好地可視化數據的分布？ 或者，是否有一種方法可以在多個圖形中繪制一定范圍的數據，然后將其粘貼在一起？

非常感謝！

Answer 1

包plotrix有一個函數gap.plot() ，但我認為它存在一些問題（請參閱如何在R中突破y-data-range繪制“多”曲線 ）。 我建議您繪制兩個圖。

  ## use small margins and relatively big outer margins (to write labels).
old.par <- par(mfrow = c(1, 2), mar = rep(0.5, 4), oma = c(4, 4, 1, 1))

plot(density(data), xlim = c(-1000, 29000), main = "", bty="c")               # diff 30000
abline(v = par("usr")[2], lty=2)          # keep the same diff of xlim to avoid misleading
plot(density(data), xlim = c(471000, 501000), main = "", yaxt ="n", bty="]")  # diff 30000
abline(v = par("usr"[1]), lty=2)

par(old.par)