Проблемы с отображением легенды с использованием geom_pointrange в ggplot2

Question

Я использую ggplot2 с geom_pointrange в R для построения чувствительности (95% ДИ) пяти групп по шести сценариям (всего 30 графиков чувствительности).Шесть сценариев получены из набора сценариев 3x2.Я могу построить график успешно, но у меня проблемы с легендой.Я использовал три разных цвета и сплошные и пунктирные линии для представления шести сценариев.Однако я не могу отобразить пунктирные линии на легенде (только три цвета, каждый из которых повторяется дважды).Изначально я использовал шесть разных цветов, но коллеги попросили меня изменить это.

Я использовал geom_pointrange в ggplot2, но у меня были проблемы с отображением легенды:

sens <- cbind.data.frame(
  group = rep(c("Grp1", "Grp2", "Grp3", "Grp4", "Grp5"), 6),
  mean  = c(0.70, 0.70, 0.65, 0.60, 0.72, 0.85, 0.84, 0.77, 0.68, 0.77,
            0.71, 0.71, 0.69, 0.65, 0.76, 0.90, 0.88, 0.82, 0.82, 0.87,
            0.78, 0.78, 0.79, 0.73, 0.83, 0.92, 0.92, 0.92, 0.90, 0.93) * 100, 
  lower = c(0.64, 0.64, 0.59, 0.55, 0.68, 0.73, 0.73, 0.65, 0.55, 0.68,
            0.66, 0.66, 0.63, 0.59, 0.72, 0.80, 0.78, 0.70, 0.70, 0.79,
            0.73, 0.72, 0.74, 0.68, 0.80, 0.84, 0.84, 0.82, 0.79, 0.87) * 100,
  upper = c(0.75, 0.75, 0.70, 0.66, 0.76, 0.92, 0.91, 0.87, 0.80, 0.84,
            0.77, 0.76, 0.74, 0.70, 0.79, 0.95, 0.94, 0.90, 0.90, 0.92,
            0.82, 0.82, 0.83, 0.79, 0.86, 0.97, 0.97, 0.97, 0.96, 0.97) * 100,
  type = rep(c("A1&B1", "A1&B2", "A2&B1", "A2&B2", "A3&B1", "A3&B2"), each=5),
  type2 = rep(c(rep("B1", 5), rep("B2", 5)), 3),
  type3 = rep(c(rep("A1", 10), rep("A2", 10), rep("A3", 10))))


ggplot(sens, aes(group, mean, colour = type, group = type, linetype = type)) +
  geom_pointrange(aes(ymin = lower, ymax = upper), position = position_dodge(width = 0.5), size=1.25) + 
  coord_flip() +  
  xlab("Group") + 
  ylab("Sensitivity (95% CI)") +
  scale_color_manual(values = c("red", "red", "dark grey", "dark grey", "blue", "blue")) + 
  scale_linetype_manual(values=c(1,2,1,2,1,2)) +
  theme(axis.text.x=element_text(size=15, vjust=0.5, color = 'black'),  # x-axis labels
        axis.text.y=element_text(size=15, vjust=0.5, color = 'black'),  # y-axis labels
        axis.title.x=element_text(size=17.5, vjust=0.1),                # x-title justification  
        axis.title.y=element_text(size=17.5, vjust=1.5),                 # y-title justification
        legend.title=element_blank(),
        legend.position= "bottom",       
        legend.text = element_text(size = 14.5) 
  )

Я хочу видеть легенду с шестью различными надписями (3 разных цвета и 2 разных типа линий).

Answer 1

Вы можете удалить точки в легенде, чтобы сделать линии более заметными и увеличить размер легенды. Ключ:

ggplot(sens, aes(group, mean, colour = type, group = type, linetype = type)) +
  geom_pointrange(aes(ymin = lower, ymax = upper), position = position_dodge(width = 0.5), size=1.25) + 
  coord_flip() +  
  xlab("Group") + 
  ylab("Sensitivity (95% CI)") +
  scale_color_manual(values = c("red", "red", "dark grey", "dark grey", "blue", "blue"), 
                     ########### CHANGE HERE ###############
                     guide = guide_legend(override.aes = list(shape = NA))) + 
  scale_linetype_manual(values=c(1,2,1,2,1,2)) +
  theme(axis.text.x=element_text(size=15, vjust=0.5, color = 'black'),  # x-axis labels
        axis.text.y=element_text(size=15, vjust=0.5, color = 'black'),  # y-axis labels
        axis.title.x=element_text(size=17.5, vjust=0.1),                # x-title justification  
        axis.title.y=element_text(size=17.5, vjust=1.5),                 # y-title justification
        legend.title=element_blank(),
        legend.key.size = unit(2.5, "lines"), ########### CHANGE HERE ###############
        legend.position= "bottom",       
        legend.text = element_text(size = 14.5) 
  )

Другой подход может заключаться в изменении формы точек:

ggplot(sens, aes(group, mean, colour = type, group = type, shape = type)) +
  geom_pointrange(aes(ymin = lower, ymax = upper), position = position_dodge(width = 0.5), size=1.25) + 
  coord_flip() +  
  xlab("Group") + 
  ylab("Sensitivity (95% CI)") +
  scale_color_manual(values = c("red", "red", "dark grey", "dark grey", "blue", "blue")) + 
  scale_shape_manual(values = c(16,17,16,17, 16,17)) +
  theme(axis.text.x=element_text(size=15, vjust=0.5, color = 'black'),  # x-axis labels
        axis.text.y=element_text(size=15, vjust=0.5, color = 'black'),  # y-axis labels
        axis.title.x=element_text(size=17.5, vjust=0.1),                # x-title justification  
        axis.title.y=element_text(size=17.5, vjust=1.5),                 # y-title justification
        legend.title=element_blank(),
        legend.position= "bottom",       
        legend.text = element_text(size = 14.5) 
  )

Это также можно комбинировать с различными типами линий, ноэто может запутаться ...