使用 ggplot2 进行分面 qqplots

2024-06-28

假设我有以下数据:

datapoints1 = data.frame(categ=c(rep(1, n), rep(2, n)), vals1=c(rt(n, 1, 2), rnorm(n, 3, 4)))
datapoints2 = data.frame(categ=c(rep(1, n), rep(2, n)), vals2=c(rt(n, 5, 6), rnorm(n, 7, 8)))

Using ggplot2,我该如何使用facet在单个命令中创建两个 QQplot 的功能,即一个包含两个t样本,另一个有两个高斯样本?


首先,合并两个数据框:

dat <- cbind(datapoints1, vals2 = datapoints2[ , 2])

然后,对数据进行排序:

dat_sort <- do.call("rbind", lapply(unique(dat$categ), FUN = function(x) {data.frame(categ = x, vals1 = sort(dat$vals1[dat$categ == x]), vals2 = sort(dat$vals2[dat$categ == x]))}))

如果两个样本向量的长度相同,那就很简单:

ggplot() + 
 geom_point(data = dat_sort, aes(x = vals1, y = vals2)) +
 facet_wrap( ~ categ, scales = "free")

一个例子n = 1000:

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

使用 ggplot2 进行分面 qqplots 的相关文章

  • knitr pandoc:“无法使用 pdf writer 生成 pdf 输出”

    Up front using pandoc 在knitr中 当尝试将 md或 Rmd编译为PDF时 它会抱怨 我正在简化流程可重复的研究 正如许多地方所记录的那样 我在用着pandoc and knitr并制作出色的文档 我还试图为一些不太
  • R 绘制一些 unicode 字符,但不绘制其他字符

    我们的系统管理员刚刚将我们的操作系统升级到 SLES12SP1 我重新安装了 Rv3 2 3 并尝试绘图 我用cairo pdf并尝试绘制 x 标签为的图 u0298即太阳能符号 但它不起作用 标签只是空白 例如 cairo pdf Rpl
  • 如何根据子字符串匹配选择 R data.table 行(类似于 SQL)

    我有一个带有字符列的 data table 并且只想选择其中包含子字符串的行 相当于SQLWHERE x LIKE substring E g gt Months data table Name month name Number 1 12
  • igraph - 将文本添加到社区图

    使用下面的代码 我可以在 igraph 提供的 空手道俱乐部 示例中生成社区图 请参见第一张图片 library igraph karate lt make graph Zachary wc lt cluster walktrap kara
  • 如何使用 ggplot2 在轴标签中使用上标

    如何在x轴上打印埃平方 我尝试如下 labs x x axis 2 y y axis 我们可以用bquote library ggplot2 ggplot mtcars aes hp mpg geom point labs x bquote
  • 从数据/系数创建 lm 对象

    有谁知道可以在给定数据集和系数的情况下创建 lm 对象的函数 我对此很感兴趣 因为我开始使用贝叶斯模型平均 BMA 并且我希望能够根据 bicreg 的结果创建 lm 对象 我想访问所有不错的通用 lm 函数 例如诊断绘图 预测 cv lm
  • 在 R 中调用“boot”时没有数据

    我有一个包含很多列的数据框 第一列包含 系统 1 系统 2 等类别 第二列包含代表 0 和 1 的数字 请看下面 例如 SYSTEM Q1 Q2 S1 0 1 S1 1 0 S2 1 1 S2 0 0 S2 1 1 我在 R 中有这段代码来
  • R 矩阵到犰狳的转换非常慢

    观察 对于中等大小的矩阵 将矩阵从 R 传递到 C 的开销要慢得多arma mat类型比NumericMatrix类型 大约需要 250 倍的时间 这是一个最小的例子 include
  • 使用同一变量的多个子集创建新数据框

    我想创建一个新的数据框 其中的列是由不同变量分割的同一变量的子集 例如 我想创建一个新的变量子集 b 其中列由不同变量 year 的子集分割 set seed 88 df lt data frame year rep 1996 1998 3
  • 根据值绘制具有条件颜色的折线图

    我想绘制折线图 根据值 它应该改变它的颜色 我发现的是 plot sin seq from 1 to 10 by 0 1 type p col ifelse sin seq from 1 to 10 by 0 1 gt 0 5 red ye
  • 为什么 NaN^0 == 1

    受到早期高尔夫代码的提示 为什么会 gt NaN 0 1 1 这非常有道理NA 0为 1 因为NA缺少数据 并且any数字提高到 0 将得到 1 包括 Inf and Inf 然而NaN应该代表非数字 那么为什么会这样呢 当帮助页面出现时
  • 是否有一个 R 函数可以获取 n 个对象取 k P(n,k) 的排列数?

    或者我必须给予 P nk lt factorial n factorial n k or P nk lt choose n k factorial k 谢谢 我不知道任何现有的功能 如果 n 较大 您的第一个建议将会失败 当编写为函数时 你
  • 在数据框中使用 Ifelse

    我正在使用的数据框是 gt df lt data frame Name c Joy Jane Jack Jad M1 c 10 40 55 90 gt df Name M1 1 Joy 10 2 Jane 40 3 Jack 55 4 Ja
  • 张量流:RStudio 中的 [NOT FOUND] 错误

    我尝试在中运行以下代码RStudio library tensorflow x data lt runif 100 min 0 max 1 y data lt x data 0 1 0 3 W lt tf Variable tf rando
  • 将嵌套 for 循环转换为 R 中的并行循环

    下面您可以在 R 中找到一段代码 我想将其转换为使用多个 CPU 作为并行进程运行 我尝试使用foreach包 但并没有走得太远 考虑到我有 3 级嵌套循环 我找不到一个很好的例子如何让它工作 我们将非常感谢您的帮助 下面的代码示例 我做了
  • 如何在 R 中计算带有变量的表达式?

    我希望这段代码能够设置plt等于 10 gt var plt gt eval paste0 var lt 10 1 plt lt 10 但相反 它返回一个字符串 I tried eval as expression paste0 var l
  • 绘图中的文本错误 | Jupyter-R

    我正在尝试在通过 Anaconda 安装的 Jupyter 笔记本上使用 R 中的 ggplot2 绘制一个带有一些数据的简单散点图 以下是我使用的代码 ggplot data df1 1 1000 geom point mapping a
  • 在 R 中将数组转换为矩阵

    我有一个数组 其中包括名为 comp 的项目 是 否 的两个熟练度变量 theta0 theta1 这需要转换为一个矩阵 有什么方法可以转换像底部那样的矩阵吗 我的数组如下所示 gt priCPT i6 comp Yes theta1 th
  • 如何引用基于 data.frame 中的变量的列表?

    我有一张简单的桌子emp id and job code 我想返回正确的payout基于job code 我已经用嵌套的 ifelse 来解决这个问题 但是如果我有更多怎么办job code s library dplyr set seed
  • update() 具有局部协变量的函数内的模型

    我需要从函数内部更新回归模型 理想情况下 该函数应该适用于任何类型的模型 lm glm multinom clm 更准确地说 我需要添加函数内部定义的一个或多个协变量 这是一个例子 MyUpdate lt function model ra

随机推荐

  • 如何按组计算日期之间的时间差

    我有一个包含日期 时间和位置的数据框 我想计算组内记录与上一条记录 根据日期排列 之间的分钟差异 并变异为新列 我已经弄清楚如何使用循环来完成此操作 但这仅适用于所有组 位置 而且我不确定如何按组执行此操作 fake data set fo
  • 检测最接近的或父级 div

    我有以下脚本 http jsfiddle net oshirowanen pALBV http jsfiddle net oshirowanen pALBV 如何使其更加动态 这样如果屏幕上的按钮数量未知 我就不必手动指定每个按钮的按钮点击
  • 从php表单获取cookie值

    对于用户可以在表单中输入的值 我该如何设置 cookie 值和名称 我该如何在第二页上显示该值 我不能不使用 cookie 所以虽然可能有更聪明的方法来做到这一点 但我只是想知道如何使用 cookie 来做到这一点 谢谢
  • Cassandra 不等于运算符

    向所有 Cassandra 专家提出问题 我有一个包含大约一百万条记录的列族 我想以这样的方式查询这些记录 我应该能够执行Not Equal To种操作 我用谷歌搜索了这个 看来我必须使用某种Map Reduce 有人可以告诉我在这方面有哪
  • 如何更改 Sublime Text 3 中的光标样式

    在某些时候 我的光标从垂直条变成了下划线 我一定是按下了某种神奇的组合键 知道如何改回来吗 Hit Insert That ll change it back You might have noticed that if you went
  • 确保数组属性的每个元素都符合 React 中的自定义形状

    我想确保数组属性的每个元素都符合特定的形状 这与中回答的问题不同用形状反应 proptype 数组 https stackoverflow com questions 32325912 react proptype array with s
  • 在哪里放置rails_admin 视图的部分内容

    我正在尝试使用部分内容来自定义rails admin 中的编辑视图 在 config initializers rails admin rb 中我有 RailsAdmin config do config config model Albu
  • 基于 Spring mvc Java 的配置不起作用。控制台显示没有错误,但我的jsp页面没有显示

    你好 我正在将我的简单演示项目从 bean 配置转换为纯基于 java 的配置 Bean 配置可以很好地创建表等 但我的 java 配置没有显示任何页面 我解决了许多错误 但现在控制台没有显示指定问题的错误 这是我的代码 请找出问题所在 或
  • 为什么从 numpy 数组继承的类的方法返回不同的东西?

    让我们看一下 numpy 数组的一些非常简单的行为 import numpy as np arr np array 1 2 3 4 5 max value arr max returns 5 现在我创建一个继承自的类np ndarray c
  • 使 XSLT 结果文档由“父”XSLT 文档递归处理

    我有文档grandchild xml 其使用granchild xsl处理后的结果必须转发到child xsl 然后它必须最终由parent xsl处理和输出 我尝试过添加一个xml stylesheet元素添加到结果文档中 期望使用引用的
  • 访问 Microsoft Graph 报告时出现 504 超时错误

    我不断地重复504 Gateway Timeout 15 秒后 尝试通过 Microsoft Graph 访问激活报告时 有时 查询会成功 但通常会在 11 12 秒后成功 随后的 CSV 数据下载可能需要近一分钟才能下载 401 个字节
  • 从名称获取 AWS EMR 集群 ID

    AWS CLI 命令aws emr list clusters返回以下 json 有没有一种方法可以通过 bash 或 groovy 使用名称来获取 Id 我不能只使用 ID 因为我正在删除集群并仅使用相同的名称重建它们 所以我知道集群的名
  • RefenceError:UET 未定义

    我很难理解为什么这个 javascript 错误不断出现在我的网站上 而且我无法重现它 我已按照 bing 文档安装了他们的跟踪像素 因此我的标题中包含以下内容
  • 安卓。是否可以在 C/C++ 上编写自定义函数并在 SQL 查询中使用它?

    SQLite3不支持存储过程 是否可以在 C C 上编写自定义函数并在 SQL 查询中使用它 类似的东西 SELECT FROM table WHERE CUST FUNCTION name1 name2 1 or SELECT CUST
  • 连接到嵌入式 Derby 数据库时出现 java.lang.NoClassDefFoundError

    我在嵌入式模式下使用 Apache NetBeans IDE 11 3 带有 JDK 13 和 Derby 连接数据库时出现以下错误 Exception in thread main java lang NoClassDefFoundErr
  • 使用 Spring Boot 外部 log4.properties 文件

    我正在使用可执行 jar 运行我的应用程序 我在 resources 文件夹中有 log4j properties 在产品中 我想覆盖它并将其放在外部目录中 我如何使用 Spring Boot 做到这一点 Dlogging config p
  • 运行 Spring 单元测试的 AOP 问题

    我有一个 Spring Web 应用程序 配置为使用 JDK 代理进行 AOP AOP 注解 例如 Transactional 在接口上声明 而不是在实现类上声明 应用程序本身工作正常 但当我运行单元测试时 它似乎尝试使用 CGLIB 来实
  • Java:最小和最大日期范围

    我从具有存储日期值的数据库中获取一些数据 并让用户选择他们想要查看数据的日期范围 我用于获取这些日期范围的所有代码都有效 除了获取涵盖所有时间的日期范围的方法之外 这将是 Java 处理的最早可能数据的起始值 到最大可能日期的结束值 我的代
  • 谷歌地图移动SDK不再免费了吗?

    根据这个链接 谷歌价格指南 https cloud google com maps platform user guide Google 从 6 月 11 日起更改了价格 这是否意味着在移动设备上使用谷歌地图 SDK 不再免费 是和不是 您
  • 使用 ggplot2 进行分面 qqplots

    假设我有以下数据 datapoints1 data frame categ c rep 1 n rep 2 n vals1 c rt n 1 2 rnorm n 3 4 datapoints2 data frame categ c rep