在 R 中获取图的二次方程项

2023-12-26

我需要找到我在 R 中绘制的图的二次方程项。 当我在 Excel 中执行此操作时,该术语出现在图表上的文本框中,但我不确定如何将其移动到单元格以供后续使用(以应用于需要校准的值),或者实际上如何在 R 中请求它。如果它在 R 中是可召唤的,它是否可以保存为一个对象以用于将来的计算?

这看起来应该是 R 中的一个简单的请求,但我找不到任何类似的问题。非常感谢任何人可以就此提供的任何帮助。


所有的答案都提供了您似乎想做的事情的各个方面,但到目前为止还没有将它们全部整合在一起。让我们考虑一下 Tom Liptrot 的答案示例:

fit <- lm(speed ~ dist + I(dist^2), cars)

这为我们提供了一个拟合线性模型,变量为二次dist。我们使用以下方法提取模型系数coef()提取器功能:

> coef(fit)
 (Intercept)         dist    I(dist^2) 
 5.143960960  0.327454437 -0.001528367

所以你的拟合方程(由于打印而进行四舍五入):

\hat{速度} = 5.143960960 + (0.327454437 * dist) + (-0.001528367 * dist^2)

(其中 \hat{speed} 是响应速度的拟合值)。

如果您想将这个拟合方程应用于某些数据,那么我们可以编写自己的函数来完成它:

myfun <- function(newdist, model) {
    coefs <- coef(model)
    res <- coefs[1] + (coefs[2] * newdist) + (coefs[3] * newdist^2)
    return(res)
}

我们可以像这样应用这个函数:

> myfun(c(21,3,4,5,78,34,23,54), fit)
[1] 11.346494  6.112569  6.429325  6.743024 21.386822 14.510619 11.866907
[8] 18.369782

对于一些新的距离值(dist),这就是你想要从 Q 中做的事情。但是,在 R 中,我们通常不会做这样的事情,因为,为什么用户必须知道如何从所有不同类型形成拟合值或预测值可以在 R 中安装的模型吗?

在 R 中,我们使用标准方法和提取器函数。在这种情况下,如果您想将 Excel 显示的“方程”应用于所有数据以获得此回归的拟合值,在 R 中我们将使用fitted()功能:

> fitted(fit)
        1         2         3         4         5         6         7         8 
 5.792756  8.265669  6.429325 11.608229  9.991970  8.265669 10.542950 12.624600 
        9        10        11        12        13        14        15        16 
14.510619 10.268988 13.114445  9.428763 11.081703 12.122528 13.114445 12.624600 
       17        18        19        20        21        22        23        24 
14.510619 14.510619 16.972840 12.624600 14.951557 19.289106 21.558767 11.081703 
       25        26        27        28        29        30        31        32 
12.624600 18.369782 14.057455 15.796751 14.057455 15.796751 17.695765 16.201008 
       33        34        35        36        37        38        39        40 
18.688450 21.202650 21.865976 14.951557 16.972840 20.343693 14.057455 17.340416 
       41        42        43        44        45        46        47        48 
18.038887 18.688450 19.840853 20.098387 18.369782 20.576773 22.333670 22.378377 
       49        50 
22.430008 21.93513

如果您想将模型方程应用于一些未用于拟合模型的新数据值,那么我们需要从模型中获取预测。这是使用以下方法完成的predict()功能。使用我插入的距离myfun如上所述,这就是我们以更以 R 为中心的方式做到这一点的方式:

> newDists <- data.frame(dist = c(21,3,4,5,78,34,23,54))
> newDists
  dist
1   21
2    3
3    4
4    5
5   78
6   34
7   23
8   54
> predict(fit, newdata = newDists)
        1         2         3         4         5         6         7         8 
11.346494  6.112569  6.429325  6.743024 21.386822 14.510619 11.866907 18.369782

首先,我们创建一个新的数据框,其组件名为"dist",包含我们想要从模型中获得预测的新距离。值得注意的是,我们在此数据框中包含一个与创建拟合模型时使用的变量同名的变量。这个新数据框必须包含用于拟合模型的所有变量,但在这种情况下我们只有一个变量,dist。另请注意,我们不需要包含任何有关dist^2。 R 将为我们处理这个问题。

然后我们使用predict()函数,给它我们的拟合模型并提供刚刚创建的新数据框作为参数'newdata',为我们提供新的预测值,该值与我们之前手动预测的值相匹配。

我掩盖的一点是predict() and fitted()实际上是一整组函数。有以下版本lm()模型,用于glm()模型等。它们被称为generic函数,与methods(如果您愿意的话,可以提供版本)适用于几种不同类型的对象。你这个用户一般只需要记住使用fitted() or predict()等等,而 R 负责针对您提供的拟合模型类型使用正确的方法。以下是基本 R 中可用的一些方法fitted()通用函数:

> methods(fitted)
[1] fitted.default*       fitted.isoreg*        fitted.nls*          
[4] fitted.smooth.spline*

   Non-visible functions are asterisked

根据您加载的其他软件包,您可能会获得更多。这*只是意味着你不能直接引用这些函数,你必须使用fitted()R 计算出要使用哪一个。请注意,没有方法lm()对象。这种类型的对象不需要特殊的方法,因此default方法将被使用并且是合适的。

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

在 R 中获取图的二次方程项 的相关文章

随机推荐

  • SQL Server 多行计算

    如果我有以下格式的数据 Account Period Values Revenue 2013 01 01 5432 Revenue 2013 02 01 6471 Revenue 2013 03 01 7231 Costs 2013 01
  • Android ListFragment更新/刷新和onItemClick

    我的 Android 项目遇到问题 它是一个音乐曲目播放器 显示 ListFragment 中的每个曲目 每个项目都有一个位图 其中示波器绘制在位图画布上 Listfragment 中的项目会定期刷新 这是通过以下方式完成的 通知数据集更改
  • git浅克隆到特定标签

    我想克隆 Linux 内核存储库 但只能从 3 0 版本开始 因为内核存储库太大了 如果我可以进行浅克隆 它会使我的版本控制工具运行得更快 我的问题的核心是 我如何告诉 git depth 参数的 n 值是什么 我希望这会起作用 克隆htt
  • 如何使用Python从最高到最低设置k-Means聚类标签?

    我有一个包含 38 个公寓的数据集以及它们早上 下午和晚上的用电量 我尝试使用 scikit learn 的 k Means 实现对该数据集进行聚类 并得到了一些有趣的结果 First clustering results 这一切都很好 对
  • 显示在合并中丢失的观察结果

    假设我想通过两列的键合并两个不同的数据帧 Dataframe One 有 70000 个 obs 共 10 个变量 Dataframe Two 有 4500 个 obs 共 5 个变量 现在 我使用此代码检查了新数据框中的观察结果是如何留下
  • View.inflate 与 LayoutInflater

    使用上主要区别是什么LayoutInflater和静态方法View inflate 使用它们中的任何一个是否有任何缺点 或者它们可能有不同的用途 如果看源码的话View inflate 我们看到这个 public static View i
  • 在我的活动过程中保持屏幕处于唤醒状态

    我的应用程序中有三项活动 我想在第二个活动时保持屏幕唤醒 除非手动按下 锁定 键 否则屏幕在我的第二个活动中不应关闭 我浏览了很多链接 但它们对我来说似乎不清楚 正如 Android 教程中所讨论的保持屏幕开启 https develope
  • 加载活动选项卡页面源的 Chrome 扩展

    我希望你可以帮助我 我正在尝试在 Chrome 中创建一个扩展 它将活动选项卡的源加载到变量中 到目前为止我有 清单 json name My Extension manifest version 2 version 0 1 descrip
  • 如何调整文本字体大小以适合文本视图

    android 有没有办法调整文本视图中的文本大小以适应其占用的空间 例如 我正在使用一个TableLayout并添加几个TextViews 到每一行 因为我不想要TextView我宁愿看到它降低了内容的字体大小 有任何想法吗 我努力了me
  • 使用 gradle-experimental:0.2.0 签署 android apk

    我已经成功修改了我的 build gradle 以与 gradle experimental 一起使用 但是当尝试添加签名配置块时它失败了 我按照以下指示进行操作http tools android com tech docs new bu
  • Objective-C 标头中的冒号对齐

    我想使用 uncrustify 来对齐冒号上的参数 但取得了一些有限的成功 使用中的信息Uncrustify 可以在 Objective C 方法调用中对齐冒号吗 https stackoverflow com questions 9387
  • 更改 vuetify 网格中的列宽

    我在行下有三列 每列的 cols 值为 4 我想减小 col1 的宽度 我尝试将 cols 值从 4 更改为 3 但列变得比我预期的要小 所以我想要一个介于 3 到 4 之间的 cols 值 考虑 cols 值为 3 5 我有办法实现这个目
  • tycho 插件 + maven-dependency-plugin:从本地项目而不是存储库复制依赖项

    主要目标 将项目部署为 jarandeclipse 插件 当前状态 项目构建良好 作为 jar 包 现在我想创建第二个项目 将 jar 项目包装为 eclipse 插件 使用 tycho maven plugin 创建 eclipse pl
  • 将 generated-sources 作为源文件夹添加到 Eclipse

    我正在使用 maven jaxb plugin 基于 xsd 文件生成类文件源
  • PHP 重定向相同地址不同端口

    现在我有这个 header Refresh 0 url http 192 168 100 100 10500 redirect2 php 我怎样才能做同样的重定向但不写地址 只写端口 这两个文件位于同一服务器上的同一文件夹中 问题是我不知道
  • 获取 cURL 以跟踪重定向

    我在 Laravel 中有一个应用程序 可以抓取网站并检索特定信息 最近我看到越来越多的网站出现 403 禁止访问 所以我决定第一次尝试 cURL 我现在遇到的唯一问题是 对于具有 301 或 302 重定向的网站 cURL 不遵循它们 这
  • 允许 UIImagePickerController 编辑视频但不能编辑图像

    上传图像或视频Whatsapp 似乎使用UIImagePicker 可以在该视图中编辑视频 但无法编辑图像 貌似在SDK中 allowsEditing属性决定是否允许对图像和视频进行编辑 我怎样才能获得像Whatsapp一样的行为 其中视频
  • CSS表格,反转显示内容的顺序

    我知道这有点前沿 但无论如何 这就是问题 Given div First Div div div Second Div div one two display table cell 这给了我一个可爱的并排排列的div one 在左边 two
  • 如何从 ASP.NET 代码隐藏中获取当前 URL

    我的应用程序托管在不同的服务器上 我想获取当前服务器上页面的 URL 如何在代码后面获得这个属性 string url HttpContext Current Request Url AbsoluteUri http thehost com
  • 在 R 中获取图的二次方程项

    我需要找到我在 R 中绘制的图的二次方程项 当我在 Excel 中执行此操作时 该术语出现在图表上的文本框中 但我不确定如何将其移动到单元格以供后续使用 以应用于需要校准的值 或者实际上如何在 R 中请求它 如果它在 R 中是可召唤的 它是