0

我使用 aa taus <- c(.05,.25,.5,.75,.95) 生成了散点图和分位数线，我对所有结果感到满意，但似乎我的代码重复了少量。当我创建手动比例时，我必须重复代码部分以清理图例。

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

![在此处输入图像描述][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, .25, .5, .75, .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, .25, .5, .75, .95),
                        values=c(1, 1, 1, 1, 1),
                        name="Quantile Factors",
                        labels=c(" 5th  Percentile", "25th Percentile","Median", "75th Percentile", "95th Percentile")) +
                     scale_linetype_manual(breaks=c(.05, .25, .5, .75, .95),
                        values=c("dotted", "dashed","solid", "dashed", "dotted"), 
                        name="Quantile Factors",
                        labels=c(" 5th  Percentile", "25th Percentile","Median", "75th Percentile", "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, 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=SSA)
gg1

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

gg_Munger <- gg_Munger + scale_colour_manual(breaks=c(.05, .25, .5, .75, .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, .25, .5, .75, .95),
                              values=c(1, 1, 1, 1, 1),
                              name="Quantile Factors",
                              labels=c(" 5th  Percentile", "25th Percentile","Median", "75th Percentile", "95th Percentile")) +
                         scale_linetype_manual(breaks=c(.05, .25, .5, .75, .95),
                              values=c("dotted", "dashed","solid", "dashed", "dotted"), 
                              name="Quantile Factors",
                              labels=c(" 5th  Percentile", "25th Percentile","Median", "75th Percentile", "95th Percentile")) 

gg_Munger



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

0 回答 0