R语言ggplot2包绘制哑铃图

数据和代码获取:请查看主页个人信息!!!

载入R包

复制代码
rm(list=ls())
library(tidyverse)
library(ggalt) # devtools::install_github("hrbrmstr/ggalt")
library(scales)

数据载入

复制代码
health <- read.csv("health.csv")
areas <- read.csv("areas.csv")

数据整理

复制代码
health <-
 health %>%
 mutate(area_id=trunc(area_id)) %>%
 arrange(area_id, pct) %>%
 mutate(year=rep(c("2014", "2013"), 26),
        pct=pct/100) %>%
 left_join(areas, "area_id") %>%
 mutate(area_name=factor(area_name, levels=unique(area_name)))
​
health <- setNames(bind_cols(filter(health, year==2014), filter(health, year==2013))[,c(4,1,5)],
        c("area_name", "pct_2014", "pct_2013"))
  • health 数据集中的 area_id 列四舍五入为整数。

  • 根据 area_idpct 列对数据集进行排序。

  • 在数据集中添加一个 year 列,其值分别为 "2014" 和 "2013",并将 pct 列的值除以 100。

  • areas 数据集与 health 数据集进行左连接,连接键为 area_id 列。

  • area_name 列转换为因子,并按唯一值的顺序设置其水平。

  • health 数据集中的 year 为 2014 和 2013 的数据分别筛选出来,并将它们合并成新的数据集。最后,将列重命名为 "area_name"、"pct_2014" 和 "pct_2013"。

可视化

复制代码
txt_col <- "black"
gg <- ggplot(health, aes(x=pct_2014, xend=pct_2013, y=area_name, group=area_name)) +
 geom_dumbbell(colour="#a3c4dc", size=1.5, colour_xend="#0e668b",
                        dot_guide=TRUE, dot_guide_size=0.15) +
 scale_x_continuous(label=percent) +
 labs(x=NULL, y=NULL) +
 theme_bw() +
 theme(plot.background=element_rect(fill="#f7f7f7"),
       panel.background=element_rect(fill="#f7f7f7"),
       panel.grid.minor=element_blank(),
       panel.grid.major.y=element_blank(),
       panel.grid.major.x=element_line(),
       axis.ticks=element_blank(),
       legend.position="top",
       panel.border=element_blank())
gg
​
ggsave('pic.png', width = 5, height = 5)
  • ggplot() 函数指定了基本图形,并设置了 xxendygroup 的美学映射。

  • geom_dumbbell() 函数创建了 dumbbell 图,用于表示两个时间点之间的变化。参数设置了线的颜色、粗细、点的颜色和指导线的长度。

  • scale_x_continuous() 函数将 x 轴标签格式化为百分比。

  • labs() 函数设置了 x 和 y 轴的标签为 NULL。

  • theme_bw() 函数设置了基本的黑白主题,并通过 theme() 函数进一步调整了图的样式,包括背景颜色、网格线、轴标记和图例位置。

  • ggsave() 函数保存了绘制的图形到名为 "pic.png" 的文件中,指定了图形的宽度和高度。

相关推荐
赵钰老师1 小时前
【结构方程模型SEM】最新基于R语言结构方程模型分析
开发语言·数据分析·r语言
vastsmile2 小时前
(R)26.04.23 hermes agent执行本地命令超级慢的原因
开发语言·elasticsearch·r语言
生信小窝14 小时前
081B 基于R包-Ecospat的生态位重叠分化迁移分析与可视化制图(单物种和多物种版本)【2027】
人工智能·python·r语言
数智工坊1 天前
R-CNN目标检测算法精读全解
网络·人工智能·深度学习·算法·目标检测·r语言·cnn
数智工坊1 天前
Faster R-CNN 全精读:实时目标检测的里程碑之作
网络·人工智能·深度学习·目标检测·r语言·cnn
新缸中之脑3 天前
TurboOCR:优化杰作
开发语言·r语言
花间相见3 天前
【MS-Swift实战】:LoRA原理+核心参数(r/alpha)调参指南(适配Qwen-1.8B医疗场景)
开发语言·r语言·swift
叛逆的小小黄3 天前
maxent建模结果中响应曲线的美化
经验分享·笔记·r语言·maxent
sghuter4 天前
AI赋能CI/CD:Gemini实战脚本生成
开发语言·人工智能·ci/cd·青少年编程·r语言
cqbzcsq4 天前
从RNA-seq原始数据开始数据分析(Salmon、tximport基因表达矩阵、DESeq2差异表达、WGCNA共表达网络)
数据挖掘·r语言·生物信息学