带有facet_grid的geom_smooth和不同的拟合函数

时间:2017-07-04 18:51:58

标签: r ggplot2 facet-grid

首先,为这个例子道歉,但我找不到更好的数据集来证明这个问题。希望它足够了。假设我试图制作一个传输的平面网格(自动与手动)和来自mtcars数据集的齿轮数,其中mpg与位移相对应,如下所示:

# Load library
library(ggplot2)

# Load data
data(mtcars)

# Plot data
p <- ggplot(mtcars,aes(x = disp, y = mpg)) + geom_point() + facet_grid(gear ~ am)
p <- p + geom_smooth()
print(p)

给出,

enter image description here

注意,我使用geom_smooth添加了趋势线,并且默认使用黄土曲线。我可以使用nls为方法拟合用户定义的函数而不是黄土曲线,然后说明公式,这很好。但是,是否可以为每个方面拟合不同的用户指定的曲线?例如,左上方面板的线性回归和右下角的衰减指数。这可能吗?还是我用锤子开螺丝?

修改: 自定义(即用户定义的)拟合函数的解决方案是here

1 个答案:

答案 0 :(得分:3)

根据here给出的建议,可能的解决方案是:

# Load library
library(ggplot2)

# Load data
data(mtcars)

# Vector of smoothing methods for each plot panel
meths <- c("loess","lm","lm","lm","lm","lm","lm")

# Smoothing function with different behaviour in the different plot panels
mysmooth <- function(formula,data,...){
   meth <- eval(parse(text=meths[unique(data$PANEL)]))
   x <- match.call()
   x[[1]] <- meth
   eval.parent(x)
}

# Plot data
p <- ggplot(mtcars,aes(x = disp, y = mpg)) + geom_point() + facet_grid(gear ~ am)
p <- p + geom_smooth(method="mysmooth")
print(p)

enter image description here