如何在 ggplot2 线图的每个方面显示每个组的误差线？

问题描述 投票：0回答：1

我有一个用 ggplot2 制作的多面折线图。每个面内有 3 条线（3 组）。我想显示每个方面每个组的误差线。但是，我的代码仅显示 1 组错误栏。

数据如下：




structure(list(Ceramide = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 
3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 
5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 
7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 
8L, 8L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 
10L, 10L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 
11L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 13L, 13L, 13L, 
13L, 13L, 13L, 13L, 13L, 13L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 
14L, 14L), levels = c("Cer.d18.0.16.0.", "Cer.d18.0.24.0.", "Cer.d18.1.17.0.", 
"GlcCer..d18.1.16.0.", "GlcCer..d18.1.22.0.", "GlcCer..d18.1.24.0.", 
"SM.d18.0.16.0.", "SM.d18.0.20.0.", "SM.d18.0.24.0.", "SM.d18.1.16.0.", 
"SM.d18.1.18.0.", "SM.d18.1.20.0.", "SM.d18.1.22.0.", "SM.d18.1.24.0."
), class = "factor"), Time = c(1, 1, 1, 2, 2, 2, 3, 3, 3, 1, 
1, 1, 2, 2, 2, 3, 3, 3, 1, 1, 1, 2, 2, 2, 3, 3, 3, 1, 1, 1, 2, 
2, 2, 3, 3, 3, 1, 1, 1, 2, 2, 2, 3, 3, 3, 1, 1, 1, 2, 2, 2, 3, 
3, 3, 1, 1, 1, 2, 2, 2, 3, 3, 3, 1, 1, 1, 2, 2, 2, 3, 3, 3, 1, 
1, 1, 2, 2, 2, 3, 3, 3, 1, 1, 1, 2, 2, 2, 3, 3, 3, 1, 1, 1, 2, 
2, 2, 3, 3, 3, 1, 1, 1, 2, 2, 2, 3, 3, 3, 1, 1, 1, 2, 2, 2, 3, 
3, 3, 1, 1, 1, 2, 2, 2, 3, 3, 3), group = structure(c("control", 
"nonsurgery", "surgery", "control", "nonsurgery", "surgery", 
"control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery", "control", "nonsurgery", "surgery", "control", "nonsurgery", 
"surgery"), levels = c("Control", "Non-Surgery", "Surgery")), 
    lipid_median = c(31.7097370252311, 34.8841161859164, 33.8290414320265, 
    34.1684491300412, 36.8599396212937, 34.6874485703282, 32.3304714410658, 
    34.3821632026185, 34.2918789867913, 145.421759786139, 163.919371708403, 
    161.414926600578, 168.052324203937, 171.299868576961, 163.482362347776, 
    156.007891775535, 150.109050938528, 140.413393151251, 7.55427982744622, 
    8.43892179554658, 8.84364339049424, 8.47618281296365, 9.73184893405914, 
    8.89943502138335, 9.54992271330453, 9.75591139150991, 8.33992912914349, 
    120.522851205812, 125.88352056183, 124.742784442186, 135.854812388421, 
    147.013486685589, 138.741281304854, 175.417239189346, 174.897434598814, 
    187.270501509993, 677.051591528187, 713.879764664933, 736.210940972631, 
    750.318094316966, 908.009357805281, 891.9355482138, 998.724623171574, 
    1049.42761667314, 1077.18803324428, 515.501136153635, 547.04680741506, 
    573.986096454809, 606.695091894266, 705.473041362558, 696.376555407304, 
    865.287364243197, 934.24691964498, 1041.58738466609, 2396.61854413292, 
    2537.68682151989, 2655.14871236125, 2709.59299639615, 2836.2032465437, 
    2740.33497762081, 2891.72482041317, 3120.80299514295, 3398.57801419228, 
    1192.04543834071, 1325.15708062365, 1373.82781244999, 1261.73343893034, 
    1493.43715612378, 1378.26319014475, 1308.33341650225, 1427.098659029, 
    1383.38634626637, 225.48902982964, 255.589558412548, 260.787456931915, 
    239.491154054235, 264.145991461667, 258.183045212106, 243.190264905069, 
    257.541556100897, 245.264138149509, 44846.5182860006, 46612.102474005, 
    49407.1947133826, 46418.2518567308, 50920.7013400283, 50964.4113345973, 
    47738.1122862251, 50623.1346625701, 53886.4226383792, 10121.4255466185, 
    10560.678173424, 11812.8685609723, 10248.8308646972, 11459.0944364258, 
    10906.0994579433, 10056.9557976663, 10349.2740180237, 9543.47413288215, 
    22840.2447470509, 24353.1897501268, 26494.3889699438, 23295.5202105924, 
    29569.7631259859, 26427.4951377757, 28332.8398285893, 29242.663675407, 
    28274.2067377423, 48590.5678321935, 57003.9794389598, 61764.7881690795, 
    52978.3341523694, 64576.4593418158, 58564.0252854514, 58321.7642765169, 
    62054.6822882946, 58373.9112859722, 21062.0272182675, 23075.0210421309, 
    25978.4554297391, 23613.9666659656, 28988.7090450242, 26758.4962229301, 
    26935.1615084363, 28585.1403999489, 28525.4576752008), error25_mean = c(24.2346184329791, 
    24.2346184329791, 24.2346184329791, 25.5670968115505, 25.5670968115505, 
    25.5670968115505, 25.096215055222, 25.096215055222, 25.096215055222, 
    105.383596099945, 105.383596099945, 105.383596099945, 111.642319476418, 
    111.642319476418, 111.642319476418, 100.924841961843, 100.924841961843, 
    100.924841961843, 4.69249765133767, 4.69249765133767, 4.69249765133767, 
    5.53075253153509, 5.53075253153509, 5.53075253153509, 5.44700156122991, 
    5.44700156122991, 5.44700156122991, 90.5082978297349, 90.5082978297349, 
    90.5082978297349, 102.435152667828, 102.435152667828, 102.435152667828, 
    128.990824500558, 128.990824500558, 128.990824500558, 366.618976742406, 
    366.618976742406, 366.618976742406, 515.694055014973, 515.694055014973, 
    515.694055014973, 584.752767225474, 584.752767225474, 584.752767225474, 
    297.695004557019, 297.695004557019, 297.695004557019, 417.989154819201, 
    417.989154819201, 417.989154819201, 533.575947699868, 533.575947699868, 
    533.575947699868, 1809.98995236381, 1809.98995236381, 1809.98995236381, 
    2184.20172139014, 2184.20172139014, 2184.20172139014, 2282.69728959617, 
    2282.69728959617, 2282.69728959617, 627.061827023858, 627.061827023858, 
    627.061827023858, 792.813146312815, 792.813146312815, 792.813146312815, 
    758.242434200142, 758.242434200142, 758.242434200142, 125.307690779964, 
    125.307690779964, 125.307690779964, 146.508496378177, 146.508496378177, 
    146.508496378177, 144.668297639273, 144.668297639273, 144.668297639273, 
    33773.788835109, 33773.788835109, 33773.788835109, 39348.7616985581, 
    39348.7616985581, 39348.7616985581, 37887.3127288562, 37887.3127288562, 
    37887.3127288562, 7665.54019614451, 7665.54019614451, 7665.54019614451, 
    8343.47012458038, 8343.47012458038, 8343.47012458038, 7632.22950854432, 
    7632.22950854432, 7632.22950854432, 12817.5016527568, 12817.5016527568, 
    12817.5016527568, 16840.1422237024, 16840.1422237024, 16840.1422237024, 
    16712.3475713006, 16712.3475713006, 16712.3475713006, 27344.676729014, 
    27344.676729014, 27344.676729014, 33688.9651052161, 33688.9651052161, 
    33688.9651052161, 30563.0891504596, 30563.0891504596, 30563.0891504596, 
    11222.4607168074, 11222.4607168074, 11222.4607168074, 14474.5924596798, 
    14474.5924596798, 14474.5924596798, 14996.2737491839, 14996.2737491839, 
    14996.2737491839), error75_mean = c(45.6274436242186, 45.6274436242186, 
    45.6274436242186, 46.3250271711141, 46.3250271711141, 46.3250271711141, 
    44.1410863901294, 44.1410863901294, 44.1410863901294, 257.298206727319, 
    257.298206727319, 257.298206727319, 236.928413792139, 236.928413792139, 
    236.928413792139, 222.166200366427, 222.166200366427, 222.166200366427, 
    13.9818618325458, 13.9818618325458, 13.9818618325458, 14.737541922349, 
    14.737541922349, 14.737541922349, 14.681766170529, 14.681766170529, 
    14.681766170529, 171.639325994704, 171.639325994704, 171.639325994704, 
    192.425598936517, 192.425598936517, 192.425598936517, 235.054195358118, 
    235.054195358118, 235.054195358118, 1345.92244450824, 1345.92244450824, 
    1345.92244450824, 1446.37417417621, 1446.37417417621, 1446.37417417621, 
    1690.12813085014, 1690.12813085014, 1690.12813085014, 988.599076943615, 
    988.599076943615, 988.599076943615, 1124.21640197936, 1124.21640197936, 
    1124.21640197936, 1421.89923746621, 1421.89923746621, 1421.89923746621, 
    3583.04960825947, 3583.04960825947, 3583.04960825947, 3565.4733498257, 
    3565.4733498257, 3565.4733498257, 4101.06275926889, 4101.06275926889, 
    4101.06275926889, 2612.34329875199, 2612.34329875199, 2612.34329875199, 
    2433.16211146142, 2433.16211146142, 2433.16211146142, 2263.42116616642, 
    2263.42116616642, 2263.42116616642, 461.090471271548, 461.090471271548, 
    461.090471271548, 442.036576511789, 442.036576511789, 442.036576511789, 
    413.885341711507, 413.885341711507, 413.885341711507, 64516.9309659138, 
    64516.9309659138, 64516.9309659138, 63219.8911218361, 63219.8911218361, 
    63219.8911218361, 65114.6601617267, 65114.6601617267, 65114.6601617267, 
    15112.1061771317, 15112.1061771317, 15112.1061771317, 14380.35473491, 
    14380.35473491, 14380.35473491, 12917.5242138469, 12917.5242138469, 
    12917.5242138469, 40726.4653152597, 40726.4653152597, 40726.4653152597, 
    42029.1504387542, 42029.1504387542, 42029.1504387542, 44774.7182898133, 
    44774.7182898133, 44774.7182898133, 100670.488661911, 100670.488661911, 
    100670.488661911, 98911.2606696033, 98911.2606696033, 98911.2606696033, 
    96209.4844634763, 96209.4844634763, 96209.4844634763, 49892.4612000841, 
    49892.4612000841, 49892.4612000841, 49214.3230551524, 49214.3230551524, 
    49214.3230551524, 49239.3503559168, 49239.3503559168, 49239.3503559168
    )), row.names = c(NA, -126L), class = c("tbl_df", "tbl", 
"data.frame"))

这是我的代码：


# Define custom labels for Ceramide names to use in the graph
ceramide_labels <- c("Cer.d18.0.16.0." = "Cer(d18:0/16:0)",
                     "Cer.d18.0.24.0." = "Cer(d18:0/24:0)",
                     "Cer.d18.1.17.0." = "Cer(d18:1/17:0)",
                     "GlcCer..d18.1.16.0." = "GlcCer (d18:1/16:0)",
                     "GlcCer..d18.1.22.0." = "GlcCer (d18:1/22:0)",
                     "GlcCer..d18.1.24.0." = "GlcCer (d18:1/24:0)",
                     "SM.d18.0.16.0."="SM(d18:0/16:0)",
                     "SM.d18.0.20.0." = "SM(d18:0/20:0)",
                     "SM.d18.0.24.0."= "SM(d18:0/24:0)",
                     "SM.d18.1.16.0." = "SM(d18:1/16:0)",
                     "SM.d18.1.18.0." = "SM(d18:1/18:0)",
                     "SM.d18.1.20.0." = "SM(d18:1/20:0)",
                     "SM.d18.1.22.0." = "SM(d18:1/22:0)",
                     "SM.d18.1.24.0." = "SM(d18:1/24:0)")

# Creating a multifaceted line plot of the nominally significant ceramides
plot<-graphs%>%
  ggplot(aes(x = Time, y = lipid_median)) +
  geom_line(aes(color = group, group = group), linewidth = 0.5) + # Set line width
  geom_point(aes(shape = group), size = 2.5) + # Set point color
  facet_wrap(~ Ceramide, scales = "free_y", labeller = labeller(Ceramide = ceramide_labels), 
              nrow = 4) + # Adjust the height and width as needed) +
  geom_errorbar(aes(ymin = error25_mean, ymax = error75_mean, color=group), width = 0.2) +  # Add error bars
  labs(x = "Time Point", 
       y = "Sphingolipid Value", 
       title = "Nominally Significant Sphingolipid Means over Time",
       color = "Group") +
  theme_minimal() +
  guides(color = "none", shape = guide_legend(title = "Group")) +  
  scale_x_continuous(breaks = c(1, 2, 3), labels = c("0", "2 Years", "12 Years")) +
  scale_y_continuous(expand = c(0, 0)) +
  theme(plot.margin = margin(20, 30, 20, 30, "pt"), 
        panel.spacing = unit(0.5, "lines"),  # Adjust panel spacing
        #axis.ticks.x = element_blank(),  # Remove x-axis ticks
        axis.text.x = element_text(color = "black", size = 6.5),  # Keep x-axis labels
        strip.text = element_text(size = 8),
        strip.background = element_rect(fill='beige'),  
        panel.background = element_rect(fill = 'white',color='black'),
        plot.background = element_rect(fill = 'white', color = 'white'),
        panel.border =element_rect(colour='black', fill=NA),
        plot.title = element_text(hjust = 0.5)) #+  # Center the plot title

如何才能显示所有 3 组的误差线？我希望误差线与不同组的线条颜色相匹配

r ggplot2
1个回答
0
投票

您需要躲避误差线，可能还有点。

graphs %>%
  ggplot(aes(x = Time, y = lipid_median)) +
  geom_line(aes(color = group, group = group), linewidth = 0.5) + # Set line width
  geom_point(aes(shape = group), size = 2.5,
             position = position_dodge(width=0.2)) + # Set point color
  facet_wrap(~ Ceramide, scales = "free_y", labeller = labeller(Ceramide = ceramide_labels), 
             nrow = 4) + # Adjust the height and width as needed) +
  geom_errorbar(aes(ymin = error25_mean, ymax = error75_mean, color=group), 
                width = 0.2,
                position = position_dodge(width=0.2)) + ... # Add error bars

最新问题
© www.soinside.com 2019 - 2024. All rights reserved.