处理 nls - R 脚本中的 0 错误

2023-12-23

有什么方法可以让我的 nls 在进行非线性拟合时具有 0 残差误差吗?我的数据中有一些情况,其中拟合应该有 0 错误,但 nls 总是失败并吐出错误。

谁能告诉我:

  1. 如何测试这是否是 nls 吐出的错误?
  2. 如何允许 0 错误情况? (完美契合)

这是我的 nls 电话:

fit <- nls(y ~ ifelse(g, m1 * (x - x0) + y0, m2 * (x - x0) + y0),
            start = c(m1 = -1, m2 = 1, y0 = 0, x0 = split),
            algorithm = "port",
            lower = c(m1 = -Inf, m2 = -Inf, y0 = -Inf, x0 = split),
            upper = c(m1 = Inf, m2 = Inf, y0 = Inf, x0 = (split+1)),
            data=data.frame(x,y))

As 提及 https://stackoverflow.com/questions/7127940/r-script-nls-not-working/7128102#7128102在之前的回答中,?nls 明确指出你不应该使用nls对于 0 错误数据。要直接引用您正在使用的函数的帮助文件:

不要对人工“零残差”数据使用 nls。

nls 函数使用相对偏移收敛准则 比较当前参数估计的数值不精确性 到残差平方和。这对表单数据表现良好

y = f(x, θ) + eps

(var(eps) > 0)。它未能表明数据的收敛性 形式

y = f(x, θ)

因为该标准相当于比较两个组成部分 舍入误差。如果您想在人工数据上测试 nls,请添加 噪声分量,如下例所示。

一个潜在危险的选择是使用warnOnly = TRUE强迫nls在收敛之前返回,仅发出警告(无错误):

x <- -(1:100)/10
y <- 100 + 10 * exp(x / 2)
nlmod <- nls(y ~  Const + A * exp(B * x),control = nls.control(warnOnly = TRUE))

上面的例子也几乎直接取自?nls.

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

处理 nls - R 脚本中的 0 错误 的相关文章

  • R 中的插补 MICE 仍不存在于数据集中

    运行 MICE 包后 5 个完整插补集中的缺失值数量从 147428 减少到 46093 但不是应该是 0 NA 吗 Thanks 这是我的 MICR 代码 imp mice newdata imputationSet1 complete
  • 将嵌入其他文本的长州名称转换为两个字母的州缩写

    我的目标是识别以具有其他文本的字符向量写出的美国各州 并将各州转换为缩写形式 例如 北卡罗来纳州 到 北卡罗来纳州 如果向量只有长形式的状态名称 那就很简单 但是 我的向量在随机位置有其他文本 如示例 states 中所示 states l
  • r dplyr - 读取文件列表并使用文件名作为变量

    我想用每个文件的名称替换 image annotate 下面的 文本 library dplyr library purrr library magick list files path pattern png full names T g
  • 当我创建新变量时出了什么问题?

    我想根据原始变量施加的条件创建一个新变量 比方说 原始变量 var 是由 1 20 中的随机样本组成的向量 并且 当原来的 var 大于10时 新变量 newvar 被设置为缺失 当 var 小于10时 新变量 newvar 被设置为等于
  • 使用 ``magrittr::`%>%` `` 时 magrittr 管道出错

    不管出于什么原因我在玩magrittr管道语法 并遇到一个奇怪的错误 当您 scope 显式限定调用时发生 gt 我知道使用下面的语法会破坏管道的用途 但我很好奇为什么会发生错误 第一次致电sum按预期工作并输出1 第二次调用会导致错误 E
  • 按钮:带滚动条的下载按钮仅下载几行

    我正在处理超过 100 000 行的表并使用DT包裹 开发版本0 1 56 在 Shiny App 中将其可视化 此外我正在使用DT扩展如下 Buttons 下载不同格式的数据 然而虽然Scroller扩展程序也已激活 我只能下载几行 不是
  • 无法安装bigrf包

    我在安装 bigrf R 软件包时遇到了一些极其令人沮丧的困难 我已经尝试过以下方法 install packages bigrf repo https github com aloysius lim bigrf git type sour
  • 使用 google 查询邮政编码距离

    我有两个邮政编码列表 R 语言 其中一个是孩子的地址及其学业成绩 另一个是学校的地址 我希望能够为每个孩子找到最近的学校 所以大概需要通过转换为长和纬度值来计算邮政编码之间的距离 然后我希望能够在谷歌地图上绘制每所学校的所有孩子 并看看住在
  • 根据 B 列中的一系列值获取 A 列的平均值

    我的数据框有几列 如下所示 df1 lt data frame A c 1 2 4 B c 1 3 1 C c 1 1 3 我有两个条件来获取 A 列的平均值 条件1 我想在B为1时获得A列的平均值 即只对row1和row2进行平均 条件2
  • 加载 plyr 包时出现问题

    我使用 R 2 13 1 但未能成功尝试在 R 中加载包 plyr 1 6 我已将其手动安装到目录 R library 中 我的代码是 libPaths R library library plyr 我收到消息 库 plyr 中的错误 pl
  • 将一个数据框中的行(带有行名称)与另一个数据框中的匹配列名称相乘

    我有两个数据框 df1 lt data frame Values c 0 01 0 05 row names c X Y df1 Values X 0 01 Y 0 05 df2 lt data frame c 0 1 1 c 1 0 0
  • 是否可以旋转 R 中的绘图(基本图形)?

    我搜索了这个 发现使用 grid 有多种方法可以旋转图像 并且对于某些绘图 您可以使用它们的旋转 例如plot x y 而不是plot y x 不过我想知道是否有R 中旋转绘图的通用方法 适用于基础图形中生成的任何绘图 您可以导出图形 将其
  • 如何使用键盘命令中断 R 中正在运行的代码?

    如何中断 R 中正在运行的脚本 假设我遇到无限循环或其他问题 并且我想停止脚本运行 有键盘命令可以做到这一点吗 如果有帮助的话 在 Python 中是 Ctrl c 我在R中尝试过 但没有成功 自我回答 几乎是其他人的评论和答案的总结 In
  • 通过另一个函数将数据和列名称传递给 ggplot

    我将直接跳到一个示例并在后进行评论 cont lt data frame value c 1 20 variable c 1 20 1 20 1 5 1 20 2 group rep c 1 2 3 each 20 value variab
  • 将 RMarkdown 文档编织为 Word 时方程式和引用丢失

    我不确定这个问题是否更适合LaTeX论坛 我将其发布在这里是因为我怀疑问题更多是关于knitr和 RMarkdown 相比于 LaTeX 我在 RStudio 中有以下 RMarkdown 文档 title Capricious Behav
  • R中的预测和预测函数之间的区别

    两者之间有什么区别吗predict and forecast R 中的函数 如果是 在哪些具体情况下应该使用它们 Intro predict 适用于多种 R 对象 模型 基础库的一部分 forecast 对于时间序列 预测包的一部分 参见示
  • 获取非零数据的列意味着

    R 可以获得数据帧非零值的 colMeans 吗 data lt data frame col1 c 1 0 1 0 3 3 col2 c 5 0 5 0 7 7 colMeans data 1 33 4 我想要这样的东西 mean dat
  • 使用两列数据分割数据帧并对结果数据帧列表应用通用转换

    我想根据两列中的值将一个大数据帧拆分为一个数据帧列表 然后 我想对结果列表中的所有数据帧应用通用数据转换 滞后转换 我知道 split 命令 但只能让它一次处理一列数据 您需要将所有想要分割的因素放在一个列表中 例如 split mtcar
  • 替换因子列中的

    我想更换
  • 将公式传递给 R 中的函数?

    对此的任何帮助将不胜感激 我正在使用 Lumley 调查包 并试图简化我的代码 但遇到了一些小障碍 在我的代码中调用包中的 svymean 函数如下 其中第一个参数是指示我想要哪些变量的公式 第二个参数是该数据集 svymean hq eh

随机推荐

  • 颤动曲线条

    我想知道是否有更好的解决方案来制作如下图所示的弯曲条 这是我的颤振代码 import package flutter web material dart class CurvedBar extends StatelessWidget con
  • 如何用我自己的 4x4 矩阵更新 opengl modelview 矩阵?

    我有 4x4 矩阵用于对象的变换 float mat44 16 但我不知道如何使用我的矩阵更新 OpenGL ModelView 矩阵 我应该使用 glTranslatef glRotatef 与矩阵中的相关值还是应该使用 glLoadMa
  • 初学者无法创建第一个 Android 应用程序项目

    我已将 ADT 捆绑包 内部版本 v22 0 1 685705 下载到我的 Macbook Pro OS X 10 8 4 上 并根据 Android 开发人员说明进行安装 6 月 19 日星期三 当我尝试使用所有默认设置创建新的应用程序项
  • 以编程方式创建并显示 UIPickerView

    我试图以编程方式创建 UIPickerView 并将其显示为文本字段的第一个响应者 但是 选择器视图没有显示 textField 连接到界面生成器中的对象 但 pickerView 是以编程方式创建的 class View UIViewCo
  • 单个整数的 Web API 响应类型

    调用仅返回单个整数的 Web API REST 服务的 正确 方法是什么 我对 XML JSON 或其他任何内容没有任何要求 对服务的调用只需要返回一个整数 我是否使用ResponseType属性在这里 我的服务返回类型为HttpRespo
  • 为什么 JTable 列中没有显示动画 .gif 图标?

    Here is the processing gif Here is initial png Here is the output 这是代码 processing gif 正在其他位置工作 例如在某个选项卡中JTabbedPane 这里在a
  • CSS3 Flexbox 项目之间的间距

    对于 Flexbox 来说有些陌生 尽管在 CSS 方面经验丰富 在我看来 我读过的大多数教程都方便地 掩盖 的一件事是 Flex 项目之间的间距 例如 被引用最多的教程之一是CSS Tricks 中的这个 https css tricks
  • 创建Web服务时Eclipse Axis错误

    每次我尝试在 Eclipse 中创建新的 Web 服务时 都会收到以下错误 IWAB0489E Error when deploying Web service to Axis runtime axis admin failed with
  • Mailkit 对 Imap 进行身份验证失败

    我正在尝试使用 MimeKit 和 MailKit 从收件箱中获取电子邮件 但登录邮件服务器Imap不起作用 我使用 993 端口和 ssl 我尝试使用 Telnet 在端口 143 上连接 登录 获取 它有效 我尝试使用 Openssl
  • 使用Reducer React方法

    我正在玩useReducerReact 的 api 并且想知道理论 文档 部分和我实现的部分之间的区别 带有 useReducer 钩子的组件的初始状态 const fields dispatch React useReducer form
  • 擦除突出显示的 JavaScript getSelection 函数

    突出显示不是问题 我的问题是 如何在 javascript 中使用 window getSelection 删除突出显示 并创建节点新跨度关闭和新跨度开口以擦除选定的突出显示区域 请参阅屏幕截图 function removeHighlig
  • 在Python中交错多个相同长度的列表[重复]

    这个问题在这里已经有答案了 在Python中 有没有一种好方法来交错两个相同长度的列表 说我被给予 1 2 3 and 10 20 30 我想把它们变成 1 10 2 20 3 30 发布问题后 我意识到我可以简单地执行以下操作 val f
  • 通过 C#.NET 将 ZPL 打印到 Zebra 打印机

    我正在建立一个网站 将 ZPL 代码打印到 Zebra 打印机 140XiIII 具体如下 打印机通过并行端口连接到计算机 打印机未连接到网络 网站管理员希望尽可能避免使用 Javascript 而是专注于 C 我尝试过使用代码将 ZPL
  • R 中的虚拟列,提供的示例[重复]

    这个问题在这里已经有答案了 我有一个看起来像这样的数据框 uid language name 1 333 English 2 333 French 3 333 Dutch 4 654 Spanish 5 546 English 6 546
  • 如何正确导入stddraw?

    编辑 IntelliJ CE 我想要的 能够写作 setCanvas 500 500 代替 StdDraw setcanvas 500 500 问题 我不知道如何正确导入 Stddraw 库 如果我只是这样做 import StdDraw
  • Yii2 Email 如何设置发件人姓名

    我使用 Mailer 发送电子邮件 所以我对发件人姓名有疑问 这是我的配置 mailer gt class gt yii swiftmailer Mailer useFileTransport gt false messageConfig
  • 如何从标准 ServiceStack 错误响应中删除堆栈跟踪

    我刚刚开始使用 ServiceStack 我想了解是否可以从标准错误响应中删除堆栈跟踪 我尝试关闭调试模式但没有任何运气 public override void Configure Funq Container container Set
  • Pandas 图,vars() 参数必须具有 __dict__ 属性?

    它之前工作得很好 但由于某种原因 现在我遇到了奇怪的错误 熊猫版本 1 2 3 matplotlib 版本 3 7 0 示例数据框 df cap Date 0 1 2022 01 04 1 2 2022 01 06 2 3 2022 01
  • 在链接点击的另一个网站上使用相同的 php 会话

    我有个问题 假设我在网站上域名1 com我点击一个链接域名1 com这让我想到域名2 com 但我想用同样的session id 设置为session start in 域名1 com on 域名2 com 我如何在 PHP 中实现这一点
  • 处理 nls - R 脚本中的 0 错误

    有什么方法可以让我的 nls 在进行非线性拟合时具有 0 残差误差吗 我的数据中有一些情况 其中拟合应该有 0 错误 但 nls 总是失败并吐出错误 谁能告诉我 如何测试这是否是 nls 吐出的错误 如何允许 0 错误情况 完美契合 这是我