如何在 ggplot 上的确切位置放置插图并设置 geom_segment 的颜色?

如何解决如何在 ggplot 上的确切位置放置插图并设置 geom_segment 的颜色?

我正在创建一个关于 loess 工作原理的插图。我的两个疑问在这个问题的末尾。一、设置:

library(tidyverse)
data(melanoma,package = "lattice")
mela <- as_tibble(melanoma)
tric = function(x) if_else(abs(x) < 1,(1 - abs(x)^3)^3,0)
scl = function(x) (x - min(x))/(max(x) - min(x))
mela1 <- mela %>%
  slice(1:9) %>%
  mutate(dist = abs(year - year[5]),scaled = scl(dist),weight = tric(scaled)
  )
mod1 <- lm(incidence ~ year,data = mela1,weights = weight) %>%
  augment(.,mela1)
mela2 <- mela %>%
  slice(10:18) %>%
  mutate(dist = abs(year - year[5]),weight = tric(scaled)
  )
mod2 <- lm(incidence ~ year,data = mela2,mela2)
mela3 <- mela %>%
  slice(19:27) %>%
  mutate(dist = abs(year - year[5]),weight = tric(scaled)
  )
mod3 <- lm(incidence ~ year,data = mela3,mela3)
mela4 <- mela %>%
  slice(28:37) %>%
  mutate(dist = abs(year - year[5]),weight = tric(scaled)
  )
mod4 <- lm(incidence ~ year,data = mela4,mela4)

主要情节:

col <- rainbow_hcl(start = 12,4,l = 20)
colB <- rainbow_hcl(start = 12,l = 100)
main <- ggplot(data = mela,aes(x = year,y = incidence)) +
# segment 1
  geom_segment(
    aes(x = 1936,xend = 1944,y = 2.115717,yend = 2.115717)) +
# segment 2
  geom_segment(
    aes(x = 1945,xend = 1953,y = 3.473217,yend = 3.473217)) +
# segment 3
  geom_segment(
    aes(x = 1954,xend = 1962,y = 1.170247,yend = 1.170247)) +
# segment 4
  geom_segment(
    aes(x = 1963,xend = 1972,y = 2.7,yend = 2.7)) +
  geom_point(data = mod1,color = col[1],shape = 1) +
  geom_point(data = mod2,color = col[2],shape = 0) +
  geom_point(data = mod3,color = col[4],shape = 5) +
  geom_point(data = mod4,color = col[3],shape = 2) +
  geom_line(data = mod1,y = .fitted),color = col[1]) +
  geom_line(data = mod2,color = col[2]) +
  geom_line(data = mod3,color = col[4]) +
  geom_line(data = mod4,color = col[3]) +
  scale_x_continuous(breaks = c(1940,1949,1958,1967))

插入

inset1 <- ggplot(data = mod1,y = weight)) +
  geom_line(color = col[1]) +
  geom_area(fill = colB[1]) +
  theme_void()
inset2 <- ggplot(data = mod2,y = weight)) +
  geom_line(color = col[12) +
  geom_area(fill = colB[2]) +
  theme_void()
inset3 <- ggplot(data = mod3,y = weight)) +
  geom_line(color = col[3]) +
  geom_area(fill = colB[3]) +
  theme_void()
inset4 <- ggplot(data = mod4,y = weight)) +
  geom_line(color = col[4]) +
  geom_area(fill = colB[4]) +
  theme_void()

问题 1: 如何放置四个 inset,使权重函数的 y = 0 位于相应 geom_segment 的高度?我希望主图坐标中的插图高度 = 2。

问题 2: 如何将每个段的颜色设置为相应插图的颜色?

解决方法

不确定我是否做对了一切。但我已经尽力了。 (; 你可以大大简化你的代码

  1. ...通过将您的模型数据绑定到一个数据框以及段的数据中。
  2. ...通过一些命名向量和scale_xxx_manual
  3. 映射美学并设置颜色和形状

对于您的插图,无需制作单独的图并尝试将它们放入主图。您可以通过额外的 geom_linegeom_ribbon 简单地添加它们。要获取线段的高度,将线段数据加入模型数据,以便您可以根据线段的 geom_ribbon 值设置 y 的起始值

library(tidyverse)
library(broom)
library(colorspace)

col <- setNames(col,c("mod1","mod2","mod4","mod3"))
colB <- setNames(colB,"mod3"))
shapes <- setNames(c(1,5,2),"mod3","mod4"))

mods <- list(mod1 = mod1,mod2 = mod2,mod3 = mod3,mod4 = mod4) %>% 
  bind_rows(.id = "mod")

# segments data
dseg <- tribble(
  ~mod,~x,~xend,~y,"mod1",1936,1944,2.115717,1945,1953,3.473217,1954,1962,1.170247,1963,1972,2.7,)

main <- ggplot(data = mela,aes(x = year,y = incidence)) +
  geom_segment(data = dseg,aes(x = x,xend = xend,y = y,yend = y,color = mod)) +
  geom_point(data = mods,aes(color = mod,shape = mod)) +
  geom_line(data = mods,y = .fitted,color = mod)) +
  scale_color_manual(values = col) +
  scale_shape_manual(values = shapes) +
  scale_x_continuous(breaks = c(1940,1949,1958,1967)) +
  guides(color = FALSE,shape = FALSE,fill = FALSE)

mods1 <- left_join(mods,select(dseg,mod,y),by = "mod")

# Add insets
main +
  geom_line(data = mods1,y = weight + y,color = mod,group = mod)) +
  geom_ribbon(data = mods1,ymin = y,ymax = weight + y,fill = mod,group = mod)) +
  scale_fill_manual(values = colB)

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


使用本地python环境可以成功执行 import pandas as pd import matplotlib.pyplot as plt # 设置字体 plt.rcParams[&#39;font.sans-serif&#39;] = [&#39;SimHei&#39;] # 能正确显示负号 p
错误1:Request method ‘DELETE‘ not supported 错误还原:controller层有一个接口,访问该接口时报错:Request method ‘DELETE‘ not supported 错误原因:没有接收到前端传入的参数,修改为如下 参考 错误2:cannot r
错误1:启动docker镜像时报错:Error response from daemon: driver failed programming external connectivity on endpoint quirky_allen 解决方法:重启docker -&gt; systemctl r
错误1:private field ‘xxx‘ is never assigned 按Altʾnter快捷键,选择第2项 参考:https://blog.csdn.net/shi_hong_fei_hei/article/details/88814070 错误2:启动时报错,不能找到主启动类 #
报错如下,通过源不能下载,最后警告pip需升级版本 Requirement already satisfied: pip in c:\users\ychen\appdata\local\programs\python\python310\lib\site-packages (22.0.4) Coll
错误1:maven打包报错 错误还原:使用maven打包项目时报错如下 [ERROR] Failed to execute goal org.apache.maven.plugins:maven-resources-plugin:3.2.0:resources (default-resources)
错误1:服务调用时报错 服务消费者模块assess通过openFeign调用服务提供者模块hires 如下为服务提供者模块hires的控制层接口 @RestController @RequestMapping(&quot;/hires&quot;) public class FeignControl
错误1:运行项目后报如下错误 解决方案 报错2:Failed to execute goal org.apache.maven.plugins:maven-compiler-plugin:3.8.1:compile (default-compile) on project sb 解决方案:在pom.
参考 错误原因 过滤器或拦截器在生效时,redisTemplate还没有注入 解决方案:在注入容器时就生效 @Component //项目运行时就注入Spring容器 public class RedisBean { @Resource private RedisTemplate&lt;String
使用vite构建项目报错 C:\Users\ychen\work&gt;npm init @vitejs/app @vitejs/create-app is deprecated, use npm init vite instead C:\Users\ychen\AppData\Local\npm-
参考1 参考2 解决方案 # 点击安装源 协议选择 http:// 路径填写 mirrors.aliyun.com/centos/8.3.2011/BaseOS/x86_64/os URL类型 软件库URL 其他路径 # 版本 7 mirrors.aliyun.com/centos/7/os/x86
报错1 [root@slave1 data_mocker]# kafka-console-consumer.sh --bootstrap-server slave1:9092 --topic topic_db [2023-12-19 18:31:12,770] WARN [Consumer clie
错误1 # 重写数据 hive (edu)&gt; insert overwrite table dwd_trade_cart_add_inc &gt; select data.id, &gt; data.user_id, &gt; data.course_id, &gt; date_format(
错误1 hive (edu)&gt; insert into huanhuan values(1,&#39;haoge&#39;); Query ID = root_20240110071417_fe1517ad-3607-41f4-bdcf-d00b98ac443e Total jobs = 1
报错1:执行到如下就不执行了,没有显示Successfully registered new MBean. [root@slave1 bin]# /usr/local/software/flume-1.9.0/bin/flume-ng agent -n a1 -c /usr/local/softwa
虚拟及没有启动任何服务器查看jps会显示jps,如果没有显示任何东西 [root@slave2 ~]# jps 9647 Jps 解决方案 # 进入/tmp查看 [root@slave1 dfs]# cd /tmp [root@slave1 tmp]# ll 总用量 48 drwxr-xr-x. 2
报错1 hive&gt; show databases; OK Failed with exception java.io.IOException:java.lang.RuntimeException: Error in configuring object Time taken: 0.474 se
报错1 [root@localhost ~]# vim -bash: vim: 未找到命令 安装vim yum -y install vim* # 查看是否安装成功 [root@hadoop01 hadoop]# rpm -qa |grep vim vim-X11-7.4.629-8.el7_9.x
修改hadoop配置 vi /usr/local/software/hadoop-2.9.2/etc/hadoop/yarn-site.xml # 添加如下 &lt;configuration&gt; &lt;property&gt; &lt;name&gt;yarn.nodemanager.res