Проблемы с добавлением легенды в 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")

введите здесь описание изображения


r ggplot2 legend ggfortify
person Forever    schedule 19.11.2017    source источник

Ответы (2)


arrow_upward
1
arrow_downward

Решение с 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")

введите здесь описание изображения

person Marco Sandri    schedule 19.11.2017

arrow_upward
0
arrow_downward

Похоже, это ошибка с ggfortify.* Однако вы можете добиться идентичных результатов, просто используя geom_line() из ggplot2:

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"))

Изображение результатов

* Связанный вопрос: Нанесение нескольких распределений плотности на один график

person momo    schedule 08.11.2020