Соответствие значениям двойной оси Y в R

Я не уверен, как вы делаете это в базе R, но вы можете выровнять вторую ось с помощью небольшой математики и проб и ошибок с помощью ggplot2

library(tidyverse)

dat |>
  mutate(SPEI = (SPEI+(25/15))*15)|>
  pivot_longer(VCI:SPEI) |>
  mutate(name = factor(name, levels = c("VCI", "TCI", "SPEI"))) |>
  ggplot(aes(Time, value, color = name))+
  geom_line()+
  geom_hline(linetype = "dotted", yintercept = 40)+ 
  geom_hline(linetype = "dashed", yintercept = 2.5)+
  scale_color_manual(values = c("Aquamarine3", "Chocolate2", "darkorchid3"))+
  scale_y_continuous(name = "VCI and TCI", 
                     breaks = seq(0, 100, by = 20),
                     limits = c(0, 100),
                     sec.axis = sec_axis(trans = ~(./15)-(25/15), name = "SPEI", 
                                         breaks = seq(-2,5, by = 1)))+
  ggtitle("Temporal trend - drought indices, growing season")+
  labs(color = "")+
  theme(panel.background = element_blank(),
        panel.border = element_rect(fill = NA, color = "black"),
        legend.key = element_blank(),
        legend.box.background = element_rect(fill = NA, color = "black"),
        legend.position = c(0.07,0.87),
        legend.title = element_blank())

ОБНОВЛЯТЬ отредактировано, чтобы сделать 40 = -1

library(tidyverse)

dat |>
  mutate(SPEI = ((SPEI+(55/15))*15))|>
  pivot_longer(VCI:SPEI) |>
  mutate(name = factor(name, levels = c("VCI", "TCI", "SPEI"))) |>
  ggplot(aes(Time, value, color = name))+
  geom_line()+
  geom_hline(linetype = "dotted", yintercept = 40)+ 
  geom_hline(linetype = "dashed", yintercept = 32.5)+
  scale_color_manual(values = c("Aquamarine3", "Chocolate2", "darkorchid3"))+
  scale_y_continuous(name = "VCI and TCI", 
                     breaks = seq(0, 100, by = 20),
                     limits = c(0, 100),
                     sec.axis = sec_axis(trans = ~(./15)-(55/15), name = "SPEI", 
                                         breaks = seq(-4,3, by = 1)))+
  ggtitle("Temporal trend - drought indices, growing season")+
  labs(color = "")+
  theme(panel.background = element_blank(),
        panel.border = element_rect(fill = NA, color = "black"),
        legend.key = element_blank(),
        legend.box.background = element_rect(fill = NA, color = "black"),
        legend.position = c(0.07,0.87),
        legend.title = element_blank())

Вопрос можно решить, вычислив максимальное значение одной из осей от другой по формуле линейной интерполяции, объяснение смотрите в Википедии.

Здесь я предпочитаю дать более компактный и общий ответ, который можно повторно использовать в подобных случаях. В приведенном ниже решении используются случайные данные и не используется маркировка для конкретной задачи.

Это работает следующим образом:

1. Создайте несколько случайных данных с одной серией для x и двумя сериями ya и yb для y.
2. Установите пределы оси Y для ya, минимум для yb и соответствующие значения для обеих осей, но не указывайте максимальное значение для второй оси.
3. Рассчитайте максимум для yb.
4. Дальше можно делать все как обычно.

## create a few random test data
set.seed(123)
x <- 1:10
ya <- runif (x, min = 20, max = 70)
yb <- runif (x, min = -1.5, max= 1)

## set limits for the two y axes, but omit yb_max
ya_min <- 0; ya_max <- 100; yb_min <- -3

## set value, where the two axes should match
ya_match <- 40; yb_match <- -1

## calculate yb_max using linear interpolation (here: extrapolation)
yb_max <- yb_min + (ya_max - ya_min) * (yb_match - yb_min) / (ya_match - ya_min)

## set graphical parameter to increase margin for right axis
par(mar=c(5, 4, 5, 5) + .1)

## plot data set "a" 
plot(x, ya, type = "l", col = "black", ylim = c(ya_min, ya_max))
abline(h = ya_match, col = "grey", lty = "dashed")

## add a second plot for data set "b", omitting axes and axis labels
par(new = TRUE)
plot(x, yb, type = "l", col = "red", ylim = c(yb_min, yb_max),
     axes = FALSE, xlab = "", ylab = "")

## add right axis and axis label
axis(4)
mtext("yb", side = 4, line = 3, col = "red")