关于折线回归

一、说明

        今天的帖子主要是关于使用折线回归找到最佳值。即将某条曲线分解成包络线段,然后用分段回归方式优化。但它也涉及使用 SAS 和 R 的剂量反应研究和样条曲线。这不是第一篇关于这些主题的文章,但我确实想在其中添加折线。只是因为它还在使用。

二、药物剂量的应用

        寻找最佳值并不是什么新鲜事。为了找到最佳值,您通常会进行剂量反应研究。剂量反应研究旨在找出物质的数量或接触量之间的关系(剂量)及其对动物的总体影响(反应)。要有效地设计剂量反应研究,您需要了解生物系统并提出一个有针对性的问题:

  1. 有剂量效应的生物学证据吗?
  2. 最佳剂量是多少?
  3. 需要什么剂量才能产生不同于对照/零反应的反应?
  4. 剂量关系的本质是什么?

        不同的变量可以通过不同的方式建立不同的关系。下面是示例(最大值、最小值、山脊和鞍点)。这些内容之前已在此处进行了展示和解释。

剂量反应关系(顺时针):最大、最小、鞍形和脊形。

        剂量反应最好以 3D 而不是 2D 方式探索,并且最好使用不同的视角来显示。正如您所看到的,两个或三个变量是可行的,但想象一下必须以图形方式表示 6 因素交互作用。如果没有某种形式的聚类,几乎不可能做到这一点。

响应面示例,在 SAS 中进行分析。

        每种剂量反应设计成功的关键是了解生物学。这听起来可能是一个悖论——如果您已经对剂量反应曲线有了很好的了解,为什么还要进行剂量反应设计呢?这是因为许多物质都遵循给定的化学式。因此,如果您已经了解了曲线,则可以最大化研究设计。这也带来了一个问题,因为方程的微小变化可能会导致剂量反应曲线的大幅波动。

        从 2D 到 3D。从剂量反应(左)到反应面(右)。响应面是一个山脊,这意味着我们无法看到整个关系,但它暗示着最小值。

有大量关于如何进行剂量反应研究的材料。

总之,如果您想找到最佳值,您需要以这样的方式设计和分析研究,以便您可以找到最佳值(如果它确实存在)。这意味着:

  1. 了解所包含变量的生物学。
  2. 了解他们是否以及如何相互作用。
  3. 创造一个可以实现最佳效果的机会之窗。
  4. 为一项研究提供动力,使机会之窗(如果它确实存在)变得可见。

三、样条曲线:别担心,我会保持简短

        这不会是我关于样条线的第一段文字,也不会是我的最后一段,但我将主要使用图形来保持简短。下面,您可以看到样条线基本上是如何由连接件组成的。该图还将帮助您理解什么是折线回归。这是一条分段线性线。

        样条曲线已经超越了许多形式的旧回归方法,例如折线或多项式回归,可以找到最佳值。这就是为什么我们在解决折线回归之前显示样条线的原因。虽然样条建模有更多的建模自由度,但它也有更多的自由度来创建一些非常奇怪的关系。

        以上三种回归形式:线性、样条、LOESS。执行样条或 LOESS 回归很容易,但过度拟合的危险是真实存在的。

        这里又举了两个例子,其中制作没有实际意义的曲线的危险是很容易完成的。请记住:我从一开始就说过,生物学、物理学、化学或任何潜在的因果机制都是至关重要的。并且优于无意识的统计拟合。

四、折线回归:以及为什么我们应该将其称为分段回归

折线回归有多个名称,其中包括:

1.分段回归

2.分段回归

3.断棍回归

4.曲棍球棒回归

也许,与其称之为折线,分段回归是一个更能说明问题的名称,因为它是该方法的核心。分段回归模型具有两个或多个子模型,每个子模型在解释变量的单独域上定义。最简单的分段回归模型假设:

  1. 当 X 小于某一阈值时,响应由一个参数模型建模。
  2. 当 X 大于阈值时,响应由一个参数模型建模。

阈值也称为断点、分割点、连接点或结点。

分段建模意味着分解数据。

        但是,正如我之前所说,还有更高级的数据建模方法。尽管如此,这不应消除使用折线回归或分段回归(或分段回归——无论你怎么称呼它)的需要。

        如今,建模变得相当容易。通过使用自动化,我们可以让算法找到最佳的切点。这种情况发生在各种回归中:多项式、分段、样条、LOESS、GAM 等。每次都会评估结果的可能性,直到找到最小值。

        但是,我们也可以手动完成。它所需要的只是了解底层数据,并寻找断点有意义的点。最后,断点是回归线最好采取不同路线以维持最小 RMSE 的点。

        您可以在下面看到,无论是在理论上还是在 SAS 中,如何使用不同的函数来拟合非线性曲线。折线是有道理的,但很多时候并不是最佳解决方案。

        正如您所看到的,折线模型包含这样的假设:末端部件遵循一条不同的路径到达末端。这就是为什么单个断点很少是不够的,并且需要添加更多断点。添加足够的量,最终会得到一条样条线。

综上所述,折线回归应该称为分段回归。分段回归只不过是:

  1. 分解数据
  2. 分别分析一下
  3. 连接数据
  4. 迭代该过程

记住!折线通常不能提供模型的最佳拟合。

五、SAS 和 SAS 中的折线回归R:为什么软件并不重要

数据并不总是适合任何形式的分段建模。明智地选择了模特。

尝试在此类数据上拟合折线是没有意义的。

        然而,这并不意味着不可能,正如您从下面的结果中看到的那样。不过,说实话,我宁愿符合大均值。

在SAS

        我们也可以在 R 中使用用于分析分段模型的包来执行相同的操作。

        在 R 中,您可以看到断点 (9.8) 和系数。

SAS 和 R 之间的结果比较。

正如您所看到的,数据来自随机块设计,因此我们可以添加随机组件 - 每个块 (ID) 一条折线。

每个 ID 的分段回归

        最后,我上面所做的大部分工作都是由 SAS 或 R 等程序自动完成的。请记住,简单的部分是建模,困难的部分是首先选择为什么要建模。如果一条断线真的有意义的话。

@marc.jacobs012

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/233504.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

使用docker-compose搭建docker私服与配置WebUI

简介 本文介绍了使用docker compose 搭建 docker私服 环境 Docker version 24.0.6, build ed223bc Docker Compose version v2.21.0 正文 一、创建registry文件夹 我的路径是/usr/loca/docker/registry 二、创建并编写docker-compose.yml version: "3.9" services…

ios备忘录怎么导入华为 方法介绍

作为一个常常需要在不同设备间切换的人,我深知备忘录的重要性。那些突如其来的灵感、重要的会议提醒、甚至是生活中的琐碎小事,我们都习惯性地记录在备忘录里。但当我决定从iPhone转向华为时,一个问题困扰了我:如何将那些珍贵的备…

智能优化算法应用:基于学校优化算法3D无线传感器网络(WSN)覆盖优化 - 附代码

智能优化算法应用:基于学校优化算法3D无线传感器网络(WSN)覆盖优化 - 附代码 文章目录 智能优化算法应用:基于学校优化算法3D无线传感器网络(WSN)覆盖优化 - 附代码1.无线传感网络节点模型2.覆盖数学模型及分析3.学校优化算法4.实验参数设定5.算法结果6.…

Axure中继器完成表格的增删改查的自定义元件(三列表格与十列表格)

目录 一、中继器 1.1 定义 1.2 特点 1.3 适用场景 二、三列表格增删改查 2.1 实现思路 2.2 效果演示 三、十列表格增删改查 3.1 实现思路 3.2 效果演示 一、中继器 1.1 定义 在Axure中,"中继器"通常指的是界面设计中的一个元素,用…

IntelliJ IDE 插件开发 | (三)消息通知与事件监听

系列文章 IntelliJ IDE 插件开发 |(一)快速入门IntelliJ IDE 插件开发 |(二)UI 界面与数据持久化IntelliJ IDE 插件开发 |(三)消息通知与事件监听 前言 在前两篇文章中讲解了关于插件开发的基础知识&…

【linux】解决matplotlib中文显示乱码

一、Debian操作系统版本号 $ cat /etc/debian_version 12.4 二、问题 matplotlib中文显示乱码。 三、解决方法: plt.rcParams["font.sans-serif"] "Noto Sans CJK JP" 四、解决过程 1.问GPT 回答了一些相关问题,但是好像用…

Python tkinter 初探Toplevel控件搭建父子窗口

目录 Toplevel控件搭建父子窗口 最简明的父子窗口框架 改进一:屏蔽和开放按钮 改进二:子窗口始终在主窗口之上 改进三:增加子窗口的关闭协议 改进四:使子窗口长获焦点 总结 Toplevel控件搭建父子窗口 最近,用P…

k8s集群1.23.0版本部署说明

1.部署 k8s1.23.0版本与1.26.0版本的部署基本差不多,只不过k8s 1.23版本不需要部署cri-docker,所以只需要在1.26.0版本部署的基础上不要cri-docker的部署即可 参考:kubeadm部署k8s 1.26.0版本高可用集群_kubeadm 高可用集群-CSDN博客 搭建…

行为型设计模式(一)模版方法模式 迭代器模式

模板方法模式 Template 1、什么是模版方法模式 模版方法模式定义了一个算法的骨架,它将其中一些步骤的实现推迟到子类里面,使得子类可以在不改变算法结构的情况下重新定义算法中的某些步骤。 2、为什么使用模版方法模式 封装不变部分:模版…

搭建知识付费平台?明理信息科技为你提供全程解决方案

明理信息科技saas知识付费平台 在当今数字化时代,知识付费已经成为一种趋势,越来越多的人愿意为有价值的知识付费。然而,公共知识付费平台虽然内容丰富,但难以满足个人或企业个性化的需求和品牌打造。同时,开发和维护…

【容器Docker】Docker学习笔记

1、什么是Docker: Docker 将程序和程序运行所依赖的所有环境都打包到镜像里。“build once, run anywhere”Docker 是容器的一种实现。 Windows 下如何安装Docker: 官方安装教程:Install Docker Desktop on Windows | Docker Docs有两种安装套装&…

使用openMVS库,在VS2022中启用c++17标准编译仍然报错

使用openMVS库,在VS2022中启用c17标准编译仍然报错 现象 项目中引用了某些开源库(例如openmvs2.1.0),编译时要求启用编译器对c17的支持。 没问题!大家都知道在下图所示的位置调整C语言标准: 但是&#…

智能优化算法应用:基于类电磁机制算法3D无线传感器网络(WSN)覆盖优化 - 附代码

智能优化算法应用:基于类电磁机制算法3D无线传感器网络(WSN)覆盖优化 - 附代码 文章目录 智能优化算法应用:基于类电磁机制算法3D无线传感器网络(WSN)覆盖优化 - 附代码1.无线传感网络节点模型2.覆盖数学模型及分析3.类电磁机制算法4.实验参数设定5.算法…

ArcGIS Pro中去除黑边方法汇总

在有些时候(比如镶嵌栅格后),我们获取到的影像数据可能会有黑边,这里为大家汇总一下在ArcGIS Pro中去除黑边的方法,希望能对你有所帮助。 数据来源 本教程所使用的数据是从水经微图中下载的影像数据,除了…

vmware离线安装docker-compose

vmware离线安装docker-compose 最近安装docker-compose,发现git取拉取,不是拒绝连接就是报443错误,或者其他错误 最后发现用包直接传上去好用,不用git拉取了 离线安装docker-compose 本文章给的docker-compose离线包,…

JNDI注入Log4jFastJson白盒审计不回显处理

目录 0x00 前言 0x01 Maven 仓库及配置 0x02 JNDI 注入简介 0x03 Java-第三方组件-Log4J&JNDI 0x04 Java-第三方组件-FastJson&反射 0x05 白盒审计 - FastJson 0x06 白盒审计 - Log4j 0x07 不回显的处理方法 0x00 前言 希望和各位大佬一起学习,如果…

多线程面试题

文章目录 1. 如何停止正在运行的线程2. 请你谈谈JMM(java内存模型)3. AQS4. ReentrantLock实现原理5. 死锁怎么检测 1. 如何停止正在运行的线程 设置一个共享变量作为线程退出的标记,当这个标记不满足时while循环,线程一直运行&a…

uniapp图片上传说明

目录 1.文件上传组件 2.单文件上传 3.多文件上传 4.注意点 1.文件上传组件 前端上传组件使用uni-file-picker,可以自行进行下载使用。默认上传到绑定的服务空间,配置属性auto-upload为false关闭自动上传,可以限定上传的是图片还是文件&am…

Databend 开源周报第 124 期

Databend 是一款现代云数仓。专为弹性和高效设计,为您的大规模分析需求保驾护航。自由且开源。即刻体验云服务:https://app.databend.cn 。 Whats On In Databend 探索 Databend 本周新进展,遇到更贴近你心意的 Databend 。 新增对 Delta 和…

使用Python编写简单网络爬虫实例:爬取图片

🍎个人主页 🏆个人专栏:日常聊聊 ⛳️ 功不唐捐,玉汝于成 目录 ​编辑 简介 步骤 1. 安装依赖库 2. 创建目录 3. 发送HTTP请求并解析页面 4. 查找图片标签并下载图片 注意事项 结语 我的其他博客 简介 网络爬虫是一种…