首先,为这个例子道歉,但我找不到更好的数据集来证明这个问题。希望它足够了。假设我试图制作一个传输的平面网格(自动与手动)和来自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)
给出,
注意,我使用geom_smooth
添加了趋势线,并且默认使用黄土曲线。我可以使用nls
为方法拟合用户定义的函数而不是黄土曲线,然后说明公式,这很好。但是,是否可以为每个方面拟合不同的用户指定的曲线?例如,左上方面板的线性回归和右下角的衰减指数。这可能吗?还是我用锤子开螺丝?
修改: 自定义(即用户定义的)拟合函数的解决方案是here。
答案 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)