Построение вторичной оси с использованием ggplot - PullRequest
Новая
       24

Построение вторичной оси с использованием ggplot

1 голос
/ 28 января 2020

Я пытаюсь построить три переменные (SA, SA1, SA2) с двумя переменными (SA & SA2) на левой оси y и одной переменной (SA1) на правой вторичной оси y. Я пытался зафиксировать пределы оси, используя limits = c(1e15,5e15) на левой оси Y, пытаясь ограничить вторичную ось между limits = c(3e17,4.2e17), но я не могу построить ось сёкондарного с моими настроенными пределами. DATA Link 

library(ggplot2)
test <- read.xlsx2("filepath/test.xlsx", 1, header=TRUE)
View(test)
test$SA=as.numeric(levels(test$SA))[test$SA]
test$SA1=as.numeric(levels(test$SA1))[test$SA1]
test$SA2=as.numeric(levels(test$SA2))[test$SA2]
g <- ggplot(test,aes(x=year, y=  SA, group = 1)) + geom_line(mapping = aes(x = test$year, y = test$SA)) 
+ geom_line(mapping = aes(x = test$year, y = test$SA2), color = "red") +  geom_line(mapping = aes(x = test$year, y = test$SA1), size = 1, color = "blue")
 g+scale_y_continuous(name = "primary axis title",
+                      sec.axis = sec_axis(~./5, name = "secondary axis title (SA1)"))

enter image description here enter image description here

Окончательное решение @ dc37 дает мне следующее результат: 

ggplot(subset(DF, Var != "SA1"), aes(x = year, y = val, color = Var))+
  geom_line()+
  scale_y_continuous(name = "Primary axis", sec.axis = sec_axis(~.*100, name = "Secondary"))

Спасибо enter image description here

1 Ответ

3 голосов
/ 29 января 2020

Аргумент sec.axis только создает новую ось, но не меняет ваши данные и не может использоваться для построения данных.

Чтобы иметь возможность наносить данные из двух групп с большим диапазоном, сначала необходимо уменьшить SA1.

Здесь я уменьшил его, разделив его на 100 (потому что отношение между максимальным значением SA1 и максимальным значением SA и SA2 близко к 100), и я также изменил ваш фрейм данных в более длинный формат, более подходящий для ggplot2: 

library(lubridate)
df$year = parse_date_time(df$year, orders = "%Y") # To set year in a date format
library(dplyr)
library(tidyr)
DF <- df %>% mutate(SA1_100 = SA1/100) %>% pivot_longer(.,-year, names_to = "Var",values_to = "val")

# A tibble: 44 x 3
    year Var         val
   <int> <chr>     <dbl>
 1  2008 SA      1.41e15
 2  2008 SA1     3.63e17
 3  2008 SA2     4.07e15
 4  2008 SA1_100 3.63e15
 5  2009 SA      1.53e15
 6  2009 SA1     3.77e17
 7  2009 SA2     4.05e15
 8  2009 SA1_100 3.77e15
 9  2010 SA      1.52e15
10  2010 SA1     3.56e17
# … with 34 more rows

Затем вы можете построить его, используя (я поднабор данных, чтобы удалить «SA1» и сохранить преобразованный столбец «SA1_100»): 

library(ggplot2)
ggplot(subset(DF, Var != "SA1"), aes(x = year, y = val, color = Var))+
  geom_line()+
  scale_y_continuous(name = "Primary axis", sec.axis = sec_axis(~.*100, name = "Secondary"))

enter image description here

Кстати, в ggplot2 вам не нужно проектировать столбец, используя $, просто напишите его имя.

Данные 

structure(list(year = 2008:2018, SA = c(1.40916e+15, 1.5336e+15, 
1.52473e+15, 1.58394e+15, 1.59702e+15, 1.54936e+15, 1.6077e+15, 
1.59211e+15, 1.73533e+15, 1.7616e+15, 1.67771e+15), SA1 = c(3.63e+17, 
3.77e+17, 3.56e+17, 3.68e+17, 3.68e+17, 3.6e+17, 3.6e+17, 3.68e+17, 
3.55e+17, 3.58e+17, 3.43e+17), SA2 = c(4.07e+15, 4.05e+15, 3.94e+15, 
3.95e+15, 3.59e+15, 3.53e+15, 3.43e+15, 3.2e+15, 3.95e+15, 3.03e+15, 
3.16e+15)), row.names = c(NA, -11L), class = c("data.table", 
"data.frame"), .internal.selfref = <pointer: 0x56412c341350>)
...