通过 ggplot 中的多个方面进行一次回归

Question

我试图通过一个有 20 个方面的图形来拟合线性回归（实际上是 9 个）。每次我拟合回归时（使用 geom_smooth，使用 method = lm），它拟合 20 条线，一条穿过每个面，但是我希望每个 ReefSpecies 组合的一条线穿过所有 20 个面。

这是我的图：

Similar Figure

这是我目前所拥有的：

Biomass <- c(20, 10, 5, 4, 5, 7, 8, 22, 13, 13, 15, 18, 2, 5, 7, 10)
Season <- c("Winter", "Spring", "Summer", "Fall")
Year <- c("1", "2", "3", "4")
ReefSpecies <- c("Admiral Ma", "Jaap Mf", "Grecian Ma", "Alligator Mf", "Jaap Mf", "Grecian Ma", "Alligator Mf", "Admiral Ma", "Grecian Ma", "Alligator Mf", "Admiral Ma", "Jaap Mf", "Alligator Mf", "Admiral Ma", "Jaap Mf","Grecian Ma")
Seasonal <- data.frame(Biomass, Season, Year, ReefSpecies)

testp <- ggplot(data = Seasonal, aes(x = Season, y = Biomass, group =        ReefSpecies, fill = ReefSpecies, colour = ReefSpecies))
testp <- testp + geom_point(stat = "identity", position="identity", inherit.aes = TRUE)
testp <- testp + facet_grid(. ~ Year, scales="fixed")
testp <- testp + theme(axis.text.x = element_text(angle = 90))
testp <- testp + theme(panel.margin.x = unit(0, "lines"))
testp <- testp + theme(legend.position = "top")
testp

Answer 1

基于 cmets，您确实不想要 place an identical smooth on each facet of a ggplot（您可以通过 setting the faceting variable to NULL in the smooth 完成。

您做想要的是在所有方面都有一个单一的回归。我认为如果没有一些黑客攻击like that shown here，这是不可能的。你可以试试。

但是，我建议您退后一步考虑一下您为什么要这样做以及平滑意味着什么。也许这意味着刻面不是正确的选择？在这种情况下，您可能会考虑定义一个 Time 变量来解释不同年份的季节并在此基础上进行回归（不包括方面）。

一个示例（使用调整后的数据，因为您的示例数据每年的观察次数不超过一次）：

Year <- sort(rep(Year, 4))
Seasonal <- data.frame(Biomass, Season, Year, ReefSpecies)
Seasonal$Time <- interaction(Season, Year)

ggplot(Seasonal, aes( Time,  Biomass, color=ReefSpecies)) + 
  geom_point() +
  geom_smooth(aes(group=ReefSpecies), method="lm")