Добавление легенды в график множественной регрессии - PullRequest
Новая
       13

Добавление легенды в график множественной регрессии

0 голосов
/ 10 мая 2018

Я пытаюсь построить 4 регрессии, все с одним и тем же предиктором, с соответствующими точками в 4 разных цветах. Мне нужна легенда, которая назначает 4 метки для 4 цветов. Просто простая таблица с 2 столбцами. Однако при попытке установить руководство с помощью scale_linetype_manual или scale_fill_manual, легенда не появляется. 

Forest_EPs_pure <- data.frame("Sand"=c(23,41,32,34,38,49,32,18,91,18,117,61,68,43,46,59,74,88,58,92,882,941,870,926,861,848,810,964,874,860,755,942),
                              "B2_fs"=c(61,68,80,20,43,72,93,60,75,75,97,83,74,51,87.5,74,97,74,71,60,0,0,0,0,0,0,17,0,0,0,0,0),
                              "H_total"=c(89.7,109.6,15.4,140.1,145.4,164,11.6,118.1,110.8,50.4,77.8,124.6,115.6,152,63.6,112.9,127.7,138.5,69.8,149.4,32,63.4,35.7,84.5,5.8,0,8.8,1.6,20.2,31.8,25.8,0.5),
                              "B1_fs"=c(0,0,0.5,0,0.5,0,3,0,0,0.5,2,10,0,0,7,0,0,20,15,18,75,71,115,72,65,95,40,35,43,96,98,95),
                              "S_total"=c(8.6,8.5,2.2,18.9,37.9,53.1,3.8,76,67.4,4.8,35.2,78,74.5,65.6,41.6,58,37.1,60.7,68,39.5,44.4,69.6,5.5,42.1,19.3,0,77.1,0.5,96.2,2.0,8.5,0.5))



    ggplot(data=Forest_EPs_pure, aes(x=Sand))+
      geom_point(aes(y=B2_fs), color="green")+
      geom_point(aes(y=H_total), color="black")+
      geom_point(aes(y=B1_fs), color="brown")+
      geom_point(aes(y=S_total), color="grey")+
      geom_smooth(aes(y=B2_fs), method="glm",color="green")+
      geom_smooth(aes(y=H_total), method="glm", color="black")+
      geom_smooth(aes(y=B1_fs), method="glm", color="brown")+
      geom_smooth(aes(y=S_total), method="glm", color="grey")+
      scale_x_continuous("Sand content in g/kg")+
      scale_y_continuous("Share of respective layer in %")

1 Ответ

0 голосов
/ 10 мая 2018

ggplot обычно ожидает данные в длинном формате. Для этого вы можете использовать gather из tidyr.

Например, с этими данными у вас будет столбец, содержащий «B1_fs», «B2_fs», «H_total» и «S_total», или любые другие столбцы, которые могут быть у вас до изменения формы данных. Это позволяет вам сопоставить эту переменную с цветом, что дает вам легенду, и установить цвета вручную с помощью именованного вектора в scale_color_manual.

Как правило, вы не хотите использовать несколько geom одного и того же типа, чтобы делать что-то похожее, например, многократно использовать geom_point для построения точек в разных цветах. 

library(ggplot2)

Forest_EPs_pure <- data.frame("Sand"=c(23,41,32,34,38,49,32,18,91,18,117,61,68,43,46,59,74,88,58,92,882,941,870,926,861,848,810,964,874,860,755,942),
                                                            "B2_fs"=c(61,68,80,20,43,72,93,60,75,75,97,83,74,51,87.5,74,97,74,71,60,0,0,0,0,0,0,17,0,0,0,0,0),
                                                            "H_total"=c(89.7,109.6,15.4,140.1,145.4,164,11.6,118.1,110.8,50.4,77.8,124.6,115.6,152,63.6,112.9,127.7,138.5,69.8,149.4,32,63.4,35.7,84.5,5.8,0,8.8,1.6,20.2,31.8,25.8,0.5),
                                                            "B1_fs"=c(0,0,0.5,0,0.5,0,3,0,0,0.5,2,10,0,0,7,0,0,20,15,18,75,71,115,72,65,95,40,35,43,96,98,95),
                                                            "S_total"=c(8.6,8.5,2.2,18.9,37.9,53.1,3.8,76,67.4,4.8,35.2,78,74.5,65.6,41.6,58,37.1,60.7,68,39.5,44.4,69.6,5.5,42.1,19.3,0,77.1,0.5,96.2,2.0,8.5,0.5))

forest_long <- tidyr::gather(Forest_EPs_pure, key = variable, value = value, -Sand)
head(forest_long)
#>   Sand variable value
#> 1   23    B2_fs    61
#> 2   41    B2_fs    68
#> 3   32    B2_fs    80
#> 4   34    B2_fs    20
#> 5   38    B2_fs    43
#> 6   49    B2_fs    72

ggplot(forest_long, aes(x = Sand, y = value, color = variable)) +
    geom_point() +
    geom_smooth(method = glm) +
    scale_color_manual(values = c(B2_fs = "green", H_total = "black", B1_fs = "brown", S_total = "grey"))

Создано в 2018-05-10 пакетом Представ (v0.2.0).

...