I am still struggling with correct labels of my plots. I managed to get this to work.
df <-data.frame(a=seq(1,20),b=seq(1,20))
ggplot(df, aes(x = a, y = b)) + geom_point() +
ylab(expression(paste(CH[4]~" formation rate"~"(",mu, mol, day^{-1},")")))
What it is missing is the operator sign " * " between mumol and day^-1. It should look like this:
Thank you for your help.