将图例添加到 ggplot2 + ggfortify 时出现问题

Problems adding legend to ggplot2 + ggfortify

我在使用

时遇到问题
scale_colour_manual 

ggplot 的函数。我试过了

guide = "legend" 

强制图例出现,但不起作用。代表代码:

library(ggfortify)
library(ggplot2)
p  <- ggdistribution(pgamma, seq(0, 100, 0.1), shape = 0.92, scale = 22, 
                     colour = 'red')
p2 <- ggdistribution(pgamma, seq(0, 100, 0.1), shape = 0.9, scale = 5, 
                     colour = 'blue', p=p)

p2 + 
theme_bw(base_size = 14) +
theme(legend.position ="top") +
xlab("Precipitación") +
ylab("F(x)") +
scale_colour_manual("Legend title", guide = "legend", 
                      values = c("red", "blue"), labels = c("Observado","Reforecast")) +
ggtitle("Ajuste Gamma")

stat_function的解决方案:

library(ggplot2)
library(scales)

cols <- c("LINE1"="red","LINE2"="blue")
df <- data.frame(x=seq(0, 100, 0.1))
ggplot(data=df, aes(x=x)) + 
stat_function(aes(colour = "LINE1"), fun=pgamma, args=list(shape = 0.92, scale = 22)) +
stat_function(aes(colour = "LINE2"), fun=pgamma, args=list(shape = 0.9, scale = 5)) +
theme_bw(base_size = 14) +
theme(legend.position ="top") +
xlab("Precipitación") +
ylab("F(x)") +
scale_colour_manual("Legend title", values=c(LINE1="red",LINE2="blue"),
                    labels = c("Observado","Reforecast")) +
scale_y_continuous(labels=percent) +
ggtitle("Ajuste Gamma")

这似乎是 ggfortify 的错误。*您可以简单地使用 ggplot2 中的 geom_line() 获得相同的结果:

library(ggplot2)

# Sequence of values to draw from dist(s) for plotting
x = seq(0, 100, 0.1)

# Defining dists
d1 = pgamma(x, shape=0.92, scale=22)
d2 = pgamma(x, shape=0.90, scale=5)

# Plotting
p1 = ggplot() +
     geom_line(aes(x,d1,colour='red')) + 
     geom_line(aes(x,d2,colour='blue')) + 
     
     theme_bw(base_size = 14) +
     theme(legend.position="top") +

     ggtitle("Ajuste Gamma") +
     xlab("Precipitación") +
     ylab("F(x)") +
     
     scale_colour_manual("Legend title", 
                         guide = "legend", 
                         values = c("red", "blue"), 
                         labels=c("Observado", "Reforecast"))


* 相关问题:Plotting multiple density distributions on one plot