绘制从一个数据集生成的分位数回归线并将这些线叠加到另一个散点图上，保持源格式

问题描述

我使用 aa taus

第一个问题：如何选择图例中的文本并为我选择的选定变量着色？ 我想把克拉夫岛和斯蒂尔顿湾染成深蓝色...

![在此处输入图片说明][1]

看来我还不能输入图像或数据的 csv...我还是发帖，抱歉缺少数据。

第二个问题：如何绘制第一个图 (gg_taus) 中的 qeom_quantile 回归线并将它们放置在新的散点图上。我当前的代码正在创建一个新的分位数回归，但我想要第一个。

![在此处输入图片说明][2]

代码如下：

    library(tidyverse)
    library(dplyr)
    library(tidyr)
    library(readr)
    library(plyr)
    library(summarytools)
    library(ggplot2)library(quantreg)

First scatter plot

gg <- ggplot(ESB_REF,aes(x=WaterDepth_m,y=S_TMI)) + 
      geom_point(aes(shape=Assessment_Area,))+
      theme(panel.border = element_rect(linetype = "solid",fill = NA),panel.grid.major = element_blank(),panel.grid.minor = element_blank(),panel.background = element_blank()) + 
      scale_y_continuous(limits = c(0,1))+
      scale_x_continuous(breaks = seq(from = 0,to = 3,by = 0.5))+
      xlab("Water Depth (meters)") + ylab("Std.Macroinvertebrete Trimetric Index (sTMI)")+
      labs(shape="Assessment Areas")
gg

Building a plot with selected taus capturing a range of quantiles in steps  
but cant change the legend variable text color

taus <- c(.05,.25,.5,.75,.95)

gg_taus <- gg + 
  geom_quantile(quantiles = taus,aes(linetype = factor(..quantile..),colour = factor(..quantile..),size = factor(..quantile..) ))

gg_taus <- gg_taus + scale_colour_manual(breaks=c(.05,.95),values=c("red","red","black","green","green"),name="Quantile Factors",labels=c(" 5th  Percentile","25th Percentile","Median","75th Percentile","95th Percentile")) +
                     scale_size_manual(breaks=c(.05,values=c(1,1,1),"95th Percentile")) +
                     scale_linetype_manual(breaks=c(.05,values=c("dotted","dashed","solid","dotted"),"95th Percentile")) 
gg_taus <- gg_taus + theme(legend.title = element_text(size = 10),legend.key=element_rect(size=5),legend.text=element_text(aes(color=fill)))
gg_taus

Trying to visualize a specific Assessment Areas observations overlaying the first plots quantile lines for a new data set.

Munger <- ESB_Rebind %>% group_by(SampleID,Assessment_Area) %>% filter(Assessment_Area == "Munger Landing")
glimpse(Munger)
SSA <- paste("SSA #",Munger$Assess_ID [1]) #Creates a label for legend for Assessment ID number

gg1 <- ggplot(Munger,y=S_TMI)) + geom_point(aes(shape=Assessment_Area))+
        theme(panel.border = element_rect(linetype = "solid",panel.background = element_blank()) + 
        scale_y_continuous(limits = c(0,1))+
        scale_x_continuous(breaks = seq(from = 0,by = 0.5))+
        xlab("Water Depth (meters)") + ylab("Std.Macroinvertebrete Trimetric Index (sTMI)")+
        labs(shape=SSA)
gg1

gg_Munger <- gg1 + 
  geom_quantile(quantiles = taus,size = factor(..quantile..)))

gg_Munger <- gg_Munger + scale_colour_manual(breaks=c(.05,"95th Percentile")) +
                         scale_size_manual(breaks=c(.05,"95th Percentile")) +
                         scale_linetype_manual(breaks=c(.05,"95th Percentile")) 

gg_Munger

  [1]: https://i.stack.imgur.com/P1dnf.png
  [2]: https://i.stack.imgur.com/Ua49B.png

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）