使用ggplot2拟合nls - 类型为“symbol”的错误对象不是子表

时间:2014-12-05 11:35:46

标签: r ggplot2 curve-fitting data-fitting nls

我正在尝试使用以下数据与ggplot配合:     df <- data.frame(t = 0:30, m = c(125.000000, 100.248858, 70.000000, 83.470795, 100.000000, 65.907870, 66.533715, 53.588087, 61.332351, 43.927435, 40.295448, 44.713459, 32.533143, 36.640336, 40.154711, 23.080295, 39.867928, 22.849786, 35.014645, 17.977267, 21.159180, 27.998273, 21.885735, 14.273962, 13.665969, 11.816435, 25.189016, 8.195644, 17.191337, 24.283354, 17.722776)

我到目前为止的代码(有点简化)是

ggplot(df, aes(x = t, y = m)) +
geom_point() +
geom_smooth(method = "nls", formula=log(y)~x)

但是我收到以下错误

Error in cll[[1L]] : object of type 'symbol' is not subsettable

我浏览过stackoverflow并发现了类似的问题,但我无法解决问题。我真的想在不改变轴的情况下绘制数据。

非常感谢任何帮助。

1 个答案:

答案 0 :(得分:5)

对于nls,您必须更仔细地指定参数。此外,您可能不想使用log(y),因为这将绘制对数而不是y。我的猜测是你想使用像y ~ exp(a + b * x)这样的东西。请参阅下面的示例。

ggplot(df, aes(x = t, y = m))+
  geom_point()+
  geom_smooth(method = "nls", formula = y  ~ exp(a + b * x), start=list(a=0, b=1), se=FALSE)

您也可以使用glm代替nls。例如,以下内容为您提供相同的解决方案。

ggplot(df, aes(x = t, y = m))+
  geom_point()+
  geom_smooth(method = "glm", formula = y  ~ x, family=gaussian(link = "log"), se=FALSE)