如何用两个 y 轴绘制这些数据？

Question

我有数据：

structure(list(Group.1 = structure(c(17897, 17928, 17956, 17987, 
18017, 18048, 18078, 18109, 18140, 18170, 18201, 18231, 18262, 
18293, 18322, 18353, 18383, 18414, 18444, 18475, 18506, 18536, 
18567, 18597), class = "Date"), Total.y = c(431L, 399L, 424L, 
421L, 390L, 383L, 397L, 403L, 476L, 507L, 505L, 556L, 604L, 618L, 
39491L, 119263L, 111533L, 98341L, 77406L, 64553L, 116157L, 130517L, 
130685L, 149184L), t_VidP = c(92.3433874709977, 94.7368421052632, 
91.2735849056604, 94.061757719715, 96.6666666666667, 93.9947780678851, 
94.7103274559194, 95.0372208436725, 94.9579831932773, 91.3214990138067, 
93.8613861386139, 94.9640287769784, 93.5430463576159, 93.3656957928803, 
73.8699956952217, 79.24251444287, 81.4395739377583, 81.8275185324534, 
83.9921969873136, 85.6954750360169, 92.2596141429272, 93.4452983136296, 
93.5248880896813, 92.8739006864007)), row.names = c(NA, -24L), class = "data.frame")

看起来像：

Group.1             Total.y   t_VidP
    1 2019-01-01     431 92.34339
    2 2019-02-01     399 94.73684
    3 2019-03-01     424 91.27358
    4 2019-04-01     421 94.06176
    5 2019-05-01     390 96.66667
    6 2019-06-01     383 93.99478

我想要一个左侧 y 轴上带有 Total.y 的图表。 以及右侧 y 轴上的 t_VidP（范围从 0 到 100）。 我希望 x 轴是月份。 我还想要一个 x 轴标题和 y 轴标题。 理想情况下，它看起来像： https://imgur.com/a/VjuLm0O

Answer 1

这是两个带有第二个 y 轴的图，其中第二个带有 log10 比例的 y 轴，因为Total.y范围很大。 另外，请参阅此 SO 帖子。

至于第二个 y 轴，诀窍是计算比例因子M 。 由于百分比列在 0-100 范围内，将其除以 100 并使用scales::percent自动格式化百分比。

library(ggplot2)
library(scales)

M <- max(df1$Total.y)

ggplot(df1, aes(Group.1)) +
  geom_line(aes(y = Total.y)) +
  geom_line(aes(y = M * t_VidP/100), linetype = "dashed") +
  scale_x_date(date_labels = "%Y %b") +
  scale_y_continuous(
    sec.axis = sec_axis(~ . / M, name = "% Vid", labels = percent)
  ) +
  ggtitle("y axis in natural scale") +
  labs(x = "Months", y = "Total") +
  theme_bw()

ggplot(df1, aes(Group.1)) +
  geom_line(aes(y = Total.y)) +
  geom_line(aes(y = M * t_VidP/100), linetype = "dashed") +already 
  scale_x_date(date_labels = "%Y %b") +
  scale_y_continuous(
    trans = "log10",
    breaks = trans_breaks("log10", function(x) 10^x),
    labels = trans_format("log10", math_format(10^.x)),
    sec.axis = sec_axis(~ . / M, name = "% Vid", labels = percent)
  ) +
  ggtitle("y axis in log10 scale") +
  labs(x = "Months", y = "Total") +
  theme_bw()

Answer 2

此外，基本 R 方法将是：

par(mar = c(5, 4, 4, 4) + 0.3)
plot(df$Group.1, df$Total.y, col = "red", type = "l", xlab = "Date", ylab = "Total")
par(new = T)
plot(df$Group.1, df$t_VidP, col = "blue", type = "l", xlab = "", ylab = "", axes = F)
axis(side = 4, at = pretty(range(df$t_VidP)))
mtext("%Vid", side = 4, line = 3)
legend(x = "left", legend=c("Total", "%Vid"),
       col=c("red", "blue"), lty=1, cex=0.8)

Answer 3

我试图让 aes 远离 y。

library(tidyverse)
library(lubridate)
df %>% 
  mutate(Total.y_transformed = Total.y/(43485.2 /90.54 )) %>% 
  mutate(month = month(ymd(Group.1), label = TRUE)) %>% 
  pivot_longer(
    cols = c(Total.y_transformed, t_VidP),
    names_to = 'names',
    values_to = 'value'
  ) %>% 
  arrange(names)

ggplot(df1, aes(x = Group.1, y = value, color=names)) + 
  geom_line(size=1) +   
  scale_color_manual(labels = c("t_VidP", "Total"), values = c("blue", "red")) + 
  scale_y_continuous(sec.axis = sec_axis(~.*(43485.2 /90.54), name = "Total")) +
  scale_x_date(date_labels = "%Y %b") +
  ggtitle("Title") +
  labs(x = "Months", y = "t_VidP",  colour = "Parameter") +
  theme(legend.position = c(0.8, 0.9)) +
  theme_classic(base_size=12)