R ggplot2 - 在所有方面的顶部和右侧添加刻度

Question

我用来ggplot2创建以下情节：

library(ggplot2)

df <- data.frame(site=rep(c("CA-Oas","US-Oho","CA-Obs","CA-TP1"), each=124),
                 exper=rep(c("bberry","medlyn"), each=62, 4),
                 sm=runif(496, min=0, max=100),
                 lh=runif(496, min=0, max=400),
                 month=rep(c("Jul","Aug"), each=31, 8))

ggplot(df, aes(x=sm, y=lh, colour=exper, shape=month)) + 
  geom_point(size=2.5, alpha=0.9) +
  facet_wrap(~site, nrow=2) +
  theme_bw(base_size=18) +
  ylab(bquote('Latent heat ('*W~m^-2*')')) +
  xlab(bquote('Soil water content (%)')) +
  scale_color_manual(values=c("bberry"="tomato2", "medlyn"="dodgerblue"),
                     labels=c("bberry"="Ball-Berry", "medlyn"="Medlyn")) +
  theme(panel.grid.minor=element_blank(),
        panel.grid.major=element_blank(),
        legend.title=element_blank(),
        legend.box="horizontal",
        legend.position=c(0.4, 0.6),
        legend.spacing = unit(-0.2,"cm"))

然而，我正在写手稿的期刊要求所有图在两个轴（上/下和左/右）上都有刻度。我正在努力实现这一目标ggplot2。

基本上，我需要的是：

1. 在所有面上的顶部和右侧轴上添加刻度（内部或外部）；
2. 仅将 x 轴标签保留在面的底部，将 y 轴标签保留在面的左侧。

这大致就是我需要实现的目标（注意轴刻度和标签）：

关于如何到达那里有任何提示吗？

Answer 1

facet_rep_wrap柠檬包中的非常适合此目的，因为它在所有方面都保留了轴线（以及可选的标签）。

有关代码的其他调整，请参阅下面的内联注释：

library(lemon)

ggplot(df, aes(x=sm, y=lh, colour=exper, shape=month)) + 
  geom_point(size=2.5, alpha=0.9) +
  facet_rep_wrap(~site, nrow = 2) +              # instead of facet_wrap
  theme_bw(base_size=18) +
  ylab(bquote('Latent heat ('*W~m^-2*')')) +
  xlab(bquote('Soil water content (%)'))  +
  scale_x_continuous(sec.axis = dup_axis()) +    # add duplicated secondary axis to get 
  scale_y_continuous(sec.axis = dup_axis()) +    # axes on top / right sides
  scale_color_manual(values=c("bberry"="tomato2", "medlyn"="dodgerblue"),
                     labels=c("bberry"="Ball-Berry", "medlyn"="Medlyn")) +
  theme(panel.grid.minor=element_blank(),
        panel.grid.major=element_blank(),
        legend.title=element_blank(),
        legend.box="horizontal",
        legend.position=c(0.4, 0.6),
        legend.spacing = unit(-0.2,"cm"),
        strip.placement = "outside",             # place facet strips outside axis
        axis.ticks.length = unit(-2.75, "pt"),   # point axis ticks inwards (2.75pt is 
                                                 # the default axis tick length here)
        axis.text.x.top = element_blank(),       # do not show top / right axis labels
        axis.text.y.right = element_blank(),     # for secondary axis
        axis.title.x.top = element_blank(),      # as above, don't show axis titles for
        axis.title.y.right = element_blank())    # secondary axis either

（注意：仅使用 ggplot2 和一些 grob hacking 就可以达到相同的结果，但这种方法相当复杂，所以我不会去那里，除非在计算环境中安装新软件包确实很麻烦......）