Question

有一个很好的解释here如何使用ggplot2创建散点图，使用nls拟合数据，并绘制拟合，所有这些都在一行中，如此

myhist = data.frame(size = 10:27, counts = c(1L, 3L, 5L, 6L, 9L, 14L, 13L, 23L, 31L, 40L, 42L, 22L, 14L, 7L, 4L, 2L, 2L, 1L) )

ggplot(data=myhist, aes(x=size, y=counts)) + geom_point() + 
    geom_smooth(method="nls", formula = y ~ N * dnorm(x, m, s), se=F, 
    start=list(m=20, s=5, N=300))

我的问题是：使用这种结构，是否可以从该调用中取出实际的nls对象？我想知道我的系数等。现在我无法弄清楚如何在不进行单独的nls调用的情况下获取它们。

Answer 1

我的问题是：使用这种结构，是否可以从该调用中取出实际的nls对象？我想知道我的系数等等。

目前ggplot2无法实现。 ggplot2函数从模型返回预测，但不返回模型对象本身。因此，您无法从nls对象中提取ggplot对象以查找系数等。

ggplot2和ggplot2-dev邮件列表中有两个相关的讨论：

https://groups.google.com/d/topic/ggplot2/7tiUB2sjCxM/discussion

https://groups.google.com/d/topic/ggplot2-dev/dLGJnzIg4ko/discussion

快速概要：

虽然许多用户要求从ggplot对象中提取统计信息，但开发人员正在考虑它，但似乎有点反对。他们希望用户使用ggplot2进行可视化，并使用适当的建模函数来探索建模参数。但是，Hadley支持实现将模型对象传递给ggplot()调用的能力的想法。因此，您不会尝试从nls对象中提取ggplot对象，而是：

mod <- nls(y ~ N * dnorm(x, m, s), se = F, start = list(m = 20, s = 5, N = 300), 
        data = myhist)
ggplot(data = myhist, aes(x = size, y = counts)) + geom_point() + 
        geom_smooth(mod)

这样，模型只需要调用一次，你可以做任何你想做的事情，而且你不必去搜索ggplot个对象来找到它。但是，我不知道何时或是否会实施。

从ggplot2创建的nls拟合中提取系数

1 个答案: