问题标签 [nls]

我正在尝试使用带有自启动选项 ssLogis 的 nls 编写一个用于拟合逻辑曲线的函数。我有一些数据——实际上将作为一个数据框出现，所以这第一步实际上只是为了获得与我将要使用的数据类似的数据设置......

因此，在 Fox 和 Weisberg 的帮助下，我编写的函数...

我什至无法摆脱困境，看看我的函数代码的其余部分是否正常......运行函数时出现错误：

然后，我收到此错误：

phi2 是一个用 ssLogis 估计的参数，所以我迷路了。这就像函数外的冠军一样运行。

我已尝试搜索有关此主题的其他线程，但没有一个修复程序对我有用。我有一个自然实验的结果，我想显示一个事件的连续发生次数符合指数分布。我的 R shell 粘贴在下面

请原谅格式不好，先在这里发帖。x 包含直方图的 bin，y 包含该直方图中每个 bin 的出现次数。dat2 在 14 处截止，因为 0 计数箱会抛出指数回归，我真的只需要适合前 14 个。那些计数超过 14 的箱我有生物学理由相信它们是特殊的。我最初遇到的问题是无穷大，因为没有一个值是 0，所以我没有得到。在按照此处不同帖子的建议给出不错的起始值后，我得到了奇异梯度错误。我看到的唯一其他帖子有更多变量，我尝试增加迭代次数，但没有成功。任何帮助表示赞赏。一个

我想将以下函数拟合到我的数据中：

f(x) = 偏移+幅度sin(频率T+相位),

或根据维基百科： f(x) = C+alpha sin(omega T+phi)

我的数据存储在两列的文件中，我通过以下方式导入它们：

并通过以下方式转换它们：

因此，我得到了两个 double 类型的变量，一个称为“T”（相当于上面等式中的时间 T）和一个称为“值”（相当于上面等式中函数 f(x) 的值）。到目前为止，一切正常！但是，当我想将 f(x) 拟合到我的数据时，我会得到不好的结果（如果有的话）。这是我的代码：

我现在的问题是，在这个简单的四行示例中（至少据我所知）唯一可能出错的是起始值。但是，它们是 Gnuplot 拟合的结果（它使用最小二乘法）并且它们看起来非常好（因为我没有“10 声望”，所以我无法附上图，抱歉）。{edit1: plot attach} 上面代码的结果：

我的后续步骤是改变起始值以获得合理的结果（也许起始值太好了......） - 没有成功！所以我现在的问题是：问题出在哪里？我想使用 R，因为我希望我可以产生一个合适的结果，即 a) 更好的 b) 比 Gnuplot 更差的起始值（至少 b）对我来说是一个非常重要的标准，这就是我不想要的原因使用 Gnuplot）。

Ps.：是的，这是我使用 R 的第一天。因此，如果您对我有解决方案/理由，如果我能理解它会很好。

PPs.：对于那些想要获得相关信息的人：我已经收集了很多链接，所以如果我得到至少“10名声望”，我可以在这里发布。{edit2：附加链接}

[R] 拟合周期性“正弦波”模型

非线性回归和非线性最小二乘

R中的非线性回归和非线性最小二乘

Nichtlineare 回归（德语）

R：非线性最小二乘

如何将曲线拟合到正弦波

官方文档 -> 手册

链接罗兰在他的回答中提供

PPs.：这里是我使用的数据（由制表符分隔，仅使用第一列（“计数器”）和第三列（“最小值”）；由于最大字符数/帖子，我删除了最后一列的一些值{edit3 ：由于附加链接而删除了其他条目}）{edit4：从“第二列（“角度”）”更改为“第一列（“计数器”）”。在代码和答案中，它的第一列也是如此。}：

波斯语中的“KAF”（ک，ك）和“YA”（ی，ي）有不同的字母，具有不同的Unicode代码和相同的拼写。

我正在使用这些命令在我的项目中设置波斯排序和比较

通过设置nls_comp和nls_sort参数，我收到这些输出

这些查询如何获得相同的输出？我可以定制nls_comp吗？

我想使用Stata通过非线性最小二乘估计以下函数：

我正在测试另一篇论文的结果，并且想使用 Stata，因为它与我正在复制的论文中使用的软件/求解器相同，并且例如，它应该比使用 GAMS 更容易做到。

我的问题是我找不到任何方法来写出上面等式的总和部分。在我的数据中，我所拥有的只是一个单独的观察值，其中 j 的值位于单独的变量中。我可以用以下方式写出整个表达式（对于三个观察/我）：

有没有一种简单的方法可以告诉 Stata 对给定数字集的表达式/变量求和，例如在 GAMS 中您可以编写：

[ORACLE DATABASE] 设置后，alter session set nls_sort = GERMAN; 更改会话集 nls_language = 德语；

当以块加载数据并使用“>”进行比较以获得不同的排序块时，显示不正确的数据（给出重复的条目）。

前任。对于不同的块，在查询中使用它... name > 'Account' 或 name > 'Project' 或 name > 'BB' 将在每个输出中给出名称 'aa' 、 'bb' 、 'cc' 。

注意，使用 nls_sort=BINARY 和 NLS_language=ENGLISH 会给出正确的输出。

我想估计以下类型的非线性回归模型：

元素说明：

最后，i=1,2 的 b[i,theta] 表示以下指数阿尔蒙多项式加权函数：

它只代表 x1 和 x2 的两个衰减权重，仅此而已。但是这些权重取决于两个参数值，它们也需要估计：theta1 和 theta2。

现在，我想使用非线性最小二乘函数 nls() 估计参数 ß0、ß1、theta1 和 theta2 的最佳（相对于 RSS 标准）值。

我尝试了以下导致错误消息的方法：

注意：为了符号简单，我预先计算了加权函数的分母值，即 1318837781。

似乎nls()将 x2 和 x2 视为参数，但它们是回归量。我在这里做错了什么，我应该如何修改代码以获得合理的结果。或者不可能用 来估计这种函数nls()？

谢谢！

我正在使用用 R 编写的包，以使实验数据适合特定模型，并且我想编写自己的模型并使用相同的包（作者声称这是可能的）。

所以我正在挖掘源文件以找到它们定义模型函数的位置。我被他们使用的方式卡住了nls。我了解nls以下案例的工作原理

创建一个数据集并在它旁边编写模型。

但是它们使用nls如下所示：

当我跑步时

查看上面示例中的公式 (Const + A * exp(B * x))。我明白了

我想了解如何查看结果公式以及nls 数据设置为 list\environment 时如何工作。有什么建议么？

这里是内容rescomp

我创建了一个目标函数 f(x0)，它通过指定初始猜测 x0 返回实际值和理论值之间的平方和。如何使用 nls() 函数指定最小化问题？

提前致谢。

我正在尝试将 sigmoid 拟合到一些剂量反应数据。当我在 ggplot 之外使用 nls 时，它工作正常。当我在 geom_smooth 中使用它时，它会引发以下错误：

示例数据如下所示：

电话是：

和：

我做错了什么吗？