CODEWITHSUNDEEP
rggplot2plottidyverseyaxis
Hydro
Hydro

Reputation: 1117

Secondary y-axis in ggplot in R?

Went through a few of secondary axis solutions proposed here but didn't get it right. I am trying to plot Elevation on the left y-axis and FlowA & Flowb on the right y-axis. My sample code will do the Elevation plotting however, struggling to get the FlowA & FlowB variables on the secondary axis. Any help would be appreciated.

library(lubridate)
library(tidyverse)

set.seed(123)

FakeData <- data.frame(Date = seq(as.Date("2001-01-01"), to= as.Date("2001-12-31"), by="day"),
                 Elevation = runif(365, 806.8,807.8),
                 FlowA = runif(365,8,15),
                 FlowB = runif(365,1,3))
ggplot(FakeData, aes(x = Date, y = Elevation))+
  geom_line()

Upvotes: 0

Views: 523

Answers (2)

Duck
Duck

Reputation: 39613

I would suggest next approach. Also to mention that the output would depend on your data. Here all Elevation values are close to 800. For second axis, you have to define a scaling factor around all variables so that they are properly showed. Next the code:

library(lubridate)
library(tidyverse)

set.seed(123)
#Data
FakeData <- data.frame(Date = seq(as.Date("2001-01-01"), to= as.Date("2001-12-31"), by="day"),
                       Elevation = runif(365, 806.8,807.8),
                       FlowA = runif(365,8,15),
                       FlowB = runif(365,1,3))
#Scale factor
scalefactor <- max(FakeData$Elevation)/max(max(FakeData$FlowA),max(FakeData$FlowB))
#Plot
ggplot(FakeData, aes(x = Date))+
  geom_line(aes(y = Elevation,group=1,color='Elevation'),show.legend = T)+
  geom_line(aes(y = FlowA*scalefactor, color = 'FlowA'))+
  geom_line(aes(y = FlowB*scalefactor, color = 'FlowB'))+
  scale_y_continuous(sec.axis = sec_axis(~./scalefactor, name = 'Flow A and Flow B'))

The output:

enter image description here

Upvotes: 1

Matthew Farant
Matthew Farant

Reputation: 82

The solution in here works just fine for me. But I would like to add some adjustments because the code will be slightly different for your data. According to the solution, we will use 3 geom objects that represent the elevation, FlowA, and FlowB. We will also make the secondary axis for FlowA and FlowB.

ggplot(FakeData, aes(x = Date))+
  geom_line(aes(y = Elevation)) +
  geom_col(aes(y = FlowA), fill="blue") +
  geom_col(aes(y = FlowB), fill='red')+
  scale_y_log10(sec.axis = sec_axis(~ .*1, labels = scales::number_format(scale=1/10),name="Flow"))

In the code above, I will show the elevation as a line plot and flows as a bar plot. Why did I use the logarithmic scale here? Because the distance of the elevation's value range(between 806.8 until 807.8) and the flows' value range is very far. If you proceed with a regular y-axis (scale_y_continuous()), you will have this plot below: enter image description here See that the plot is not so meaningful. You can't see clearly how the flows change over time. Here's what it looks in logarithmic scale: enter image description here I use the logarithmic scale for the left y-axis and regular scale on the right y-axis. Now we can see clearly the changes in the flows over time. The Elevation will definitely be a straight line because you set it to be a random uniform distribution.

However, personally, I don't suggest you use a double y-axis because it confuses the plot user. I suggest you split the plot into two different plots.

Upvotes: 0

Related Questions