基于2个条件ggplot2格式化线型

Question

我正在为 2 个站点绘制基于相互移植（原点 -> 新，原点 -> 原点）的值。我想根据历史记录为这些值设置线型格式，其中 origin -> origin 是实线，但 origin -> new 是虚线。对于 2 个站点，"origin" 和 "new" 的值不同。以下是数据的一个子集：

structure(list(FragID = c("1004", "1006", "1015", "1038", "1087", 
"1089", "1107", "1116"), ParentID = c("166", "166", "166", "166", 
"144", "144", "144", "144"), ParentSite.x = c("Inner Bay", "Inner 
Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay", "Outer 
Bay", "Outer Bay"), FragSite.x = c("Inner Bay", "Outer Bay", "Inner 
Bay", "Outer Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer 
Bay"), TotalSA = c(0.0171970755726674, 0.0338197513082082, 
0.0215722601402604, 0.030712272182997, 0.027529366126288, 
0.029650482611575, 0.0316984120058258, 0.0356299659679559), T03SA = 
c(0.000709604810935872, 0.00148788795124323, 0.00109901406229665, 
0.000966259734683879, 0.000701528253168926, 0.000828107993427705, 
0.00079488114602085, 0.000904344998291552)), .Names = c("FragID", 
"ParentID", "ParentSite.x", "FragSite.x", "TotalSA", "T03SA"), 
class = "data.frame", row.names = c(2L, 3L, 5L, 9L, 21L, 22L, 28L, 
29L))

我想要的是内湾 -> 内湾和外湾 -> 外湾为实线，内湾 -> 外湾和外湾 -> 内湾为虚线。到目前为止我所拥有的如下：

GPBWM03 = ggplot(M3, aes(x = TotalSA, y = T03SA, 
                 group = interaction(FragSite.x, ParentSite.x), color = FragSite.x, 
                 linetype = ParentSite.x)) + 
  geom_point(alpha = 0.5, stroke = 0) +
  geom_smooth(method = "lm", formula = y~x, se = FALSE, fullrange = 
TRUE) + 
  scale_color_manual(values = c("Inner Bay" = "coral2", 
"Outer Bay" = "skyblue4"))

Answer 1

这就是你想要的吗？

ggplot(M3,
       aes(
           x = TotalSA,
           y = T03SA,
           color = FragSite.x,
           linetype = interaction(ParentSite.x, FragSite.x)
       )) +
    geom_smooth(
        method = "lm",
        formula = y ~ x,
        se = FALSE,
        fullrange =
            TRUE
    ) +
    scale_linetype_manual(values = c(1, 2, 2, 1)) +
    scale_color_manual(values =
                           c("Inner Bay" = "coral2", "Outer Bay" = "skyblue4"))

要获得如下输出：

这是使用下面的数据集制作的，因为上面原始 dput 输出中的某些值有换行符：

structure(list(FragID = c("1004", "1006", "1015", "1038", "1087", 
"1089", "1107", "1116"), ParentID = c("166", "166", "166", "166", 
"144", "144", "144", "144"), ParentSite.x = c("Inner Bay", "Inner Bay", 
"Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay", "Outer Bay", 
"Outer Bay"), FragSite.x = c("Inner Bay", "Outer Bay", "Inner Bay", 
"Outer Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay"
), TotalSA = c(0.0171970755726674, 0.0338197513082082, 0.0215722601402604, 
0.030712272182997, 0.027529366126288, 0.029650482611575, 0.0316984120058258, 
0.0356299659679559), T03SA = c(0.000709604810935872, 0.00148788795124323, 
0.00109901406229665, 0.000966259734683879, 0.000701528253168926, 
0.000828107993427705, 0.00079488114602085, 0.000904344998291552
)), class = "data.frame", row.names = c(2L, 3L, 5L, 9L, 21L, 
22L, 28L, 29L))

编辑：对于组合图例，除了线型之外，手动分配颜色可能是最简单的方法：

ggplot(M3,
       aes(
           x = TotalSA,
           y = T03SA,
           color = interaction(ParentSite.x, FragSite.x),
           linetype = interaction(ParentSite.x, FragSite.x)
       )) +
    geom_smooth(
        method = "lm",
        formula = y ~ x,
        se = FALSE,
        fullrange =
            TRUE
    ) +
    scale_linetype_manual(name = "Combination", values = c(1, 2, 2, 1)) +
    scale_color_manual(name = "Combination", values = c("coral2", "coral2", "skyblue4", "skyblue4"))

输出：

或者，您可能会发现按照此处其他一些答案的建议创建交互变量更容易。

Answer 2

另一种方法是创建一个额外的变量来编码 linetype。结果是您得到一个独立显示两个因素的图例：

M3 <- structure(list(FragID = c("1004", "1006", "1015", "1038", "1087", "1089", "1107", "1116"), 
                     ParentID = c("166", "166", "166", "166", "144", "144", "144", "144"), 
                     ParentSite.x = c("Inner Bay", "Inner Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay", "Outer Bay", "Outer Bay"), 
                     FragSite.x = c("Inner Bay", "Outer Bay", "Inner Bay", "Outer Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay"), 
                     TotalSA = c(0.0171970755726674, 0.0338197513082082, 0.0215722601402604, 0.030712272182997, 0.027529366126288, 0.029650482611575, 0.0316984120058258, 0.0356299659679559), 
                     T03SA = c(0.000709604810935872, 0.00148788795124323, 0.00109901406229665, 0.000966259734683879, 0.000701528253168926, 0.000828107993427705, 0.00079488114602085, 0.000904344998291552)), 
                     .Names = c("FragID", "ParentID", "ParentSite.x", "FragSite.x", "TotalSA", "T03SA"), 
                  class = "data.frame", row.names = c(2L, 3L, 5L, 9L, 21L, 22L, 28L, 29L))
M3$in.out <- interaction(M3$ParentSite.x, M3$FragSite.x)
M3$comb <- 'same'
M3$comb[which(M3$ParentSite.x!=M3$FragSite.x)] <- 'different'
GPBWM03 = ggplot(M3, aes(x = TotalSA, y = T03SA, 
                         group = interaction(FragSite.x, ParentSite.x), color = FragSite.x, 
                         linetype = comb)) + 
  geom_point(alpha = 0.5, stroke = 0) +
  geom_smooth(method = "lm", formula = y~x, se = FALSE, fullrange = 
                TRUE) + 
  scale_color_manual(values = c("Inner Bay" = "coral2", 
                                "Outer Bay" = "skyblue4"))
GPBWM03

Answer 3

您的 dput 格式不正确。我使用以下方法手动创建它：

df <- data.frame(FragID = c("1004", "1006", "1015", "1038", "1087", "1089", "1107", "1116"),
                 ParentID = c("166", "166", "166", "166", "144", "144", "144", "144"),
                 ParentSite.x = c("Inner Bay", "Inner Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay", "Outer Bay", "Outer Bay"),
                 FragSite.x = c("Inner Bay", "Outer Bay", "Inner Bay", "Outer Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay"),
                 TotalSA = c(0.0171970755726674, 0.0338197513082082, 0.0215722601402604, 0.030712272182997, 0.027529366126288, 0.029650482611575, 0.0316984120058258, 0.0356299659679559),
                 T03SA = c(0.000709604810935872, 0.00148788795124323, 0.00109901406229665, 0.000966259734683879, 0.000701528253168926, 0.000828107993427705, 0.00079488114602085, 0.000904344998291552))

然后，我在您的数据中添加了一个新列，它使用 tidyverse 包根据 ParentSite.x 和 FragSite.x 的交互确定行类型：

df <- df %>%
  mutate(linetype = ifelse(ParentSite.x=='Inner Bay'&FragSite.x=='Inner Bay' | ParentSite.x=='Outer Bay'&FragSite.x=='Outer Bay'
                             , 'Matching', 'Not Matching'))

然后，我们绘制：

ggplot(df, aes(x = TotalSA, y = T03SA, color = FragSite.x, linetype = linetype)) + 
  geom_point(alpha = 0.5, stroke = 0) +
  geom_smooth(method = "lm", formula = y~x, se = FALSE, fullrange = TRUE) + 
  scale_color_manual(values = c("Inner Bay" = "coral2", 
                                "Outer Bay" = "skyblue4")) +
  labs(linetype='Line Type')

基于2个条件ggplot2格式化线型

Format linetype based on 2 conditions ggplot2

grouping

r

line

ggplot2