ggplot2——自动放大geom_smooth（使用coord_cartesian）

Question

geom_smooth 很棒，很大程度上是因为它平均了很多变化。然而，正因为如此，当它缩小时，很难看出它在 x 轴上的变化。我正在生成大约 1000 个图表，我需要通过 coord_cartesian 放大 ggplot2。然而，每个图表都有不同的缩放限制。有没有办法让 ggplot2 放大以适应平滑？我对仅放大 geom_smooth 线和 geom_smooth 线加 SE 阴影区域的解决方案感兴趣。

例如，我有兴趣知道如何将这个：

ggplot(data=mtcars, aes(y=qsec,x=wt)) + geom_point() + geom_smooth()

变成这样的：

ggplot(data=mtcars, aes(y=qsec,x=wt)) + geom_point() + geom_smooth() + coord_cartesian(ylim = c(15,20))

而不明确指定限制。

Answer 1

手动拟合平滑模型可以让您更加灵活地实现这种自定义和其他类型的自定义。对于大多数项目，我一开始使用内联界面，但当我需要其他调整时，通常最终不得不切换到手动计算。

另请参阅 Hadley 书中的 §9.1.1。

require(ggplot2)

# Fit smooth manually
fit  = loess(qsec ~ wt, data=mtcars)
newx = data.frame(wt=with(mtcars, seq(min(wt), max(wt), len=100)))
pred = predict(fit, newdata=newx, se=T)
pred = cbind(wt=newx, qsec=pred$fit, se=pred$se.fit)

# Calculate limits based on extent of smooth geom
ylims = with(pred, c(floor(min(qsec-se)), ceiling(max(qsec+se))))

# Plot
dev.new(width=5, height=4)
ggplot(data=mtcars, aes(y=qsec, x=wt)) + 
  geom_point() +
  geom_smooth(aes(ymax=qsec+se, ymin=qsec-se), data=pred, stat='identity') +
  coord_cartesian(ylim = ylims)

但是，这仍然不适用于构面，因为您只能指定，例如 scales='free'，而不是直接的实际限制。