将点连接到图中的回归线

2023-11-21

考虑以下数据ISLR 书:

Income <- structure(list(X = 1:30, Education = c(10, 10.4013377926421, 
10.8428093645485, 11.2441471571906, 11.6454849498328, 12.0869565217391, 
12.4882943143813, 12.8896321070234, 13.2909698996656, 13.7324414715719, 
14.133779264214, 14.5351170568562, 14.9765886287625, 15.3779264214047, 
15.7792642140468, 16.2207357859532, 16.6220735785953, 17.0234113712375, 
17.4648829431438, 17.866220735786, 18.2675585284281, 18.7090301003344, 
19.1103678929766, 19.5117056856187, 19.9130434782609, 20.3545150501672, 
20.7558528428094, 21.1571906354515, 21.5986622073579, 22), Income = c(26.6588387834389, 
27.3064353457772, 22.1324101716143, 21.1698405046065, 15.1926335164307, 
26.3989510407284, 17.435306578572, 25.5078852305278, 36.884594694235, 
39.666108747637, 34.3962805641312, 41.4979935356871, 44.9815748660704, 
47.039595257834, 48.2525782901863, 57.0342513373801, 51.4909192102538, 
61.3366205527288, 57.581988179306, 68.5537140185881, 64.310925303692, 
68.9590086393083, 74.6146392793647, 71.8671953042483, 76.098135379724, 
75.77521802986, 72.4860553152424, 77.3550205741877, 72.1187904524136, 
80.2605705009016)), .Names = c("X", "Education", "Income"), class = "data.frame", row.names = c(NA, 
-30L))

我想重现他们的图,它将点连接到回归线来表示误差项,即:

enter image description here

很容易绘制点和回归线ggplot2:

ggplot(Income, aes(Education, Income)) + geom_point(color="red")+geom_smooth(se=FALSE)

但我想不出像他们那样将点连接到回归线的简单方法。

有没有一种简单的方法可以做到这一点(无论是base, ggplot2 or lattice)?


初步代码

require("ggplot2")

mod <- loess(Income ~ Education, data = Income)
Income <- transform(Income, Fitted = fitted(mod))

ggplot 版本

ggplot(Income, aes(Education, Income)) + 
  geom_point(color="red") + 
  geom_smooth(se=FALSE, method = "loess") +
  geom_segment(aes(x = Education, y = Income,
                   xend = Education, yend = Fitted))

基础图形版本

plot(Income ~ Education, data = Income, type = "p", col = "red",
     cex = 1.25)
points(Fitted ~ Education, data = Income)
lines(Fitted ~ Education, data = Income, col = "blue")
with(Income, segments(Education, Income, Education, Fitted))

相应调整以满足您的审美需求。

enter image description here

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

将点连接到图中的回归线 的相关文章

  • 从 R 中的向量中选择所有可能的元组

    我正在尝试用 R 编写一个程序 当给定一个向量时 将返回所有可能的tuples http en wikipedia org wiki Tuples该向量中的元素 例如 元组 c a b c c a b c 出租车 c a c c b c c
  • R::bigmemory - 如何创建角色big.matrix?

    我尝试使用bigmemory封装在R我一开始就陷入困境 我愿意 temp lt matrix paste a 1 10 5 2 并得到一个字符矩阵 没关系 但后来我尝试 x lt as big matrix temp type char 我
  • .wav 文件长度/持续时间,无需读入文件

    有没有办法提取有关 wav 文件长度 持续时间的信息 而无需在 R 中读取文件 我有数千个这样的文件 如果我必须阅读每个文件才能找到其持续时间 那将需要很长时间 Windows 文件资源管理器为您提供了打开 长度 字段的选项 并且您可以查看
  • 将年月格式转换为 POSIXct [重复]

    这个问题在这里已经有答案了 我有一些年月形式的数据 我想将其格式化以用于绘图ggplot date lt c 2016 03 2016 04 2016 05 2016 06 2016 07 2016 08 2016 09 2016 10 2
  • LDA with topicmodels,如何查看不同文档属于哪些主题?

    我正在使用 topicmodels 包中的 LDA 我已经在大约 30 000 个文档上运行它 获取了 30 个主题 并获得了主题的前 10 个单词 它们看起来非常好 但我想看看哪些文档属于哪个主题的概率最高 我该怎么做 myCorpus
  • 是否可以通过扫描从控制台读取而不回显字符?

    这是一个示例函数 passwordEntry lt function cat Enter your password pwd lt scan n 1 what character quiet TRUE invisible pwd 并测试该功
  • 使用 broom 和 tidyverse 总结 r 平方游戏

    我发布了一个问题here https stackoverflow com questions 48627287 getting adjusted r squared value for each line in a geom smooth
  • 如何从 R 运行带有特定模块的 perl 脚本?

    我可以从终端运行 perl 脚本 myperlscript pl 没有任何问题 但是 如果我尝试从 RStudio 中运行相同的 perl 脚本 则会出现以下错误 command lt myperlscript pl outputfile
  • 删除ggplot2中的负图区域[重复]

    这个问题在这里已经有答案了 如何删除 ggplot2 中 x 轴和 y 轴下方的绘图区域 请参见下面的示例 我尝试了几个主题元素 panel border panel margin plot margin 但没有任何运气 p lt ggpl
  • 计算每个唯一值出现的次数

    假设我有 v rep c 1 2 2 2 25 现在 我想计算每个唯一值出现的次数 unique v 返回唯一值是什么 但不返回它们的数量 gt unique v 1 1 2 我想要一些能给我的东西 length v v 1 1 25 le
  • R中的一元加/减是什么?

    来自 R 的详细信息部分Syntax http stat ethz ch R manual R patched library base html Syntax html帮助页面 定义了以下一元和二元运算符 他们被列出 在优先级组中 从最高
  • 我无法下载 R 中的 reshape2 包 [关闭]

    Closed 这个问题是无法重现或由拼写错误引起 help closed questions 目前不接受答案 我在尝试安装 R 包时收到此响应 gt installed packages reshape2 Package LibPath V
  • 使用 pracma::findpeaks 识别持续峰值

    我的语法有问题peakpat内的选项findpeaks内的函数pramcaR 包 v 2 1 1 我使用的是 R 3 4 3 x64 Windows 我希望该函数能够识别可能有两个重复值的峰值 并且我相信该选项peakpat这就是我能做到的
  • 通过间接引用列来修改数据框中的某些值

    我正在整理一些数据 我们将失败的数据分类到垃圾箱中 并按批次计算每个分类箱的有限产量 我有一个描述排序箱的元表 这些行按升序测试顺序排列 一些排序标签带有非语法名称 sort tbl lt tibble tribble weight lab
  • API 请求和curl::curl_fetch_memory(url, handle = handle) 中的错误:SSL 证书问题:证书已过期

    几天前 我运行了代码几个月 没有任何问题 GET url myurl query 今天我遇到一个错误 Error in curl curl fetch memory url handle handle SSL certificate pro
  • 将数据框中重叠的范围合并到唯一的组中

    我有一个 n 行 3 的数据框 df lt data frame start c 178 400 983 1932 33653 end c 5025 5025 5535 6918 38197 group c 1 1 2 2 3 df sta
  • 在 r 中的 group_by 之后建模后取消列表列的嵌套

    我想对所有组进行线性回归group by 将模型系数保存在列表列中 然后使用 unnest 扩展列表列 这里我用的是mtcars以数据集为例 注 我想用do here becausebroom tidy 不适用于所有型号 mtcars gt
  • 相当于 min() 的 rowMeans()

    我在 R 邮件列表上多次看到这个问题 但仍然找不到满意的答案 假设我有一个矩阵m m lt matrix rnorm 10000000 ncol 10 我可以通过以下方式获得每行的平均值 system time rowMeans m use
  • 旋转 Markdown 的表格 pdf 输出

    我想将 pdf 上的表格输出旋转 90 度 我正在使用 Markdown 生成报告并kable循环显示表格 如果可以的话我想继续使用kable因为还有很多其他依赖于它的东西我没有包含在这个 MWE 中 这是一个简单的例子 使用iris数据集
  • 有效地绘制大时间序列(matplotlib)

    我正在尝试使用 matplotlib 在同一轴上绘制三个时间序列 每个时间序列有 10 6 个数据点 虽然生成图形没有问题 但 PDF 输出很大 在查看器中打开速度非常慢 除了以栅格化格式工作或仅绘制时间序列的子集之外 还有其他方法可以获得

随机推荐

  • 使用浏览器作为桌面 UI

    如何使用浏览器作为桌面应用程序的 UI 到目前为止我想出的方法是 使用所有 HTML Javascript 问题 无法访问文件系统或其他任何内容 在使用应用程序时运行本地网络服务器 问题 用户完成后如何杀死它 我的用户的技术水平不足以使用
  • 通过 .net 会员提供商进行编程登录

    我正在尝试对一段需要当前登录用户进行测试的代码进行单元测试 使用 Net 2 0 Membership Provider 我如何以编程方式作为用户登录以进行此测试 if Membership ValidateUser user1 P ssw
  • 更改 jQuery Mobile 的默认数据主题

    通过 jQuery Mobile 我可以使用自定义主题创建页面 div 现在这可以工作 但要求我在每个页面中以及每次添加新页面时添加此行 我尝试添加data theme s 到 body 标签 但这没有影响 除了每页手动设置之外 还有什么方
  • 确定 matplotlib 轴大小(以像素为单位)

    给定一组axes在 matplotlib 中 有没有办法确定其大小 以像素为单位 我需要根据调整更大或更小的数字来缩放事物 特别是我想更改线宽 使其与轴尺寸成比例 这给出了以英寸为单位的宽度和高度 bbox ax get window ex
  • 如何从对象数组列表中找到最大元素?

    Collections max arraylist 不起作用 并且常规for循环也不起作用 我所拥有的是 ArrayList
  • 保护我的代码免遭逆向工程

    正如类似问题中所讨论的here and here我想保护我的代码免受逆向工程的影响 我的情况是这样的Simucal在他的 优秀 答案中描述here 基本上 归结为 你唯一的机会是 源盗窃的目标是如果您 有一些非常具体 难以 工程师 与您相关
  • set 和 get 方法与公共变量的优点[重复]

    这个问题在这里已经有答案了 可能的重复 为什么要使用 getter 和 setter 在类中创建方法来访问私有变量而不是使变量成为公共变量有什么好处吗 例如 第二种情况比第一种情况更好吗 Case 1 public class Shoe p
  • Rails:HABTM - 查找所有没有关联的记录

    我有 2 个模型 锻炼 设备 它们属于多个关系 如果我使用Workout find all joins gt equipment conditions gt equipment id 5 它有效 但如果我使用Workout find all
  • 根据高度保持div长宽比

    我需要维护width元素占其元素的百分比height 因此 随着高度的变化 宽度也会更新 相反 可以通过使用 padding top 的 值来实现 但 padding left 作为百分比将是对象宽度的百分比 而不是其高度 因此 使用这样的
  • 高效附加到可变长度字符串容器(Golang)

    问题 我需要将多个正则表达式应用于大日志文件的每一行 例如几GB长 收集非空匹配并将它们全部放入一个数组中 用于序列化并通过网络发送 如果回答以下问题 切片没有多大帮助这个问题 holds 如果切片没有足够的容量 append 将需要分配新
  • 在 SQL Server 查询中将 NULL 替换为 0

    我开发了一个查询 在前三列的结果中我得到NULL 我怎样才能替换它0 Select c rundate sum case when c runstatus Succeeded then 1 end as Succeeded sum case
  • Java SSL:如何禁用主机名验证

    标准 java SSL 套接字是否可以通过属性禁用 ssl 连接的主机名验证 到目前为止我发现的唯一方法是编写一个始终返回 true 的主机名验证器 Weblogic 提供了这种可能性 可以使用以下属性禁用主机名验证 Dweblogic s
  • 仅允许使用“sudo su - user”的 Capistrano

    我正在尝试对远程服务器进行简单的 Rails Capistrano 部署 不幸的是我无法得到sudo开箱即可正确运行 我需要在这里部署 drwxr xr x 2 user www 4096 Sep 28 15 05 my app dir 并
  • 空闲 PostgreSQL 连接是否有超时?

    1 S postgres 5038 876 0 80 0 11962 sk wai 09 57 00 00 00 postgres postgres my app 1 45035 idle 1 S postgres 9796 876 0 8
  • Python 中的模块有标准别名吗?

    遵循中提出的指南这个帖子 我正在改变所有 from module import function function agt by import module as mdl mdl function agt 在我的代码中 我尝试使用常用的别名
  • 有没有办法在不支持的浏览器中渲染 SVG 图标?

    截至目前 似乎唯一支持它们的浏览器是 Firefox 显然 Opera 曾经支持它 但后来放弃了它 也许是 JavaScript 垫片 重用过程中与序列化无关的部分如何将 SVG 元素设置到页面的图标 创建一个
  • 不区分大小写的字符串比较

    我想比较两个变量以查看它们是否相同 但我希望此比较不区分大小写 例如 这将区分大小写 if var1 var2 但我希望它不区分大小写 我该如何处理呢 这相当简单 你只需要打电话strtolower 在两个变量上 如果需要处理 Unicod
  • Mono 缺少 System.Numerics.BigInteger.Parse(string,IFormatProvider)

    我目前使用以下错误Json NET Newtonsoft JSON图书馆 Missing method System Numerics BigInteger Parse string IFormatProvider in assembly
  • 名称中包含破折号的自定义 AngularJS 指令不起作用

    我编写了以下 Angular 指令 它将向所有子项添加 required 属性 directive requireall function compile return restrict A only want it triggered f
  • 将点连接到图中的回归线

    考虑以下数据ISLR 书 Income lt structure list X 1 30 Education c 10 10 4013377926421 10 8428093645485 11 2441471571906 11 645484