问题标签 [nls]

我正在尝试重现以下代码（nls 函数表现不佳），但通过额外的实现，使用 for 循环、sprintf 和 as.formula()，根据给定光谱中的峰值数量添加变量。为了在峰之间更加连贯，我对每个峰的变量名称进行了矢量化，因此峰号 1 具有与之相关的“alfa [1]”、“峰 [1]”和“高度 [1]”。

到目前为止，我得到了预期的公式：

尽管如此，当我尝试为par生产线复制相同的系统时，我遇到了一些问题。这应该显示：

但是，当我折叠所有字符串并as.formula随后使用该命令时，我得到了：

如果我打印折叠的字符串，则字符行是预期的，所以我认为它将以某种方式链接到 as.formula 命令（即它可能不是适当的命令）

我正在尝试使用 ggplot2 绘制一些数据并拟合非线性曲线。我想将 stat_function 与我已经定义的 nls 对象一起使用，但结果会产生一条截断的曲线。我查看了以下页面，但到目前为止我还没有找到解决方案：

R - ggplot2 在线性区域外推回归线

使用 stat_function() 在 R 中的 ggplot 中绘制大量自定义函数

http://docs.ggplot2.org/0.9.3/stat_function.html

等效于 ggplot 的 curve()

我对ggplot2不是很有经验，所以如果我遗漏了一些简单的东西，我深表歉意，我很感激任何帮助。这是我的示例数据：

到目前为止，这是我用 ggplot2 制作情节的内容：

所有这些都很好，但是当我制作一个 nls 对象并尝试将它与 stat_function 一起使用时，我得到一条被截断的曲线。我想将曲线向下延伸到 x 轴。还会产生错误。

我想要的是更像以下使用基本 r 功能的东西，但使用 ggplot2 提供的所有花里胡哨，以及每种处理（Tx）的最终一条 nls 曲线。

这是迄今为止我在 R 中所做的最具挑战性的事情，因为 nls 和 LPPL 对我来说都是相当新的。

下面是我一直在使用的脚本的一部分。df 是一个数据框，由两列组成，日期和 Y，它们是标准普尔 500 指数的收盘价。我不确定它是否相关，但日期从 01-01-2003 到 12-31-2007。

当它运行时，我收到以下消息：

LPPL 的公式可以在此 pdf 文件的第 5 屏幕上找到，http://www.chronostraders.com/wp-content/uploads/2013/08/Research_on_LPPL.pdf

你知道我哪里错了吗？这适用于不同的模型，我更改了新方程的代码。这篇文章中的这段代码归功于 jlhoward，在 R 中使用 nls 重新创建研究。

谢谢您的帮助。

根据 jlhoward 的评论，可以在此处下载 df.rda：https ://drive.google.com/file/d/0B4xAKSwsHiEBb2lvQWR6T3NzUjA/edit?usp=sharing

我想用来y=a^(b^x)拟合下面的数据，

当我使用非线性最小二乘法时，

它会产生一个错误：初始参数估计处的奇异梯度矩阵。结果大致是 a = 1.1466，b = 0.6415，所以初始参数估计应该没有问题，因为我将它们定义为 a=1，b=0.5。

我在其他主题中阅读过修改曲线很方便。我在考虑类似的东西log y=log a *(b^x)，但我不知道如何处理函数规范。任何想法？

我想通过 r 平方值比较三个模型的曲线拟合。我使用nlsanddrc包运行模型。但是，这些软件包似乎都没有计算 r 平方值。他们给出了“残余标准误差”和“残余平方和”。

这两个可以用来比较模型拟合吗？

我正在尝试使用 nls() 函数来拟合信息。公式为：

但是R会生气:)

以下代码为(for )nls引发以下错误。但是类似数据集的相同代码可以正常工作（for ）。Convergence failure: singular convergence (7)fm2Data2fm1Data1

适用于该数据集

不适用于此数据集

我有我想使用 R 拟合以下等式的数据：

Z(u,w)=z0*F(w)*[1-exp((-b*u)/F(w))]

其中 z0 和 b 是常数，F(w), w=0,...,9 是一个递减阶跃函数，它取决于 w 且 F(0)=1 和 u=1,...,50。

Z(u,w) 是以 50x10 矩阵的形式观察到的一组数据（u=50,...,1 在行的一侧，w=0,...,9 沿列）。例如，由于我没有很好地解释，Z(42,3) 将是第 9 行向下和第 4 列中的元素。

使用 F(0)=1 我能够仅使用第一列（即 w=0）和代码来获得 b 和 z0 的估计值：

然后，我通过遍历每一列并使用我找到的 b 和 z0 的值，找到了 w=1,...,9 的 F(w)。

但是，我想找到一种方法来一次估计所有 12 个参数（b、z0 和 F(w) 的 10 个值），因为 b 和 z0 应该适合所有数据，而不仅仅是第一列。

有谁知道这样做的任何方法？所有帮助将不胜感激！

谢谢詹姆斯

R中有没有办法找出图中多条线的“一般”方程？下图是我想要实现的示例。我只显示了 t 的主线，但在这些主线之间也有隐藏线。我想要的是得到线条的“一般”方程。我该怎么做呢？

有人知道如何解决这个 R 问题吗？它是在关系中找到一个变化点，比如下面数据中的 x=5。

SAS 运行良好；

R失败；

由于这是一个信号清晰的简单案例，因此 R 应该可以工作。我想知道我在 R 中做错了什么（我使用了来自https://stats.stackexchange.com/questions/7527/change-point-analysis-using-rs-nls的一些代码）。有人建议/解决方案吗？

谢谢！

威廉