R levelplot 调整轴

2024-02-09

我想使用“lattice”包的 levelplot 函数将概率分布函数(PDF)绘制为 R 中的热图。我将 PDF 实现为函数,然后使用值范围和外部函数的两个向量生成水平图的矩阵。 我希望轴显示我的问题是无法在两个轴上添加适当间隔的刻度线,分别显示两个实际值范围而不是列数或行数。

# PDF to plot heatmap
P_RCAconst <- function(x,tt,D)
  {
    1/sqrt(2*pi*D*tt)*1/x*exp(-(log(x) - 0.5*D*tt)^2/(2*D*tt))
  }

# value ranges & computation of matrix to plot
tt_log <- seq(-3,3,0.05)
tt <- exp(tt_log)
tt <- c(0,tt)
x <- seq(0,8,0.05)
z <- outer(x,tt,P_RCAconst, D=1.0)
z[,1] <- 0
z[which(x == 1),1] <- 1.5
z[1,] <- 0.1

# plot heatmap using levelplot
require("lattice")
colnames(z) <- round(tt, 2)
rownames(z) <- x
levelplot(z, cex.axis=1.5, cex.lab=1.5, col.regions=colorRampPalette(c("blue", "yellow","red", "black")), at=seq(0,1.9,length=200), xlab="x", ylab="time t", main="PDF P(x,t)")

Without assigning names to columns and rows I receive the following plot where the tick marks are naturally spaced (as used from other R-functions) but the values are the row & column numbers: Without assigning names to columns and rows I receive the following plot where the tick marks are naturally spaced (as used from other R-functions) but the values are the row & column numbers

通过为列和行分配名称,我收到以下图,其中刻度线根本不可读,但至少对应于实际值:

我已经在这个看似微不足道的问题上花费了太多时间,因此非常感谢您的帮助!


也许这个基本的例子可以有所帮助,

d = data.frame(x=rep(seq(0, 10, length=nrow(volcano)), ncol(volcano)), 
               y=rep(seq(0, 100, length=ncol(volcano)), each=nrow(volcano)), 
               z=c(volcano))

library(lattice)

# let lattice pick pretty tick positions
levelplot(z~x*y, data=d)

# specific tick positions and labels
levelplot(z~x*y, data=d, scales=list(x=list(at=c(2, 5, 7, 9), 
                                            labels=letters[1:4])))

# log scale
levelplot(z~x*y, data=d, scales=list(y=list(log=TRUE)))

这一切都在?xyplot,尽管不可否认这是一页很长的文档。

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

R levelplot 调整轴 的相关文章

  • 在shiny中过滤传单地图数据

    我在用传单地图设置这个闪亮的东西时遇到了麻烦 我的原帖 https stackoverflow com questions 50111566 applying leaflet map bounds to filter data within
  • 如何计算R中移动窗口内的平均斜率

    我的数据集包含2个变量y 和 t 05s y 每 05 秒测量一次 我正在尝试计算移动中的平均坡度20秒窗口 即计算第一个 20 秒斜率值后 窗口向前移动一个时间单位 05 秒 并计算下一个 20 秒窗口 在以下位置生成连续 20 秒斜率值
  • 如何使用 usmap 标记数字而不是名称?

    我知道 usmap 有一个选项label in plot usmap 我想标记一些数字 而不是状态名称 我想 usmap 中应该有与州质心坐标相关的数据 但我不知道如何找到它 如果我能得到 坐标然后我可以用它来标记数字geom text 这
  • 绘制点之间的所有线

    我有以下 R 代码 x lt c 0 01848598 0 08052353 0 06741172 0 11652034 y lt c 0 4177541 0 4042247 0 3964025 0 4074685 d lt data fr
  • twitterR 和 ROAuth R 软件包安装

    我在安装 CRAN 上的 twitteR 和 RAOuth 软件包时遇到一些问题 我尝试了几种不同的方法 在 Windows 下使用源代码 在 Ubuntu 下使用 RStudio 我尝试了以下命令 sudo apt get install
  • 从命令行运行 R 代码 (Windows)

    我在名为 analysis r 的文件中有一些 R 代码 我希望能够从命令行 CMD 运行该文件中的代码 而无需通过 R 终端 并且我还希望能够传递参数并在我的代码中使用这些参数 例如就像下面的伪代码 C gt execute r scri
  • 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
  • 在 R 中创建虚拟变量,排除某些情况为 NA

    我的数据看起来像这样 V1 V2 A 0 B 1 C 2 D 3 E 4 F 5 G 9 我想创建一个虚拟变量R where 0 1 1 2 3 4 and NA 0 5 9 应该很简单 有人可以帮忙吗 我们可以转换V2 into a fa
  • 在 R 格子包中微调点图

    我正在尝试为不同的数据集和不同的算法绘制一堆 ROC 区域 我有三个变量 方案 指定所使用的算法 数据集 是正在测试算法的数据集 以及 Area under ROC 我正在 R 中使用lattice库 命令如下 点图 方案 Area und
  • 以引用透明的方式从函数的省略号参数中提取符号

    事情又发生了 我正要按下发布答案按钮的问题被删除了 我正在寻找一种方法来从函数的省略号参数中提取绑定到符号的对象的值以及符号 也就是说 我试图以引用透明的方式从省略号中提取符号 我尝试过使用替代品和lazy dots 但没有成功 funct
  • 在 R 中使用 lapply 绘制多个数据帧

    我正在尝试使用 lapply 函数绘制多个数据帧 每个数据帧一个图 但是尽管有关此主题的所有帖子我都找不到答案 因为我不断收到错误 图的输出列表为空 我的数据结构如下 df1 lt mtcars gt group by cyl gt tal
  • 只读取选定的列

    谁能告诉我如何仅读取下面每年数据的前 6 个月 7 列 例如使用read table Year Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 2009 41 27 25 31 31 39 2
  • 基于时间窗口的不规则时间序列的优化滚动函数

    有没有办法使用 rollapply 来自zoo包或类似的东西 优化功能 rollmean rollmedian等 使用基于时间的窗口计算滚动函数 而不是基于大量观察的函数 我想要的很简单 对于不规则时间序列中的每个元素 我想计算一个具有 N
  • 将阴影区域添加到五分位数之间的直方图中

    All 我有一个包含 2 个直方图的图表 其中我还绘制了代表第 20 40 60 和 80 个百分位数的线条 下面的代码使用虚拟数据重现了类似的图表 data lt rbind data frame x rnorm 1000 0 1 g o
  • 在 r 中的 group_by 之后建模后取消列表列的嵌套

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

    我正在尝试创建一个新变量 其中包含来自两个不同列的字符串值的唯一计数 所以我有这样的东西 例如 A tibble 4 x 2 names partners
  • python 相当于 R 中的 get() (= 使用字符串检索符号的值)

    在 R 中 get s 函数检索名称存储在字符变量 向量 中的符号的值s e g X lt 10 r lt XVI s lt substr r 1 1 X get s 10 取罗马数字的第一个符号r并将其转换为其等效整数 尽管花了一些时间翻
  • applyStrategy 错误

    我是R新手 最近运行后遇到以下错误applyStrategy函数来自quantstrat包裹 Error in eval expr envir enclos object signal not found Error in colnames
  • 如何将plot中的单变量列表图表转换为ggplot2格式?

    我正在搜索 但仍然找不到一个非常简单的问题的答案 我们如何使用 R 中的 ggplot2 生成一个变量的简单线图 我正在分析时间序列数据 并且想要对图表进行更复杂的操作 我认为如果我使用 ggplot2 代替会更好plot It works
  • 如何使用 dplyr 独立过滤每列的行

    我有以下内容 library tidyverse df lt tibble tribble gene colB colC a 1 2 b 2 3 c 3 4 d 1 1 df gt A tibble 4 x 3 gt gene colB c

随机推荐

  • React 类对象方法声明之间的区别?

    我已经看到了在类组件中声明方法的两种方法React Method 1 class someComp extends React Component handleResetShapes e gt e preventDefault this s
  • 布局中的 Zend_Dojo_Form

    我有一个 Zend Dojo Form 我已将其从我的视图 工作正常 移至我的布局 因为它在每个页面上都有用 然而 在布局中 表单不再起作用 没有任何 dijit 元素出现 它的行为就像普通的 HTML 表单一样 这是我的引导程序的相关部分
  • 更新 D3 圆形包布局

    我正在尝试使用我收到的 json 数据动态更新 d3 圆形包布局 每秒我都会调用 d3 json 来获取新的 json 我的实现只是在旧可视化下创建一个新可视化 而不是更新现有可视化 我想动态更新现有布局
  • 具有 .compact 风格的 UIDatePicker 不尊重内容拥抱优先级

    我认为新UIDatePicker iOS 14 与 compact样式无法与自动布局正确配合 制作一个简单的布局UILabel and UIDatePicker哪里的UIDatePicker具有更高的内容优先级应该导致布局UILabel如果
  • Git 预接收钩子

    当您为 git 存储库启用预接收挂钩时 它不需要参数 但对于要更新的 每个引用 它会在标准输入上接收一行格式 SP SP OF 其中 是存储在 ref 中的旧对象名称 是要存储在 ref 中的新对象名称 并且是 ref 的全名 创建新的 r
  • 如何对使用静态类的方法进行单元测试,而该静态类又使用 ConfigurationElementCollection?

    public class ConfigSection ConfigurationSection public static ConfigSection GetConfigSection return ConfigSection System
  • 这个字符串是什么类型? a:1:{s:2:"en";}

    a 1 s 2 en s 39 Sample Text 我的问题是 您对上述格式有了解吗 我想按照上面的格式 如 JSON 存储长站点设置 settings array slug gt test title gt John Doe desc
  • CloudKit:“内部错误”(1/4000); “无法获取签名证书”>

    我在模拟器上运行 该应用程序工作正常 但一旦在 iPhone 设备上 它就会崩溃 DiscoverUserInfo 中发生错误 discoveryUserInfo 是从下面的 Apple CloudKit 示例代码中复制的代码 这是苹果的事
  • Grails 3 schemaExport 包含查找 sitemesh.xml 的 FileNotFoundException 警告

    When 模式导出使用 Grails 3 3 应用程序的 Gradle 执行 但日志中存在以下警告ddl sql被建造 根据 Grails 3 文档 sitemesh xml 已被删除 因此该文件自然不可用 我错过了什么吗 工具和版本 Gr
  • 告诉 SCons 不要自动创建目录?

    我正在尝试让 SCons 查看我需要的 git 存储库 并希望使该存储库保持最新 问题是我必须告诉它 git 存储库包含哪些文件才能在构建中使用它们 如果我这样做 SCons 将在尝试克隆存储库之前创建存储库 例如 假设我想克隆 GStre
  • 我可以在虚拟机中运行 Docker 吗?

    我尝试在虚拟机上运行 Docker Host MacBook VM Parallels Windows 7 并且出现错误 是否可以 如果虚拟机是 Linux 那么您可以毫无问题地执行此操作 在 Linux 上 Docker 本质上是一个运行
  • 在 CMake 中设置 MSVC 运行时

    我按照 CMake 常见问题解答条目中的说明进行操作 如何使用静态运行时构建我的 MSVC 应用程序 https gitlab kitware com cmake community wikis FAQ how can i build my
  • Pandas:读取 CSV:ValueError:无法将字符串转换为浮点数

    我正在尝试使用 pandas read csv 读取大型且复杂的 CSV 文件 确切的命令是 pd read csv 文件名 quotechar low memory True dtype data types usecols 列 true
  • 如何将值从类传递到活动 - Android

    我有一个关于课程 任务 活动的新手问题 我来自 C 所以我不知道这是否是满足我需要的好方法 我创建了一个类 public class UDPServer private MyDatagramReceiver myDatagramReceiv
  • 当我们在 logback-test.xml 中将日志级别设置为“DEBUG”运行时,出现“Java Heap”内存问题

    Getting Java Heap当我们在 logback test xml 中运行日志级别为 DEBUG 的功能时出现内存问题 我们已经运行了大约 400 个功能 你能帮忙解决这个问题吗 它的 Maven 项目 请升级到0 9 0版本
  • HttpRoutePlanner - 它如何与 HTTPS 代理一起工作

    我设置了 HTTPS 代理 以便 HTTP 客户端可以安全地将纯 HTTP 请求发送到代理 例如 客户端可以向代理发送加密的 HTTP GET 请求 代理将删除加密并将纯 HTTP GET 请求发送到最终站点 我了解到这不是常见的设置 只有
  • AspNet Core 集成测试,将参数传递给 WebApplicationFactory

    我有设计为独立的应用程序aspnet core webapi自托管可执行文件 要启动可执行文件 必须将配置文件路径作为命令行参数传递 例如MyServer exe config path to config file 我想通过集成测试来测试
  • 识别属性的值是一个数组

    我有一个 JSON 文件 abn 63119059513 acn 119059513 business structure Private Company ngr number 1231231 cbh number 1231231 main
  • 在 R 中快速重命名相似名称的多个列

    当我创建数据框时 我有一组 43 个变量 列 默认命名为 X1 X2 X3 X46 我需要用 Q 替换 X 在 R 中是否有更快的方法来执行此操作 或者我必须在每个 X 上单独使用 rename 请帮助 我想学习实现这一目标的有效方法 我是
  • R levelplot 调整轴

    我想使用 lattice 包的 levelplot 函数将概率分布函数 PDF 绘制为 R 中的热图 我将 PDF 实现为函数 然后使用值范围和外部函数的两个向量生成水平图的矩阵 我希望轴显示我的问题是无法在两个轴上添加适当间隔的刻度线 分