CODEWITHSUNDEEP
rplotggplot2facet
Andrew Robinson
Andrew Robinson

Reputation: 141

How do I place an identical smooth on each facet of a ggplot2 object?

Here's an example: 

eg <- data.frame(x = c(1:50, 50:1),  
                 y = c(1:50, 1:50) + rnorm(100),  
                 g = rep(c("a","b"), each=50))  

qplot(x, y, data = eg) +  
  facet_wrap(~ g) +  
  geom_smooth()

I'd like to be able to plot the overall smooth on both facets as well as having the facet-specific smooths.

Edit: here's one way.

my.smooth <- gam(y ~ s(x), data = eg)
my.data <- data.frame(x = 1:50)                                           
my.data$y <- predict(my.smooth, newdata = my.data) 

qplot(x, y, data = eg) + 
    facet_wrap(~ g) + 
    geom_smooth() + 
    geom_smooth(data = my.data)

Thanks for any help!

Andrew

Upvotes: 4

Views: 2302

Answers (1)

Jon Olav Vik
Jon Olav Vik

Reputation: 1461

Clever trick: setting the faceting variable to NULL

library(ggplot2)
eg <- data.frame(x = c(1:50, 50:1),  
                 y = c(1:50, 1:50) + rnorm(100),  
                 g = rep(c("a","b"), each=50))  

p <- qplot(x, y, data = eg) +  
  facet_wrap(~ g) +  
  geom_smooth()

p + geom_smooth(data=within(eg, g <- NULL), fill="red")

Or if you prefer, use facet_grid(..., margins=TRUE):

p + facet_grid(.~g, margins=TRUE)

Upvotes: 6

Related Questions