处理 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 错误 的相关文章

  • 查找嵌套列表中元素的索引?

    我有一个类似的列表 mylist lt list a 1 b list A 1 B 2 c list C 1 D 3 是否有一种 无循环 方法来识别元素的位置 例如如果我想用 5 替换 C 的值 并且在哪里找到元素 C 并不重要 我可以这样
  • 如何在R中得到一个大的稀疏矩阵? (> 2^31-1)

    我使用一些 C 代码从数据库中获取文本文件 并从该文件创建 dgcMatrix 类型稀疏矩阵Matrix包裹 我第一次尝试构建一个具有超过 2 31 1 个非稀疏成员的矩阵 这意味着稀疏矩阵对象中的索引向量也必须比该限制长 不幸的是 向量似
  • 简单的数据框重塑

    我刚刚从长时间的写作中断中回到 R 并且在记住如何重塑数据方面遇到了一些实际问题 我知道我想做的事情很容易 但出于某种原因 我今晚很愚蠢 并且将自己与融化和重塑混淆了 如果有人能快速指出我正确的方向 我将不胜感激 我有一个这样的数据框 pe
  • 将密度曲线拟合到 R 中的直方图

    R中有没有可以将曲线拟合到直方图的函数 假设您有以下直方图 hist c rep 65 times 5 rep 25 times 5 rep 35 times 10 rep 45 times 4 看上去很正常 但其实是歪曲的 我想拟合一条倾
  • 将 JSON URL 转换为 R 数据帧

    我在将 JSON 文件 从 API 转换为 R 中的数据帧时遇到问题 例如 URL 我尝试了 S O 的一些不同建议 包括将json数据转换为R中的数据框 https stackoverflow com questions 28683769
  • 使用 R 的 flextable 包时,有没有办法将传递给 add_header_lines() 的字符串部分加粗

    我正在使用我喜欢的 flextable 包为 Word 文档创建几个表格 但是 我在将表格标题中的部分文本加粗时遇到了一些麻烦 例如 我希望标题为 Table 1 我的表格标题的其余部分 而不是 表 1 我的表格标题的其余部分 I 找到这个
  • 在函数内部调用 clusterApply 时,性能会下降

    我遇到了一个奇怪的问题clusterApply 我已经能够尽可能地隔离它 如下所示 首先 我从全局环境运行以下代码 require parallel cl lt makeCluster rep localhost 20 SOCK xl lt
  • 限制数据框中所有单元格的字符串长度?

    您好 有没有一种方法可以限制 data frame 中所有列的字符串文本大小 而不必循环遍历每一列并一次使用 str trunc 之类的东西 例如下面的数据框 我可以将所有文本大小限制为仅 5 个字符 而不必一次只执行一列吗 如果有 50
  • 如何获得属于五分位数的x?

    我正在大学学习使用 R 进行计量经济学项目 所以请原谅我的笨拙 基本上 使用并给出 一个矩阵 股票价格 行 天 列 公司股票价格 另一个矩阵 市值 行 天 列 公司市值 我必须收集第三个矩阵每天观察的属于市值分布第一五分位数的股票价格 然后
  • 错误:“rjags”的包或命名空间加载失败

    在终端的 conda 环境之一中 我能够成功安装包 rjags 但是 当我在该环境中运行 R 并运行库 rjags 时 出现以下错误 加载所需的包 coda 错误 rjags 的包或命名空间加载失败 rjags 的 loadNamespac
  • 当每个记录都是一个段落并且某些记录有 4 个字段而其他记录有 6 个字段时,如何将文本文件读入 R

    如何读取文本文件 其中每条记录都是一个段落 每个换行符表示单独的字段 复杂的是 有些记录有 4 行 有些记录有 6 行 当字段数量的差异为 1 时 DWin 解决了我的问题 但当字段数量差异为 2 时 一切都崩溃了 你可以有一个在这里看看他
  • 在另一个 Rmd 中运行选定的块

    我已经在源 Rmd 文件中运行了分析 并且希望仅使用few来自源的块 我已经看到了一些关于从源 Rmd 中提取所有块的答案来自另一个 Rmd 中的 Rmd 文件的源代码 https stackoverflow com questions 4
  • 将 sf voronoi 多边形裁剪到边界框时出错

    我正在尝试将 voronoi polygons 使用 sf package 创建 剪辑 到边界框 但它引发了我无法定义的错误 我对 R 的空间世界不太有经验 感谢所有帮助 样本数据 stations lt structure list ST
  • 如何更改 Shiny 中 navbarPage 折叠的断点

    我想用shiny navbarPage collapsible TRUE 当在小屏幕上查看我的 Shiny 应用程序时 将导航元素折叠到菜单中 默认情况下 当浏览器宽度小于 940 像素时会触发折叠 有什么方法可以改变这一点 以便在稍大的浏
  • fread 将空导入为 NA

    我正在尝试导入带有空白的 csv 读取为 不幸的是他们都读作 NA now 为了更好地演示问题 我还展示了如何NA NA and 都映射到同一事物 除了最底部的示例 这将妨碍简单的解决方法dt is na dt lt gt write cs
  • glmnet 未从 cv.glmnet 收敛 lambda.min

    我跑了20倍cv glmnet套索模型以获得 lambda 的 最佳 值 但是 当我尝试重现结果时glmnet 我收到一个错误 内容如下 Warning messages 1 from glmnet Fortran code error c
  • 如何计算嵌套函数中的粘合表达式?

    我正在尝试嵌套一个函数 该函数将两个字符串粘合在一起 该函数使用组合字符串来命名数据帧的列 然而 问题似乎是粘合表达式没有足够早地评估为字符串 我可以 并且应该 强制在将表达式作为参数传递给另一个函数之前对其进行求值吗 library ti
  • ggplot2、R 中的单条形条形图

    我有以下数据和代码 gt ddf var1 var2 1 aa 73 2 bb 18 3 cc 9 gt gt dput ddf structure list var1 c aa bb cc var2 c 73L 18L 9L Names
  • 具有 dplyr、tidyverse 和 broom 的相关矩阵 - P 值矩阵

    全部 我想使用以下方法从相关矩阵中获取 p 值dplyr 和 或扫帚包 并同时测试多个变量 我知道其他方法 但 dplyr 对我来说似乎更简单 更直观 此外 dplyr 需要关联每个变量以获得特定的 p 值 这使得该过程更容易 更快 我检查
  • 非闪亮上下文中的反应式对象绑定

    实际问题 你怎样才能近似反应性环境 行为 http shiny rstudio com tutorial lesson6 建立者shiny http shiny rstudio com函数 或者甚至可能在一个函数中使用这些函数无光泽上下文以

随机推荐

  • 颤动曲线条

    我想知道是否有更好的解决方案来制作如下图所示的弯曲条 这是我的颤振代码 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 错误情况 完美契合 这是我