4

我想知道是否有一种“直接”的方式将 ggplot facet 面板中的回归线的斜率链接到该面板的背景颜色（即在视觉上将大网格中的正斜率与负斜率分开）。

我了解如何在 GGplots 中添加回归线 - 正如在 R 中使用 qplot 添加回归线到 facet_grid 中所解释的那样

如果您之前已将此信息添加到原始数据框中，我也了解如何更改背景 - 如使用 facet_grid 有条件地更改面板背景中所述？

但是 - 有没有办法在“geom_rect”公式中执行此操作，而不必例如单独运行回归，将它们绑定到原始数​​据帧，然后将其用作 geom_rect() 的变量？geom_rect() 有没有办法使用来自 stat_smooth() 的信息？

沃特

早期问题的简单回归线图的好例子：

library(ggplot2)
x <- rnorm(100)
y <-  + .7*x + rnorm(100)
f1 <- as.factor(c(rep("A",50),rep("B",50)))
f2 <- as.factor(rep(c(rep("C",25),rep("D",25)),2))
df <- data.frame(cbind(x,y))
df$f1 <- f1
df$f2 <- f2

ggplot(df,aes(x=x,y=y))+geom_point()+facet_grid(f1~f2)+stat_smooth(method="lm",se=FALSE)
4

1 回答 1

4

这不完全是一个解决方案，而是一种解决方法。但结果似乎很好。您链接到的两个帖子都有解决方案的每个部分。James 的解决方案在这里告诉您如何从stat_smooth. Joran 的解决方案在这里讲述了如何使用geom_rect来填充背景。

# generating data: Usage of set.seed for reproducibility 
# also I changed the multiplication constant to 0.1 to have 
# at least one negative slope.

require(ggplot2)
set.seed(12)
x <- rnorm(100)
y <-  + .1*x + rnorm(100)
f1 <- as.factor(c(rep("A",50),rep("B",50)))
f2 <- as.factor(rep(c(rep("C",25),rep("D",25)),2))
df <- data.frame(cbind(x,y))
df$f1 <- f1
df$f2 <- f2

# first generate your plot in this manner and run it
# from James' post, the part outfit=fit<<-..y.. will store 
# the output of fitted values in "fit"

g <- ggplot(df,aes(x=x,y=y)) + geom_point()+facet_grid(f1~f2) 
g <- g + stat_smooth(aes(outfit=fit<<-..y..), method="lm",se=FALSE)
# now run g to generate "fit"
g

# now extract the slope for each facet and 
# construct the data.frame for geom_rect (as per Joran's post)
# Edit: Just to add more info about "fit". By default it contains
# 80 values per facet. Hence the 80*4 = 320

slopes <- fit[seq(2, 320, by = 80)] - fit[seq(1, 320, by = 80)]
tp <- unique(df[, c('f1', 'f2')])
tp <- transform(tp, slopes=slopes, x=1, y=1)
tp$pos_neg <- ifelse(slopes > 0, 1, 0)
tp$pos_neg <- factor(tp$pos_neg)

# now plot again (but with geom_rect)
g <- ggplot(df,aes(x=x,y=y)) 
g <- g + geom_rect(data = tp, aes(fill = pos_neg), xmin = -Inf, xmax = Inf, ymin = -Inf, ymax = Inf, alpha = 0.5) 
g <- g + geom_point() + facet_grid(f1~f2) + stat_smooth(method = "lm",se = FALSE)
g

输出看起来像这. 我不确定这是否是您所期望的。严格来说，您确实计算了两次拟合值，但两次您都使用stat_smooth. 就像我说的，它只是一种解决方法。

于 2012-12-27T23:03:04.520 回答