r - 如何在密度图中打破 x 轴

Question

我想绘制一个分布和一个值（带有 abline），它比我的分布中的最小值非常小，所以 abline 不会出现在图中。如何将它们绘制在同一绘图中以操纵 x 轴刻度或插入中断？

data <- rnorm(1000, -3500, 27)
estimate <- -80000

plot(density(data))
abline(v = estimate)

Answer 1

这是一个粗略的解决方案，它不是特别漂亮：

library(plotrix)

d <- density(data)
gap.plot(c(-8000,d$x), c(0,d$y), gap=range(c(-7990,-3620)), 
         gap.axis="x", type="l", xlab="x", ylab="Density", 
         xtics=c(-8000,seq(-3600,-3300,by=100)))
abline(v=-8000, col="red", lwd=2)

Answer 2

不完全清楚需要什么，但这可能是进步：

plot(density(data), xlim=range(c(data, estimate+10) ) )
abline(v = estimate, col='red')

在 package:plottrix 中有断轴绘图功能。