ggplot2 中两个不同组的不同调色板

2023-11-24

我正在尝试创建一个图，显示多个地点的观测数据和建模数据的每月土壤湿度垂直剖面。

到目前为止，我只能绘制一组值，无论是观察值还是建模值，如下例所示：

library(ggplot2)
library(RColorBrewer)

# Create customized color palette
mypal <- colorRampPalette(brewer.pal(6,"PuBu"))

ggplot(df1, aes(x=value, y=depth, colour=as.factor(month))) +
  geom_path() +
  facet_wrap(~ site) +
  scale_y_reverse() +
  scale_colour_manual(values=mypal(12)) +
  theme_bw(base_size=18) +
  ylab("Depth") + xlab(bquote('Soil moisture (' ~m^3~m^-3*')')) +
  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())

这是重现它的数据：

df1 <- structure(list(site = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 
6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 
6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 
6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 
6L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 
7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 
7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 
7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 
8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 
8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 
8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 
8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 
9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 
9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 
9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 
9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 
10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 11L, 
11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 
11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 
11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 
11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 
11L, 11L, 11L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L, 12L, 12L, 
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 
12L, 12L), .Label = c("IL_Shabbona_5_NNE", "ME_Limestone_4_NNW", 
"ME_Old_Town_2_W", "MI_Chatham_1_SE", "MI_Gaylord_9_SSW", "MN_Goodridge_12_NNW", 
"MN_Sandstone_6_W", "NH_Durham_2_N", "NH_Durham_2_SSW", "NY_Ithaca_13_E", 
"NY_Millbrook_3_W", "WI_Necedah_5_WNW"), class = "factor"), month = c(1L, 
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 
5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 
8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 
11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 
1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 
4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 
7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 
10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 
6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 
6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 
6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 
6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 
6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 
6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 
6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 
6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 
6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 
6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L), depth = c(5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 10, 10, 10, 
10, 10, 10, 10, 10, 10, 10, 10, 10, 20, 20, 20, 20, 20, 20, 20, 
20, 20, 20, 20, 20, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 
50, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 10, 10, 10, 10, 10, 10, 10, 
10, 10, 10, 10, 10, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 
20, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 100, 100, 
100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 5, 5, 5, 5, 
5, 5, 5, 5, 5, 5, 5, 5, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 
10, 10, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 50, 50, 
50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 100, 100, 100, 100, 100, 
100, 100, 100, 100, 100, 100, 100, 5, 5, 5, 5, 5, 5, 5, 5, 5, 
5, 5, 5, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 20, 
20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 50, 50, 50, 50, 50, 
50, 50, 50, 50, 50, 50, 50, 100, 100, 100, 100, 100, 100, 100, 
100, 100, 100, 100, 100, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 
10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 20, 20, 20, 20, 
20, 20, 20, 20, 20, 20, 20, 20, 50, 50, 50, 50, 50, 50, 50, 50, 
50, 50, 50, 50, 100, 100, 100, 100, 100, 100, 100, 100, 100, 
100, 100, 100, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 10, 10, 10, 
10, 10, 10, 10, 10, 10, 10, 10, 10, 20, 20, 20, 20, 20, 20, 20, 
20, 20, 20, 20, 20, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 
50, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 10, 10, 10, 10, 10, 10, 10, 
10, 10, 10, 10, 10, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 
20, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 100, 100, 
100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 5, 5, 5, 5, 
5, 5, 5, 5, 5, 5, 5, 5, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 
10, 10, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 50, 50, 
50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 100, 100, 100, 100, 100, 
100, 100, 100, 100, 100, 100, 100, 5, 5, 5, 5, 5, 5, 5, 5, 5, 
5, 5, 5, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 20, 
20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 50, 50, 50, 50, 50, 
50, 50, 50, 50, 50, 50, 50, 100, 100, 100, 100, 100, 100, 100, 
100, 100, 100, 100, 100, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 
10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 20, 20, 20, 20, 
20, 20, 20, 20, 20, 20, 20, 20, 50, 50, 50, 50, 50, 50, 50, 50, 
50, 50, 50, 50, 100, 100, 100, 100, 100, 100, 100, 100, 100, 
100, 100, 100, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 10, 10, 10, 
10, 10, 10, 10, 10, 10, 10, 10, 10, 20, 20, 20, 20, 20, 20, 20, 
20, 20, 20, 20, 20, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 
50, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 10, 10, 10, 10, 10, 10, 10, 
10, 10, 10, 10, 10, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 
20, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 100, 100, 
100, 100, 100, 100, 100, 100, 100, 100, 100, 100), clm = c(0.324, 
0.321, 0.334, 0.33, 0.285, 0.277, 0.264, 0.261, 0.261, 0.273, 
0.303, 0.339, 0.348, 0.349, 0.356, 0.333, 0.289, 0.281, 0.269, 
0.267, 0.266, 0.277, 0.303, 0.352, 0.345, 0.349, 0.355, 0.337, 
0.295, 0.288, 0.277, 0.275, 0.274, 0.284, 0.299, 0.344, 0.299, 
0.306, 0.314, 0.331, 0.312, 0.305, 0.295, 0.293, 0.292, 0.3, 
0.299, 0.295, 0.274, 0.278, 0.284, 0.308, 0.302, 0.298, 0.289, 
0.287, 0.286, 0.292, 0.291, 0.281, 0.477, 0.475, 0.475, 0.518, 
0.379, 0.353, 0.352, 0.341, 0.347, 0.373, 0.46, 0.499, 0.399, 
0.401, 0.399, 0.435, 0.346, 0.326, 0.325, 0.316, 0.32, 0.339, 
0.37, 0.409, 0.339, 0.34, 0.337, 0.362, 0.313, 0.294, 0.293, 
0.286, 0.289, 0.306, 0.317, 0.345, 0.256, 0.259, 0.26, 0.274, 
0.281, 0.262, 0.26, 0.254, 0.255, 0.271, 0.263, 0.256, 0.243, 
0.242, 0.242, 0.252, 0.281, 0.267, 0.266, 0.261, 0.26, 0.274, 
0.27, 0.25, 0.497, 0.467, 0.461, 0.434, 0.377, 0.361, 0.341, 
0.314, 0.334, 0.372, 0.394, 0.514, 0.363, 0.35, 0.344, 0.366, 
0.34, 0.327, 0.311, 0.289, 0.304, 0.336, 0.344, 0.367, 0.308, 
0.294, 0.286, 0.318, 0.304, 0.293, 0.279, 0.26, 0.272, 0.3, 0.309, 
0.308, 0.232, 0.228, 0.237, 0.281, 0.271, 0.263, 0.25, 0.235, 
0.239, 0.265, 0.277, 0.25, 0.235, 0.233, 0.24, 0.28, 0.273, 0.265, 
0.254, 0.242, 0.241, 0.264, 0.277, 0.258, 0.448, 0.446, 0.454, 
0.49, 0.341, 0.31, 0.297, 0.298, 0.325, 0.359, 0.427, 0.478, 
0.351, 0.352, 0.359, 0.403, 0.312, 0.287, 0.276, 0.276, 0.298, 
0.325, 0.333, 0.355, 0.305, 0.309, 0.312, 0.338, 0.281, 0.26, 
0.25, 0.249, 0.266, 0.29, 0.288, 0.3, 0.233, 0.233, 0.237, 0.265, 
0.256, 0.238, 0.229, 0.227, 0.236, 0.256, 0.252, 0.239, 0.231, 
0.229, 0.228, 0.253, 0.262, 0.247, 0.237, 0.235, 0.24, 0.255, 
0.257, 0.24, 0.326, 0.323, 0.341, 0.343, 0.268, 0.236, 0.209, 
0.213, 0.228, 0.27, 0.301, 0.34, 0.311, 0.31, 0.321, 0.322, 0.254, 
0.228, 0.206, 0.208, 0.22, 0.254, 0.278, 0.317, 0.285, 0.285, 
0.295, 0.302, 0.243, 0.22, 0.2, 0.2, 0.211, 0.241, 0.253, 0.283, 
0.218, 0.22, 0.222, 0.259, 0.236, 0.215, 0.197, 0.193, 0.198, 
0.225, 0.232, 0.22, 0.19, 0.188, 0.187, 0.214, 0.225, 0.205, 
0.189, 0.182, 0.183, 0.201, 0.213, 0.2, 0.437, 0.435, 0.437, 
0.456, 0.391, 0.367, 0.358, 0.339, 0.347, 0.374, 0.468, 0.468, 
0.402, 0.401, 0.406, 0.442, 0.39, 0.369, 0.362, 0.345, 0.351, 
0.364, 0.41, 0.425, 0.386, 0.386, 0.389, 0.415, 0.375, 0.355, 
0.349, 0.336, 0.34, 0.349, 0.377, 0.401, 0.318, 0.319, 0.32, 
0.343, 0.341, 0.32, 0.316, 0.307, 0.307, 0.313, 0.31, 0.321, 
0.302, 0.308, 0.312, 0.323, 0.335, 0.313, 0.309, 0.304, 0.303, 
0.307, 0.303, 0.299, 0.39, 0.387, 0.391, 0.405, 0.347, 0.335, 
0.328, 0.317, 0.327, 0.339, 0.387, 0.411, 0.387, 0.386, 0.39, 
0.392, 0.336, 0.325, 0.319, 0.31, 0.318, 0.326, 0.367, 0.399, 
0.367, 0.368, 0.373, 0.38, 0.332, 0.322, 0.317, 0.308, 0.315, 
0.321, 0.342, 0.374, 0.299, 0.302, 0.305, 0.337, 0.318, 0.307, 
0.304, 0.297, 0.3, 0.306, 0.298, 0.3, 0.265, 0.27, 0.274, 0.293, 
0.293, 0.281, 0.28, 0.275, 0.276, 0.281, 0.275, 0.265, 0.588, 
0.591, 0.555, 0.362, 0.367, 0.355, 0.317, 0.228, 0.161, 0.28, 
0.381, 0.412, 0.38, 0.389, 0.399, 0.301, 0.304, 0.294, 0.265, 
0.199, 0.149, 0.234, 0.315, 0.303, 0.302, 0.296, 0.323, 0.277, 
0.279, 0.27, 0.245, 0.19, 0.148, 0.211, 0.288, 0.262, 0.173, 
0.168, 0.188, 0.202, 0.196, 0.184, 0.171, 0.134, 0.113, 0.116, 
0.192, 0.184, 0.151, 0.136, 0.142, 0.18, 0.176, 0.164, 0.155, 
0.131, 0.119, 0.115, 0.165, 0.166, 0.588, 0.591, 0.555, 0.362, 
0.367, 0.355, 0.317, 0.228, 0.161, 0.28, 0.381, 0.412, 0.38, 
0.389, 0.399, 0.301, 0.304, 0.294, 0.265, 0.199, 0.149, 0.234, 
0.315, 0.303, 0.302, 0.296, 0.323, 0.277, 0.279, 0.27, 0.245, 
0.19, 0.148, 0.211, 0.288, 0.262, 0.173, 0.168, 0.188, 0.202, 
0.196, 0.184, 0.171, 0.134, 0.113, 0.116, 0.192, 0.184, 0.151, 
0.136, 0.142, 0.18, 0.176, 0.164, 0.155, 0.131, 0.119, 0.115, 
0.165, 0.166, 0.472, 0.46, 0.457, 0.442, 0.387, 0.382, 0.371, 
0.37, 0.377, 0.394, 0.432, 0.497, 0.452, 0.45, 0.45, 0.416, 0.367, 
0.362, 0.353, 0.352, 0.358, 0.372, 0.39, 0.453, 0.414, 0.418, 
0.42, 0.4, 0.354, 0.35, 0.342, 0.34, 0.346, 0.359, 0.363, 0.408, 
0.331, 0.341, 0.345, 0.362, 0.336, 0.332, 0.325, 0.322, 0.327, 
0.341, 0.337, 0.327, 0.296, 0.301, 0.306, 0.329, 0.326, 0.321, 
0.314, 0.311, 0.314, 0.328, 0.327, 0.309, 0.451, 0.44, 0.458, 
0.441, 0.344, 0.315, 0.303, 0.314, 0.324, 0.349, 0.396, 0.47, 
0.334, 0.327, 0.349, 0.372, 0.302, 0.279, 0.269, 0.278, 0.285, 
0.305, 0.31, 0.34, 0.287, 0.284, 0.293, 0.33, 0.279, 0.259, 0.248, 
0.256, 0.263, 0.281, 0.281, 0.288, 0.181, 0.182, 0.181, 0.226, 
0.212, 0.194, 0.184, 0.185, 0.192, 0.208, 0.209, 0.188, 0.174, 
0.175, 0.177, 0.208, 0.209, 0.193, 0.185, 0.184, 0.189, 0.201, 
0.203, 0.186, 0.365, 0.364, 0.367, 0.372, 0.309, 0.303, 0.288, 
0.292, 0.299, 0.308, 0.351, 0.382, 0.371, 0.372, 0.379, 0.376, 
0.311, 0.305, 0.291, 0.295, 0.301, 0.309, 0.343, 0.383, 0.365, 
0.366, 0.372, 0.376, 0.318, 0.313, 0.3, 0.303, 0.309, 0.316, 
0.33, 0.371, 0.33, 0.338, 0.342, 0.359, 0.33, 0.325, 0.313, 0.315, 
0.32, 0.327, 0.317, 0.324, 0.291, 0.296, 0.3, 0.319, 0.317, 0.312, 
0.303, 0.304, 0.307, 0.314, 0.304, 0.293)), row.names = c(NA, 
-720L), class = "data.frame", .Names = c("site", "month", "depth", 
"value"))

然而，我想在同一个图中再添加一个“组”。最终的绘图将具有额外的垂直线“集”，但使用不同的调色板（例如使用红色而不是蓝色）。

我尝试创建此图的失败尝试使用以下代码：

# Create some random data based on the first dataset
df1 <- cbind(df1, type='observed')
df2 <- df1[1:4]
df2$value <- df2$value * 1.5
df2 <- cbind(df2, type='modeled')
df3 <- rbind(df1,df2)

# Try to plot it
ggplot(df3, aes(x=value, y=depth, group=type, colour=as.factor(month))) +
  geom_path() +
  facet_wrap(~ site) +
  scale_y_reverse() +
  theme_bw(base_size=18) +
  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())

这显然没有显示出我想要的。为什么每个面只显示两条线？其他垂直线在哪里？如何为这两个组配置不同的调色板？

在 ggplot2 中如何做到这一点？

您可以使用interaction结合type and month: color = interaction(as.factor(month), type))。代替：group=type, colour=as.factor(month).

要创建红色和蓝色托盘，请使用两个mypal功能：

mypal <- colorRampPalette(brewer.pal(6, "PuBu"))
mypal2 <- colorRampPalette(brewer.pal(6, "YlOrRd"))

Code:

library(ggplot2)
library(RColorBrewer)

mypal <- colorRampPalette(brewer.pal(6, "PuBu"))
mypal2 <- colorRampPalette(brewer.pal(6, "YlOrRd"))

ggplot(df3, 
       aes(value, depth, color = interaction(as.factor(month), type))) +
    geom_path() +
    facet_wrap(~ site) +
    labs(title = "Soil moisture by depth and site",
         subtitle = "Observed and expected data",
         x = bquote('Soil moisture (' ~m^3~m^-3*')'),
         y = "Depth") +
    scale_y_reverse() +
    scale_colour_manual(values = c(mypal(12), mypal2(12))) +
    theme_classic() + 
    theme(legend.position = "none")

Plot:

本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)

r

plot

ggplot2

Colors

ggplot2 中两个不同组的不同调色板的相关文章

列值的切换功能

我有一个缩写变量名称错误地分散在整个列表中请参阅下面的示例结构 ID lt c SPW SM DLS SJ joe schmoe CEJ teddy roos GVF MJC LH sally fields Full names sho
将儒略日期转换为 PosixCt 日期

我发现自己在解决这个问题我需要将 R 中的儒略日期转换为正常日期 YYYY MM DD 我知道我可以指定as Date julian date origin 但我不知道应该提供哪个来源我的朱利安日期类似于 2458010 2458011
重叠并固定小平面中的 x 尺度（尺度=“自由”）以制作电影

在另一篇文章中 Justin 的帮助下我用以下方法绘制了模拟结果facet gridggplot2 中的选项包含数据和答案的帖子在这里使用facet grid选项通过ggplot2绘制数据框的列 https stackoverflow
使用 data() 的 R 包命名空间问题 -- 找不到数据集

我在尝试在我自己的包中导入包即 robfilter 时遇到了问题我尝试使用它的方法之一 adore filter 在这一行失败 data critvals 出现错误未找到数据集 critvals 如果我通过 require robfi
Sweave + RweaveHTML：cat 输出未出现在输出中

我对 Sweave RweaveHTML 有疑问我希望 cat 的输出最终出现在正在生成的 html 文件中我有一个案例它没有我不明白为什么 test function bla bla cat Result is 然后在 Rnw 文
如何合并具有相同列名的数据框

我有一个数据框如下所示 structure list Variables structure list ADA ADA LEAD LEAD BIG4 BIG4 LOGMKT LOGMKT LEV LEV ROA ROA ROAL ROAL
在前两个冒号上分割字符串

我想在前两个冒号上拆分一列字符串但不在任何后续冒号上拆分 my data lt read table text my string some data 123 34 56 78 100 87 65 43 21 200 a4 b6 c888
R data.table：在当前测量之前对出现次数进行计数

我有一组在几天内进行的测量结果测量次数通常为 4 任何测量中可以捕获的数字范围为 1 5 在现实生活中给定测试集范围可能高达 100 或低至 20 我想每天计算每个值在当天之前发生的次数让我用一些示例数据来解释 test data
将整数“20160119”转换为“日”“年”“月”的不同列

如何将一列整数转换为日期 DATE PRCP 1 19490101 25 2 19490102 5 3 19490118 18 4 19490119 386 5 19490202 38 到这样的表 days month years PRCP
R 下与 Julia 的密切联系

In my R代码我使用MixedModels Julia包裹我整合Julia in R使用JuliaCall包裹我使用非常大的数据集 1 GB 4x10 6观察和建模步骤 mixed models 分配了大量 RAM 80 GB
R 如何按行值进行分组、拆分或子集

这是上一个问题的延续R 如何按行值分组分裂 https stackoverflow com questions 64602607 r how to group by row value split 输入数据帧的变化是 id str c x
使用 R 中“rpart”包中的生存树来预测新的观察结果

我正在尝试使用 R 中的 rpart 包来构建生存树并且我希望使用这棵树来对其他观察结果进行预测我知道有很多涉及 rpart 和预测的问题但是我还没有找到任何解决我认为特定于将 rpart 与 Surv 对象一起使用的问题的方法
在 R 中，如何让 PRNG 在平台之间给出相同的浮点数？

在 R 4 1 1 中运行以下代码会在平台之间产生不同的结果 set seed 1 x lt rnorm 3 3 print x 22 0 83562861241004716 intel windows 0 8356286124100471
使用示例代码继续在 ggplot2 中遇到错误“loop_apply”未从当前命名空间（plyr）解析”

我今天一直遇到这个错误我已经从 github 下载了 plyr 但它仍然不起作用安装 plyr 后我重新启动了 R studio 甚至我的电脑看来问题可能是由于 R 解析对外部 DLL 的引用的方式发生了变化正如线程中途提到的he
如何使用PIL将灰度图转为伪彩色？

我似乎不知道如何获取我的灰度函数并将其更改为给我假颜色我知道我需要将每种颜色 R G B 分成范围然后根据每种颜色的范围分配颜色有谁知道这是如何运作的 def grayscale pic width height pic size f
按元素名称组合/合并列表

我有两个列表其元素的名称部分重叠我需要将其逐个元素合并组合成一个列表 gt lst1 lt list integers c 1 7 letters letters 1 5 words c two strings gt lst2 lt
ggplot多个时间不等的时间序列

我知道有一些与时间序列和多个数据帧相关的已回答问题但我似乎无法弄清楚这一点我想绘制 4 个不同压力传感器与时间的时间戳数据 pa 列我有来自同一实验的 4 个带时间戳的压力读数 dfs 然而由于传感器故障和数据中的其他信号传感器收
将相同的 Patch 实例添加到 matplotlib 中的多个子图中

我正在尝试将补丁的相同实例添加到 matplotlib 中的多个轴这是最小的例子 import matplotlib pyplot as mpl plt import matplotlib patches as mpl patches f
如何更改 R Markdown HTML 文档中目录的颜色和属性？

我花了很多时间谷歌搜索这个但似乎无法弄清楚我正在使用 R Markdown 制作 HTML 文档文档在这里 http rmarkdown rstudio com html document format html http rmark
包含多行 LaTeX 方程的 R 帮助页面

我正在编写 R 包文档roxygen2 我想将以下多行方程插入到帮助页面中但我的 LaTeX 代码没有被渲染 hello2 description deqn F t begin cases alpha t f L t 1 alpha t

随机推荐

重置用户密码

我正在尝试找到一种通过非交互式登录在 Azure Active Directory 中重置用户密码所有用户而不仅仅是经过身份验证的用户的解决方案目前看来这只能通过 powershell 的 MSOnline 获得Set AzureA
Android 模拟器替代品

我对 Android 开发完全陌生但我刚刚拥有一台 HTC Hero 想为其开发一些应用程序然而我使用笔记本电脑作为我的开发机器并且模拟器非常慢启动大约需要 10 15 分钟虽然我可以让它保持打开状态但在使用其他应用程序如
LINQ 按空列排序，其中顺序为升序，空值应该在最后

我正在尝试按价格对产品列表进行排序结果集需要按列按价格从低到高列出产品LowestPrice 但是该列可以为空我可以按降序对列表进行排序如下所示 var products from p in context Products whe
如何从Application.Path获取UNC路径？

我想获取 vba 代码中活动工作簿的路径 ActiveWorkbook Path做这个 BUT 我需要它来检索这样的东西 MachineName ShareFolder ETC ETC2 NOT S ETC ETC2 Where S 映射到
为什么在访问模型时，backbone.js 返回一个空数组？

我有一个路由器访问其集合我的 for 循环没有迭代模型因此我尝试记录集合以查看它返回的内容事实证明当我直接记录集合时我会按预期看到所有模型但是如果我尝试记录集合的 models 属性我会得到一个空数组这没有道理这些线直接
C# 中 async/await 的这种用法以前被发现过吗？ [关闭]

很难说出这里问的是什么这个问题模棱两可含糊不清不完整过于宽泛或言辞激烈无法以目前的形式合理回答如需帮助澄清此问题以便重新打开访问帮助中心在我之前在 stackoverflow 上提出了关于 async await 的问题之后
以编程方式查明进程是否需要用户输入

我如何以编程方式在 C 中确定另一个外部应用程序本机 java NET 或其他当前是否需要用户输入这可以在托管代码中完全完成吗我正在寻找的是实施 static Boolean IsWaitingForUserInput Stri
Google Maps API 自动完成同一页面上的第二个地址字段

我在我的页面上使用 Google Maps API 该页面要求用户填写您的当前地址和新地址我可以让自动完成功能在第一个地址上工作但它不适用于第二个地址我做了很多研究并查看了 stackoverflow 上的类似帖子但我找不到任
在 try/catch 块中等待两个承诺会导致“未处理的承诺拒绝”[重复]

这个问题在这里已经有答案了我想等待两个并行运行的承诺我不想连续等待每个承诺这有效但速度较慢出于这个原因我认为我可以首先创建两个承诺来让它们滚动比如说两个网络请求然后等待它们并能够在 catch 块中捕获错误这个假设似乎是不正
Django - 按模板中的某个字段对查询集进行分组

我有一张桌子Events 按字段排序date 我想打印模板中的事件但为每个日期使用单独的 div 例如 div class content h1 December 30th h1 div div class content h1 Dece
#1025 - mysql 中重命名错误（errno：150）

我试图在一个表 misc 中删除一个外键 id 它是表 main 中的主键 id 数据库名称 xxx alter table misc drop FOREIGN KEY id 我收到这个错误 1025 将 interview sql edc
如何使用 dotnet CLI 构建 .NET Framework 4.8 应用程序？

我继承了一个 NET Framework 4 8 应用程序我可以使用 Visual Studio 对其进行编译但是当尝试使用dotnetCLI 应用程序对于我从外部依赖项使用的所有类型我都会收到 CS0246 错误原来的错误消息
Spring授权服务器：如何使用托管在单独应用程序上的登录表单？

我正在使用 Spring Security 和 Spring Authorization Server 并尝试创建身份验证服务器我有一个基本流程允许我使用预先构建的登录页面从拜尔东指南这是我正在处理的代码我假设这个登录页面表单来自
学习 LINQ [关闭]

Closed 此问题正在寻求书籍工具软件库等的推荐不满足堆栈溢出指南目前不接受答案 Overview 我在这个网站上问过很多次的事情之一是LINQ 我提出的问题广泛而多样而且往往没有太多背景因此为了巩固我在 Linq 上获得的
ASP.NET 传递 Windows 身份验证凭据

我有一个使用 Windows 身份验证的 ASP NET Web 应用程序此应用程序需要连接到另一个 ASP NET Web 服务也使用 Windows 身份验证并使用从用户浏览器收到的相同凭据这可以做到吗如何做到我不相信没有
C++ 错误：转换为执行字符集

问题是我试图打印 ASCII 之外的一些字符例如德语元音变音字符等这些字符不适合普通的 char 变量因此显然我尝试将它们放入 wchar t 中并使用 L 初始化字符串但是每次这个字符串包含上面的一个字符时我都会收到上述错误
从 CakePHP 2.1.2 控制台 Shell 创建完整 url

我正在尝试通过控制台 shell 从 CakePHP 2 1 2 发送电子邮件最终通过 cron 作业我发送的视图是一个日历其中包含返回应用程序网页的链接我发现的问题是网址不包含正确的路径从我读到的内容来看这是因为我使用控制台后
GIF/JPEG 文件可以包含可运行的 PHP 代码吗？

测试 Web 应用程序我可以上传 GIF JPEG 文件并且我知道在该 GIF Jpeg 中正确包含 PHP 代码时可能存在威胁因为它是使用图像创建真彩色 and 图像jpeg 我正在寻找包含简单 PHP 代码如 phpinfo 或
GroovyWS 和复杂请求

我遇到了使用 GroovyWS 发送复杂请求的问题这是由soapUI生成的示例请求
ggplot2 中两个不同组的不同调色板

我正在尝试创建一个图显示多个地点的观测数据和建模数据的每月土壤湿度垂直剖面到目前为止我只能绘制一组值无论是观察值还是建模值如下例所示 library ggplot2 library RColorBrewer Create cust

ggplot2 中两个不同组的不同调色板

ggplot2 中两个不同组的不同调色板 的相关文章

随机推荐

热门标签

ggplot2 中两个不同组的不同调色板的相关文章