R 中的堆叠条形图与线条图

2024-02-22

我想将以下堆积条形图放在一起 与直线曲线。

但为什么下面的代码不起作用呢? 正确的做法是什么?

x<-c(0,1,2,3,4,5,6,7,8,9,10,11,12);

    # To plot line 
    emp_dens <- c(0.107,0.184,0.205,0.185,0.138,0.091,0.049,0.023,0.01,0.0028,0.0012,0.00023,0.00013);
    dat_dens <- as.matrix(cbind(x,emp_dens));


    # To plot stack bar
    dens_all_k <- c(1.15e-01, 1.89e-01, 2.05e-01, 1.82e-01,1.36e-01,8.68e-02,4.71e-02,2.21e-02,9.17e-03,3.37e-03,1.11e-03,3.37e-04,9.31e-05)

    # Each k0..5 compose the stack 
    # For example
    # dens_k0[1] + .... dens_k5[1] ~= dens_all_k[1]

    dens_k0 <-c(2.52e-02,8.38e-02,1.38e-01, 1.53e-01,1.27e-01,8.44e-02, 4.66e-02, 2.21e-02, 9.16e-03, 3.37e-03,1.11e-03, 3.37e-04, 9.31e-05)
    dens_k1 <- c(6.75e-02, 8.91e-02, 5.86e-02, 2.51e-02, 8.59e-03,2.25e-03, 4.90e-04,9.35e-05, 1.55e-05,2.21e-06,2.99e-07, 3.55e-08,3.92e-09)
    dens_k2 <- c(1.70e-02,1.64e-02,7.95e-03, 2.56e-03,6.20e-04,1.20e-04, 1.93e-05, 2.67e-06, 3.23e-07,3.47e-08,3.36e-09, 2.95e-10,2.38e-11)
    dens_k3 <- c(0.005124596,0,0,0,0,0,0, 0, 0, 0, 0,0,0)
    dens_k4 <- c(0.0004270497, 0, 0,0,0, 0, 0, 0, 0,0,0, 0, 0)
    dens_k5 <- c(2.760725e-05, 0, 0, 0,0,0, 0, 0, 0, 0,0, 0,0)


    barplot(cbind(0:max(x),dens_all_k),xlim=c(0,max(x)),ylim=c(0,max(emp_dens)),,space=0.1,lwd=5,xlab="Value of X",ylab="Densities",font.main=1);
    lines(dat_dens,lty=1,col="red");

dens_all_k 仅包含聚合摘要。如果您正在寻找图片中的条形图,您需要提供所有 k 信息。尝试以下代码。

dens_kall<-rbind(dens_k0,dens_k1,dens_k2,dens_k3,dens_k4,dens_k5)
ltext<-c("K0","K1","K2","K3","K4","K5")
colnames(dens_kall)<-0:12
barplot(height=dens_kall,xlim=c(0,max(x)),ylim=c(0,max(emp_dens)),,space=0.1,lwd=5,xlab="Value of X",ylab="Densities",font.main=1
        ,legend.text =ltext,
        args.legend = list(x = "topright")
        );
lines(dat_dens,lty=1,col="red");

The output isenter image description here

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

R 中的堆叠条形图与线条图 的相关文章

  • 为什么 sapply 的缩放速度比样本大小的 for 循环慢?

    假设我想采用向量 X 2 1 N 并将 e 计算为每个元 素的指数 是的 我认识到最好的方法就是通过向量化 exp X 但这样做的目的是将 for 循环与 sapply 进行比较 我通过逐步尝试三种方法 一种使用 for 循环 两种以不同方
  • 更新 R6 对象实例中的方法定义

    如何更新 R6 类实例的方法定义 正如我所期望的 S3 使用当前的方法定义 对于 R5 参考类 我可以使用 myInstance myInstance copy 在 R6 中 我尝试了 myInstance myInstance clone
  • 从 data.frame 在 ggplot 图例中添加信息

    我想在图例中添加信息 哪个传感器具有该值 这是我的代码 z lt data frame a c sensor 1 sensor 2 sensor 3 sensor 4 sensor 5 sensor 6 sensor 7 sensor 8
  • 在 R 中提取 data.frames 列表的名称以及 data.frame 中的值

    在下面的代码中 j是 data frames 的命名列表 我想知道是否有办法 a 提取变量的数值 即one short and one long 在 data frames 内并附加它们的相关名称 即 AAA or BBB or CCC 到
  • data.table 抛出“找不到对象”错误[重复]

    这个问题在这里已经有答案了 我有一个数据表 library data table mydt lt data table index 1 10 当我在全局环境中尝试它时 我可以让它工作 但当我在调试器中或在包测试中使用它时却无法工作 问题是我
  • 为什么这个 R ggplot2 代码会显示一个空白的显示设备?

    虽然 SO 通常不用于帮助解决错误 但这个显示了特别简单且特别烦人的行为 如果你是一个ggplot2用户 您可以在 10 秒或更短的时间内重现它 正如这个 GitHub 问题 ggplot gtable 创建空白显示 https githu
  • 增加雷达图中长轴标签的空间

    我想创建一个雷达图ggirahExtra ggRadar 问题是我的标签很长并且被剪掉了 我想我可以通过添加在标签和绘图之间创建更多空间margin margin 0 0 2 0 cm to element text in axis tex
  • R中IF函数的使用

    我正在短跑ifR 中的函数 但收到以下警告消息 In if runif 50 0 1 lt 0 69 the condition has length gt 1 and only the first element will be used
  • 要在子集中显示的非数字条目的维恩图

    我有以下数据框 SET1 SET2 SET3 par1 par2 par1 par2 par3 par2 par3 par4 par5 我想制作一个维恩图 其中所有这些 parX 元素都显示在各自的子集中 即作为标签 而不仅仅是重叠元素的数
  • 如何根据 ggplot2 中的汇总数据创建堆积条形图

    我正在尝试使用 ggplot 2 创建堆积条形图 我的宽格式数据如下所示 每个单元格中的数字是响应的频率 activity yes no dontknow Social events 27 3 3 Academic skills works
  • case_when 与部分字符串匹配和 contains()

    我正在使用一个数据集 其中有许多名为 status1 status2 等的列 在这些列中 它表示某人是否豁免 完整 注册等 不幸的是 豁免投入并不一致 这是一个示例 library dplyr problem lt tibble perso
  • 如何在 R 中的 dygraph 标题中使用 UTF-8 字符

    使用 Rstudio Windows8 当我使用 dygraph 函数绘制时间序列时 在尝试在主标题中使用 UTF 8 字符时遇到问题 library dygraphs dygraph AirPassengers main T tulo 这
  • 扩展数据框以使其具有与原始行中两列的范围一样多的行[重复]

    这个问题在这里已经有答案了 我有一个数据框如下 structure list symbol c u n v i a start c 9L 6L 10L 8L 7L end c 14L 15L 12L 13L 11L Names c symb
  • 如何将旋转的 NetCDF 转换回正常的纬度/经度网格?

    我有一个带有旋转坐标的 NetCDF 文件 我需要将其转换为正常的纬度 经度坐标 经度为 180到180 纬度为 90到90 library ncdf4 nc open dat nf 对于尺寸 它显示 1 5 variables exclu
  • R - 基于列名称的子集

    我的数据框有超过 120 列 变量 我想根据列名称创建子集 例如 我想创建一个子集 其中列名称包含字符串 心情 这可能吗 我一般用 SubData lt myData grep whatIWant colnames myData 我很清楚
  • 我如何查看 quantmod 包中所有可用的数据系列?

    如何显示可用的所有报价 数据系列的列表 例如使用雅虎的 getSymbols 我不知道有什么办法 TTR包有一个功能 stockSymbols 下载 NYSE AMEX 和 NASDAQ 的所有当前代码 它试图将它们采用雅虎可接受的格式 但
  • 使用data.table进行聚合

    经过 SO 用户的多次建议后 我终于尝试将我的代码转换为使用data table library data table DT lt data table plate paste0 plate rep 1 2 each 5 id rep c
  • 如何在R中分离两个图?

    每当我运行这段代码时 第一个图就会简单地覆盖前一个图 R中有没有办法分开得到两个图 plot pc title main abc xlab xx ylab yy plot pcs title main sdf xlab sdf ylab x
  • 斯皮尔曼相关性和联系

    我正在一小组配对排名上计算斯皮尔曼的 rho 斯皮尔曼因处理领带不当而闻名 例如 取2组8个排名 即使两组中有6个是平局 相关性仍然很高 gt cor test c 1 2 3 4 5 6 7 8 c 0 0 0 0 0 0 7 8 met
  • 闪亮井板宽度

    library shiny library shinydashboard ui lt dashboardPage dashboardHeader dashboardSidebar dashboardBody wellPanel tags d

随机推荐

  • sveltejs 中的响应式全宽画布

    我对 svelte 很陌生 我正在尝试使用 svelte 在全屏上渲染画布 听起来很容易做到 但我无法让它正常工作 我正在绑定一个width and height变量为clientWidth clientHeight父级并使用这些变量来设置
  • VB.NET 中 IsNumeric() 的错误结果

    我在 VB NET 中有一个函数 它循环遍历值并尝试将其转换为小数 如果IsNumeric是真的 Dim Value As String If IsNumeric Value True Then Rate CType Value Decim
  • 在 excepthook 中打印原始异常

    我正在设置系统异常钩子这样我就可以记录发生的每个异常 让我们使用以下示例 而不是写入日志 def excepthook self type value traceback print n print type print value pri
  • 无法在opencv中加载haarcascadeshaarcascade.xml

    我尝试在 android 中使用 opencv 进行人脸检测 但无法加载对象检测 xml 文件 代码如下 static CvHaarClassifierCascade cascade 0 CvMemStorage storage 0 LOG
  • Excel 删除 Python 插入的有效公式

    我正在尝试使用 Python 生成 Excel xlsx 报告openpyxl 但有一个问题 Excel 不断从我的工作表中删除公式在一些特定的细胞中 我打印了正在插入的内容 它似乎绝对有效 如果我将其从控制台粘贴到单元格 它会显示确切的结
  • 我该如何:Git 单向推送/不可编辑的工作目录?

    有没有办法拥有两台机器 假设A和B A 是我的本地计算机 我在其中开发代码并推送它 B 是一个服务器 永远不需要编辑代码 但需要一个工作目录 是否有可能做到这一点without做任何事git pull是来自 B 的什么吗 基本上 我只想以一
  • 为什么 tokio::spawn 即使使用 .clone() 也会抱怨生命周期?

    我试图编译以下看似简单的代码 但出现错误 use std io Error derive Debug struct NetworkConfig bind String node key file String async fn networ
  • 检测视图控制器何时进入后台并恢复

    我正在实现一个具有以下要求的 ViewController 如果用户在 15 分钟后返回应用程序 则视图应重新加载数据 我正在考虑使用viewDidDisappear当应用程序进入后台时保存时间戳viewDidAppear用于检查以前保存的
  • 更新复选框并从颤动对话框中返回值

    我正在尝试将一些城市列表添加到带有复选框的对话框中 以便我需要实现多次单击项目 我正在尝试做的事情如下所示 onPressed从按钮调用 Rest Service 并在成功结果时我只显示一个对话框 void showCityDialog B
  • 其余 api 设计和上传图像的工作流程。

    我想设计一个API 允许客户端上传图像 然后应用程序创建图像的不同变体 例如调整大小或更改图像格式 最后应用程序将每个变体的图像信息存储在数据库中 当我尝试确定执行此任务的正确策略时 就会出现问题 以下是我能想到的一些不同策略 策略一 发送
  • 让实体框架使用 Contains 而不是 Like 并解释“ESCAPE ~”

    我在 EF 中使用了一行 LINQ 它基本上是在做myTable Where c gt c Contains mystring 这是生成的代码 SELECT TOP 300 Extent1 ID AS ID Extent1 FKFishEn
  • Zend Framework 2 的注释命名空间未加载 DoctrineMongoODMModule

    我已经加载了 zf2 的 Doctrine MongoODM 模块 我的控制器中有文档管理器 一切都很顺利 直到我尝试保留文档 它因以下错误而失败 语义错误 SdsCore Document User 类中的注释 Document 从未导入
  • 使用龙目岛项目安全吗? [关闭]

    Closed 这个问题是基于意见的 help closed questions 目前不接受答案 如果你不知道龙目岛计划 http projectlombok org帮助解决 Java 的一些烦恼 比如生成带注释的 getter 和 sett
  • Heroku 应用程序崩溃并显示“libruby.so.1.9:无法打开共享对象文件”[已关闭]

    这个问题不太可能对任何未来的访客有帮助 它只与一个较小的地理区域 一个特定的时间点或一个非常狭窄的情况相关 通常不适用于全世界的互联网受众 为了帮助使这个问题更广泛地适用 访问帮助中心 help reopen questions 我在部署这
  • 如何使用 jquery 添加一个元素始终作为最后一个元素?

    无论如何 我希望将某个 div 添加为列表中的最后一个元素 有什么方法可以明确指定这一点吗 list append div div 会将其附加到 list 的最后 如果你想将它附加到最后一个div 以防万一之后还有其他元素 那么你可以使用
  • 在 matplotlib 中使用 SymLogNorm 标准化格式化颜色条刻度标签

    TL DR 你怎么 强制LogFormatter对每个标签都使用科学记数法 现在它使用它的值小于0或大于1000 似乎并没有暴露任何set powerlimit我也能找到的方法 有什么方法可以做到正确或者您应该使用不同的格式化程序吗 哪一个
  • 如何在mysql的现有表中添加另一个字段

    Okay I have these table 我这里还有另一张桌子 那么 idAlbum 在 gallery 表中 对此的查询是什么 我想从一个表中的一个字段中选择数据并将其插入到另一表中的字段中 如果您想在 PhpMyAdmin 中的现
  • 递归地从同一个表父子表中获取最后一条记录

    我想从同一张表中获取最后一个参考 ID 我有下表 ID UserId DelegatedToUserId 1 100 101 2 101 102 3 102 103 4 103 NULL 5 104 109 我就是无法思考 我知道我想要什么
  • 查找所有显示器的编号和分辨率

    如何轮询窗口来查看连接的显示器以及它们运行的 分辨率 In C Screen Class http msdn microsoft com en us library system windows forms screen aspx表示单个系
  • R 中的堆叠条形图与线条图

    我想将以下堆积条形图放在一起 与直线曲线 但为什么下面的代码不起作用呢 正确的做法是什么 x lt c 0 1 2 3 4 5 6 7 8 9 10 11 12 To plot line emp dens lt c 0 107 0 184