带有 ggplot 的残差图，X 轴为“排序”残差

Question

我正在尝试在 ggplot: 中重新创建这样的情节。

此图从回归输出中获取残差，并按顺序绘制它们（X 轴是残差的等级）。

我对此的最佳尝试如下：

library(ggplot2)
library(modelr)

d <- d %>% add_residuals(mod1, var = "resid")
d$resid_rank <- rank(d$resid)

ggplot(data = d, aes(x = resid_rank, y = resid)) +
  geom_bar(stat="identity") +
  theme_bw()

但是，这会产生一个完全空白的图表。我试过这样的事情：

ggplot(data = d, aes(x = resid_rank, y = resid)) +
  geom_segment(yend = 0, aes(xend=resid)) +
  theme_bw()

但这会产生方向错误的细分。这样做的正确方法是什么，并通过第三个因素为这些线条着色？

假数据集：

library(estimatr)
library(fabricatr)

#simulation
dat <- fabricate(
  N = 10000,
  y =  runif(N, 0, 10),
  x = runif(N, 0, 100)
)

#add an outlier
dat <- rbind(dat, c(300, 5))
dat <- rbind(dat, c(500, 3))

dat$y_log <- log(dat$y)
dat$x_log <- log(dat$x)
dat$y_log_s <- scale(log(dat$y))
dat$x_log_s <- scale(log(dat$x))

mod1 <- lm(y_log ~ x_log, data = dat))

Answer 1

我使用来自lm() 帮助页面的内置数据集来创建此示例。我也只是直接使用resid() 来获取残差。目前还不清楚彩色条在哪里/为什么会有所不同，但基本上您需要在 data.frame 中添加一列，说明它们为什么是红色或蓝色，然后将其传递给fill。

library(ggplot2)
#> Warning: package 'ggplot2' was built under R version 3.4.4
#example from lm
ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
weight <- c(ctl, trt)
lm.D9 <- lm(weight ~ group)
resids <- data.frame(resid = resid(lm.D9))
#why are some bars red and some blue? No clue - so I'll pick randomly
resids$group <- sample(c("group 1", "group 2"), nrow(resids), replace = TRUE)
#rank
resids$rank <- rank(-1 * resids$resid)


ggplot(resids, aes(rank, resid, fill = group)) +
  geom_bar(stat = "identity", width = 1) +
  geom_hline(yintercept = c(-1,1), colour = "darkgray", linetype = 2) +
  geom_hline(yintercept = c(-2,2), colour = "lightgray", linetype = 1) +
  theme_bw() +
  theme(panel.grid = element_blank()) +
  scale_fill_manual(values = c("group 1" = "red", "group 2" = "blue"))

^{由reprex package (v0.2.1) 于 2019-01-24 创建}