私はまだプロットの正しいラベルに苦労しています。私はこれを機能させることができました。
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},")")))
欠落しているのは、mumol と day^-1 の間の演算子記号 " * " です。次のようになり
ます。ご協力ありがとうございます。