r语言 - ggplot2 syntax for geom_smooth

Question

我有以下函数，我用它来拟合一组数据的非线性曲线。我一直收到以下错误:

Error in do.call("layer", list(mapping = mapping, data = data, stat = stat,  : 
 object 'y' not found 

错误来自于geom_smooth拟合。你能给点建议吗?

下面是一些测试数据和函数:

library("ggplot2") 
DF <- data.frame(RFP_fold=1:20, E=rnorm(20)) 

flimPlot <- function(data) {
    ggplot(data, aes(x=RFP_fold, y=E)) +
    geom_point(shape=1) +
    geom_smooth(method=nls,                     # Add non linear regression fit
            formula='y ~ a * (x / (x + K))',    # Forumla for fit
            start=list(a=max(y), K = 0.1),  # set the parameters
            se = FALSE,         # No shaded CI
            fullrange=TRUE)    # Extend regression line
   }
flimPlot(DF)
# Error in do.call("layer", list(mapping = mapping, data = data, stat = stat,  : 
object 'y' not found 

Answer 1

我不知道为什么变量没有被传递给geom_smooth，所以我只是在函数中添加了更多的变量。

flimPlot <- function(data) {
    y <- data$E
    x <- data$RFP_fold
    a <- max(dat$E)
    st <- list(a=a, K=0.1)
    fit <- nls(formula="E ~ a * (RFP_fold / (RFP_fold + K))", dat, start=st )
    K <- round(coef(fit))
    ggplot(data, aes(x=RFP_fold, y=E)) +
    geom_point(shape=1) +
    #geom_text(aes(3,15, label = paste(expr, K, sep=""))) +
    #annotate("text", x = 3, y = 15, label = paste(expr, K, sep="") ) +
    stat_smooth( geom = "smooth",
                 method=nls,   # Add non linear regression fit
                 formula="y ~ a * (x / (x + K))",    # Forumla for fit
                 start=st,  # set the parameters
                 #se = FALSE,         # No shaded CI
                 fullrange=TRUE) +   # Extend regression line
    scale_y_continuous(limits = c(0, 20)) +
    theme(plot.background = element_blank(),
          panel.grid.major = element_blank(),
          panel.grid.minor = element_blank(),
          panel.border = element_blank(), 
          panel.background = element_blank(), 
          axis.line = element_line(color = 'black'),
          text = element_text(size=12),
          #axis.title.x = element_blank(),
          #axis.title.y = element_blank(),
          legend.justification=c(1,0)) +
    ylab("FRET E%") +
    xlab("acceptor intensity")
   }