如何在密度图中打破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')

在包中：plotrix有断轴绘图功能。