ggplot2双纵坐标轴

library("scales")
library(tidyverse)
df4 <- tibble(
  year = c("2013", "2014", "2015", "2016", "2017", "2018", "2019", "2020"),
  ResearchInput = c(1.56, 2.54, 2.99, 5.63, 11.08, 12.31, 16.77, 25.92),
  ProportioninOperatingIncome = c(6.8, 6.9, 5, 4.9, 6.77, 5.6, 5.1, 4.75)
)

df4$year <- as.factor(df4$year)
df4$ProportioninOperatingIncome1 <- rescale(df4$ProportioninOperatingIncome, c(4.75, 28))
gg5 <- ggplot() +
  geom_col(data = df4, mapping = aes(x = year, y = ResearchInput), fill = '#1ec0ff', width = 0.6) +
  geom_line(data = df4, mapping= aes(x = year, y = ProportioninOperatingIncome1, group = 1), size = 1.2, colour = "#6441A4", linetype="longdash") +
  geom_point(df4, mapping = aes(x = year, y = ProportioninOperatingIncome1), size = 2.5, colour = "#4F86C6") +
  scale_y_continuous(
    breaks = seq(0, 28, 4),
    sec.axis = sec_axis(~rescale(.,c(4.5,7.5)),
                        labels = sprintf("%g%%",seq(4.5, 7.5, 0.5)))
  )

gg6 <- gg5 +
  geom_text(data = df4, mapping = aes(x = year,y = ResearchInput + 0.6,label = ResearchInput), color = "#1ec0ff") +
  geom_text(data = df4, mapping = aes(x = year,y = ProportioninOperatingIncome1 + 1.5,label = sprintf("%g%%", ProportioninOperatingIncome)), colour = "#4F86C6")

gg7 <- gg6 + theme(axis.ticks = element_line(linetype = "blank"),
    axis.text = element_text(size = 12),
    panel.background = element_rect(fill = NA)) +labs(x = NULL, y = NULL) +
  geom_hline(yintercept = 0, colour = "#DCDDD8")