在 R 中为这两个直方图添加图例

Question

我发现 R 和 ggplot 非常不直观。因此，任何指导将不胜感激。我已经尝试将我想要制作直方图的两个向量放入一个数据框中，但我永远无法让每个直方图的 alpha 值不同。

所以现在我使用代码：

ggplot(mydata) + geom_histogram(aes(x=mydata$D2prediction75At,y = ..density..,color = "blue"),color="darkblue", fill="blue",bins=41, alpha=0.75) 
 + geom_histogram(aes(x=mydata$D2y,y = ..density..,color = "red"),color="darkblue", fill="red",bins=41, alpha=0.5) 

还有一些其他选项，希望不相关：

+ theme_pander() + ggtitle("Histograms") +xlab("Subjective Elicitations")+ylab("Density")+theme(  panel.border = element_blank(),  panel.grid.major = element_blank(),  panel.grid.minor = element_blank(), axis.line = element_line(colour = "black")) + theme(plot.title = element_text(hjust = 0.5))  

生产：

我只想添加一个图例，显示红色直方图代表数据 A，蓝色代表数据 B。

这在 R 中似乎是一项非常困难的任务。

Answer 1

ggplot 设计用于处理长格式数据 - 在这种情况下，图例会自动生成。我建议将您的数据转换为 ggplot 喜欢的格式。您尚未共享任何数据，但根据您的 ggplot 调用，这是我的最佳猜测。如果您需要调试帮助，请分享您的数据的可重现样本。

mydata_long = pivot_longer(mydata, cols = c(D2prediction75At, D2y))

ggplot(mydata_long, aes(x = value, fill = name)) +
  geom_histogram(alpha = 0.65, bins = 41, color = "darkblue") +
  scale_fill_manual(values = c(D2prediction75At = "blue", D2y = "red") +
  ... # your other options

Answer 2

我们也可以

 mydata_long <- melt(mydata,  measure = c("D2prediction75At", "D2y"))
ggplot(mydata_long, aes(x = value, fill = variable)) +
  geom_histogram(alpha = 0.65, bins = 41, color = "darkblue") +
   scale_fill_manual(values = c(D2prediction75At = "blue", D2y = "red") +
   ... # your other options

Answer 3

除了 @GregorThomas 的好建议之外，您还可以尝试下一种方法。首先是使用iris 数据集的示例，您需要在aes() 中传递fill 选项，然后使用scale_fill_manual()：

library(ggplot2)
#Example
ggplot(iris)+
  geom_histogram(aes(x=Sepal.Length,y = ..density..,color = "blue",fill='B'),
                 color="darkblue",bins=41, alpha=0.75)+
  geom_histogram(aes(x=Sepal.Width,y = ..density..,color = "red",fill='A'),
                 color="darkblue",bins=41, alpha=0.5)+
  scale_fill_manual(values=c("blue","red"),'Var')

输出：

对于您的代码，试试这个（未测试，因为没有共享数据）：

#Code OP
ggplot(mydata) + 
  geom_histogram(aes(x=D2prediction75At,y = ..density..,
                     color = "blue",fill='B'),
                 color="darkblue",bins=41, alpha=0.75) +
  geom_histogram(aes(x=D2y,y = ..density..,color = "red",fill='A'),
                 color="darkblue",bins=41, alpha=0.5)+
  scale_fill_manual(values=c("blue","red"),'Data')+
  theme_pander() + 
  ggtitle("Histograms") +
  xlab("Subjective Elicitations")+
  ylab("Density")+
  theme(  panel.border = element_blank(),
          panel.grid.major = element_blank(),
          panel.grid.minor = element_blank(),
          axis.line = element_line(colour = "black")) +
  theme(plot.title = element_text(hjust = 0.5))