如何在R>中使用双Y轴将图例添加到ggplot

debugcn 发表于 Dev

王爱伦

我正在尝试找到一种将图例添加到ggplot图中的方法。我尝试了多种方法，但是没有一种方法可以让图例出现。这是我正在使用的数据集。以下是我到目前为止编写的代码。

library(dplyr)
library(tidyr)
library(ggplot2)
library(lubridate)
library(ggthemes)

stock %>%
  drop_na(disease, US_Stock, China_Stock) %>%
  filter(disease == "Ebola") %>%
  ggplot(aes(x = date)) + 
    geom_line(aes(y = US_Stock), size = 0.5, col = "dark green") +
    geom_line(aes(y = China_Stock/5+1500), size = 0.5, col = "red") +
    scale_x_date(name = "Date") +
    scale_y_continuous(name = "US Stock Market (S&P 500)", 
                       sec.axis = sec_axis(~(.-1500)*5, name = "China Stock Market (CSI300)")) +
    labs(title = "Figure 3: Stock Market under Ebola", caption = "Death Rate: 50%-90%")  +
    theme_stata()

这是我数据中的一些观察结果。

爱德华

我更喜欢将所有股票重塑为一列，因为这有助于传说和其他美学。我还使用了平均值来确定比率，这使得图形更加有趣。（但我不喜欢双轴）。

library(dplyr)
library(tidyr)
library(ggplot2)
library(lubridate)
library(ggthemes)
ratio = mean(stock$China_Stock)/mean(stock$US_Stock)

stock %>%
  mutate(China_Stock=China_Stock / ratio) %>%
  pivot_longer(cols=ends_with("Stock"), names_pattern="(.+)_(Stock)",
               names_to = c("Country", ".value")) %>%
  ggplot(aes(x = date, y=Stock, col=Country)) + 
  geom_line(size = 0.5) +
  scale_color_manual(values=c("darkgreen","red")) +
  labs(title = "Figure 3: Stock Market under Ebola", caption = "Death Rate: 50%-90%") +
  scale_y_continuous(name = "US Stock Market (S&P 500)", 
                       sec.axis = sec_axis(~.*ratio, name = "China Stock Market (CSI300)")) +
  theme_minimal()

数据：

stock <- structure(list(ID = 131:140, date = structure(c(11878, 11879, 
                                                        11880, 11883, 11884, 11885, 11886, 11887, 11890, 11891), class = "Date"), 
                       US_Stock = c(920.47, 927.37, 921.39, 917.93, 901.05, 906.04, 
                                    881.56, 847.76, 819.85, 797.7), China_Stock = c(1392.9, 1390.93, 
                                                                                    1389.45, 1379.38, 1382.41, 1393.02, 1397.41, 1403.25, 1380.28, 
                                                                                    1375.61), disease = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 
                                                                                                                    1L, 1L, 1L), .Label = "SARS", class = "factor")), class = "data.frame", row.names = c(NA, 
                                                                                                                                                                                                          -10L))

本文收集自互联网，转载请注明来源。

如有侵权，请联系[email protected] 删除。