如何绘制双对数 R 图的线性回归?

2023-12-11

我有以下数据:

someFactor = 500
x = c(1:250)
y = x^-.25 * someFactor

我以双对数图显示:

plot(x, y, log="xy")

现在我使用线性模型“找出”数据的斜率:

model = lm(log(y) ~ log(x))
model

这使:

Call:
lm(formula = log(y) ~ log(x))

Coefficients:
(Intercept)       log(x)  
      6.215       -0.250  

现在我想将线性回归绘制为红线,但 abline 不起作用:

abline(model, col="red")

将回归线添加到图中的最简单方法是什么?


lines(log(x), exp(predict(model, newdata=list(x=log(x)))) ,col="red")

在对数刻度上绘制的 x 值范围和用作自变量的 log(x) 值范围实际上有很大不同。这将为您提供全方位:

lines(x, exp(predict(model, newdata=list(x=x))) ,col="red")

enter image description here

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

如何绘制双对数 R 图的线性回归? 的相关文章

  • 如果在循环中调用summary()命令,如何获得它的输出?

    Suppose Z是特征名称的向量 如何让以下 Rscript 中的摘要命令实际打印 for var in Z cat i form paste crim var lm fit lm form data Boston summary lm
  • 聚合时间

    我的数据包含 1 分钟测量值 我希望将其聚合为 5 分钟和 10 分钟 如何在 10 分钟内将 日 列的数据分组为 5 样本数据2 structure list Day c Monday Monday Monday Monday Monda
  • 获取非零数据的列意味着

    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
  • R:从本地 tar.gz 安装 R 软件包时找不到“make”

    R 包ConvCalendar不再位于 Cran 存储库中 请参阅here https cran r project org web packages ConvCalendar index html 然而 因为我在之前的项目中大量使用了这个
  • 难以理解 R 中双括号和单括号子集之间的区别[重复]

    这个问题在这里已经有答案了 我很难理解双括号子集和单括号子集之间的区别 我在开源编程方面相当陌生 我很难理解 R 中的 help 函数 因为考虑到我目前对 R 的理解 其中的一些信息对我来说太技术性了 我尝试过谷歌搜索差异 虽然它给了我一个
  • 如何绘制对数似然函数图

    我想绘制 pi 和 pi 之间的对数似然函数 对数似然函数 llh lt function teta x sum log 1 cos x teta 2 pi x c 3 91 4 85 2 28 4 06 3 70 4 04 5 46 3
  • 仅当 ggplot 中 y 轴的下限设置为 0 时才会出现图条[重复]

    这个问题在这里已经有答案了 我正在尝试创建一个条形图 当我将限制设置为 0 7 时 就会出现条形 但是 我希望下限为 1 而不是 0 当我将下限设置为 1 时 条形图不会出现 我收到以下错误消息 Removed 8 rows contain
  • magrittr 管道中的 WOE

    如何将下面的证据代码权重放入 magrittr 管道中 df gt 我尝试过的一切似乎都不起作用 df library Information library magrittr df a c aa bb cc aa aa aa bb cc
  • 将公式传递给 R 中的函数?

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

    标题基本上说明了一切 如果我这样做 makeActiveBinding x function runif 2 GlobalEnv x 1 0 7332872 0 4707796 x 1 0 5500310 0 5013099 那我有什么办法
  • Pandas 中每列的曲线拟合 + 外推值

    我有一个包含大约 300 列的数据集 每一列都与深度相关 Pandas DataFrame 的简化版本看起来像这样 import matplotlib pyplot as plt import numpy as np import pand
  • 拆分并保存在新的 data.frames 中

    我有一个大 data frame 144 列 我想将其分成每组 3 列 子文件或子 data frame 然后将子 data frame 保存在单独的文件中 换句话说 file1 将包含从 1 到 3 的列 file2 将包含从 6 到 9
  • R + ggplot2 - 无法分配大小为 128.0 Mb 的向量

    我有一个 4 5MB 9 223 136 行 的文件 其中包含以下信息 0 0 0 0147938 3 67598e 07 0 0226194 7 35196e 07 0 0283794 1 10279e 06 0 033576 1 470
  • 如何计算满足条件的行数

    假设我有以下数据框 Data1 X1 X2 1 15 1 2 3 1 3 7 0 4 11 1 5 1 0 6 9 0 7 18 0 8 6 1 9 3 1 我想知道如何找到观察的总数X1大于 9 并且X2等于1 我想我需要使用sum 但我
  • 如果我有一个包含计数的表格,如何获得直方图/摘要统计数据? [复制]

    这个问题在这里已经有答案了 我习惯了像这样的观察向量c 1 1 1 3 4 我可以在上面看到汇总统计数据和绘制直方图等 但现在我发现自己拥有相同的数据 但采用如下压缩格式 value numObs 1 3 3 1 4 1 当数据以data
  • 列值的切换功能

    我有一个缩写变量 名称错误地分散在整个列表中 请参阅下面的示例结构 ID lt c SPW SM DLS SJ joe schmoe CEJ teddy roos GVF MJC LH sally fields Full names sho
  • 使用 ggplot 将条形图的列与线图的点对齐

    当线图的点与条形图的条具有相同的 x 轴时 有什么方法可以使用 ggplot 将它们对齐 这是我尝试使用的示例数据 library ggplot2 library gridExtra data data frame x rep 1 27 e
  • 如何删除括号内的值的行?

    我正在使用以下数据框 Name Height Eric 64 Joe 67 Mike 66 Nick 72 Dave 69 Steve 73 我想删除 名称 列以左括号 开头的所有行 因此最终的数据框如下所示 Name Height Eri
  • 将 dplyr 中的starts_with与部分列名称的向量一起使用

    我想使用 dplyr 选择与字符串向量匹配的某些列 one lt seq 1 10 two lt rnorm 10 three lt runif 10 1 2 four lt 10 1 df lt data frame one two th
  • Sweave + RweaveHTML:cat 输出未出现在输出中

    我对 Sweave RweaveHTML 有疑问 我希望 cat 的输出最终出现在正在生成的 html 文件中 我有一个案例 它没有 我不明白为什么 test function bla bla cat Result is 然后在 Rnw 文

随机推荐

  • BufferedReader 用于大 ByteBuffer?

    有没有一种方法可以使用 BufferedReader 读取 ByteBuffer 而不必先将其转换为 String 我想将相当大的 ByteBuffer 作为文本行读取 出于性能原因 我想避免将其写入磁盘 在 ByteBuffer 上调用
  • 仅当行不存在时才在 PL/pgSQL 中运行 SQL 语句

    我想在 Postgres 9 6 的 PL pgSQL 函数中做这样的事情 INSERT INTO table1 id value VALUES 1 a ON CONFLICT DO NOTHING If the above stateme
  • Sqlalchemy 的问题并将 jsonb 数组插入到 postgresql

    所以我试图将 jsonb 值数组插入到我的数据库中 但我似乎无法正确格式化它 这是我的代码 updated old passwords append index 1 password hashed password user old pas
  • 如何为express.js 服务器设置 SSL 证书?

    之前 在旧版本的 Express 中 我可以这样做 express createServer key keyFile cert certFile 然而 在较新版本的 Express 中 这不再有效 var app express 我应该打电
  • 产品口味本地化

    我有 3 种产品风格 每种风格都进行了调试和发布 我已经成功地为每种风格赋予了独特的字符串和图标 现在我正在准备本地化 这需要针对每种产品口味使用不同的字符串 这是我当前的文件夹 目录设置 myApp src main res values
  • 如何在运行 CLI 和 Apache2Handler 时将系统环境变量获取到 PHP 中?

    我的系统是Ubuntu我已经设置了我的环境变量 etc environment 如果我在跑步PHP脚本使用CLI 环境变量来自 etc environment被认可 但是 如果我去执行PHP脚本通过http domain test php
  • Python电子邮件引用-可打印编码问题

    我使用以下命令从 Gmail 中提取电子邮件 def getMsgs try conn imaplib IMAP4 SSL imap gmail com 993 except print Failed to connect print Is
  • 使用 PyQt4 进行 OpenCV 视频捕获

    我一直在尝试使用 PyQt4 GUI 和 OpenCV 捕获视频 我创建了一些按钮 如 开始 结束 等来控制捕获 开始很好 但我无法停止捕捉 要停止捕获 我需要中断 while 循环开始捕获 下面的功能 我无法实现 现在 结束捕获 毁坏了窗
  • 在贪婪重复中回溯平衡组可能会导致不平衡?

    作为出于此问题的目的而一般编写的示例 我的目的是匹配一些数量a的 那么同样数量的b的 再加上一个b 检查此片段中展示的两种模式 也在 ideone com 上 var r1 new Regex xn a a a
  • 如何使CardView具有可点击和可检查的效果,以及如何使其成为深色主题?

    背景 在引入CardView之前 我在上面做了一些选择器my app模仿卡片 并让用户选择应用程序使用的主题 有些人更喜欢深色主题 问题 我想让它看起来和工作起来更原生 所以我尝试使用 CardView 遗憾的是 我无法理解如何设置 Car
  • Math.cos 不准确

    alert Math cos Math PI 2 为什么结果不精确为零 这是不准确的 还是一些实施错误 Math PI 2是实际值的近似值pi 2 取该近似值的精确余弦值不会产生零 您获得的值是该精确值的近似值 最高可达基础浮点数据类型的精
  • 恢复已删除的 Eclipse 项目

    我想在 eclipse 中创建一个新的 git 存储库 当我删除旧的存储库时 不幸的是整个项目已从工作区中删除了 有什么办法可以恢复项目吗 我将不胜感激你的回答 如果您的工作有另一个存储库 例如中央 git 服务器 另一个开发人员 另一台计
  • 在应用程序之间共享文件

    我可以与另一个应用程序共享一个应用程序相关的数据吗 假设我在 apk2 的 resources raw 文件夹中有一个音乐文件 我可以在 apk1 1 中使用相同的文件吗 thx 如果您的应用程序使用相同的证书进行签名并具有相同的andro
  • 为什么 Vue 无法解析本地主机的图像?

    作为 vue js 的后端 我使用 laravel 端口 8000 在我的数据库中 我有用户及其个人资料照片的名称 this user photo 所以 我想展示这张照片 img alt Profile Photo 当我去http loca
  • MySQL:事务与锁定表

    我对事务与锁定表有点困惑 以确保数据库完整性并确保 SELECT 和 UPDATE 保持同步并且没有其他连接干扰它 我需要 SELECT FROM table WHERE LIMIT 1 if condition passes Update
  • Objective-C:`@synthesize fooBar;` 与`@synthesize fooBar=_fooBar;` [重复]

    这个问题在这里已经有答案了 可能的重复 带下划线前缀的综合属性和变量 这是什么意思 我在代码中见过这两个 有什么不同 synthesize fooBar synthesize fooBar fooBar synthesize fooBar
  • 处理IE浏览器中的ctrl+按键事件

    I m using hotkeys Ctrl key in my flex application getting problem when my app is running in IE when I press Ctrl D im ge
  • 使用 Doctrine 构建通用的 OO ACL

    我正在寻求设计一个以学说为基础的 ACL 系统供我自己使用 尽管我在一些最初的设计考虑因素上遇到了困难 现在我正在考虑根据类和唯一标识符来制作它 并将它们存储在表中 如下所示 Table ACL ResourceClass Resource
  • 子进程打开('source venv/bin/activate'),没有这样的文件?

    我想进入 python 文件中的虚拟环境 但它没有引发这样的文件 import subprocess subprocess Popen source Users XX Desktop mio worker venv bin activate
  • 如何绘制双对数 R 图的线性回归?

    我有以下数据 someFactor 500 x c 1 250 y x 25 someFactor 我以双对数图显示 plot x y log xy 现在我使用线性模型 找出 数据的斜率 model lm log y log x model