具有复杂功能的stat_smooth

Question

我想在ggplot中使用非线性回归，这样我可以绘制数据并拟合曲线到反应速率r =（k1 * PA * PB ^ 0.5）/（1 + k2 * PA ^ 0.5）其中PA和PB按反应速率r列表。我想找到常量k1和k2。我已经能够使用nls（）进行回归，但我需要在ggplot数据之上绘制曲线。

1. r      PA      PB
2. 0.02   0.002   0.005
3. 0.05   0.004   0.009
4. 0.09   0.005   0.019

plot <- qplot(data = melt.data[which(melt.data$variable ==     c('PA','PB')),], 
              x = Experimental.rate,
              y = value, 
              color = variable)

plot <- plot + stat_smooth(method = nls, 
              formula = with(data, 'y ~ k1 * PA * PB^0.5 / (1 + k2 * PA^0.5 )'), 
              start = list(k1 = 0.1, k2 = 0.1))

这种类型的东西......例如，数据只有30行连续排列。让我知道我还能提供什么来使这个更容易理解 -