9

我正在尝试使用两个叠加的密度图创建直方图。问题：我希望一个密度是一条虚线，它工作得很好，但在图例中，虚线不会出现，如下例所示

x<-sort(rnorm(1000))
data<-data.frame(x=x,Normal=dnorm(x,mean(x),sd=sd(x)),Student=dt(x,df=3))

ggplot(data,aes(y=x))+geom_histogram(aes(x=x,y=..density..),
color="black",fill="darkgrey")+geom_line(aes(x=x,y=Normal,color="Normal"),size=1,
linetype=2)+ylab("")+xlab("")+labs(title="Density estimations")+geom_line(aes(x=x,y=Student,color="Student"),size=1)+
scale_color_manual(values=c("Student"="black","Normal"="black"))

任何想法如何获得图例中的虚线？

非常感谢你！

雷纳

示例图

4

2 回答 2

6

“ggplot”方式通常喜欢数据采用“长”格式，并带有单独的列来指定每种美学。在这种情况下，线型应该被解释为一种美学。处理这个问题的最简单方法是使用reshape2包将数据准备成适当的格式：

library(reshape2)
data.m <- melt(data, measure.vars = c("Normal", "Student"), id.vars = "x")

然后修改您的绘图代码，使其看起来像这样：

ggplot(data,aes(y=x)) +
  geom_histogram(aes(x=x,y=..density..),color="black",fill="darkgrey") +
  geom_line(data = data.m, aes(x = x, y = value, linetype = variable), size = 1) +
  ylab("") +
  xlab("") +
  labs(title="Density estimations")

结果是这样的：

在此处输入图像描述

于 2012-11-26T20:52:42.733 回答
1

您想将其重塑为长格式...使其更简单

x<-sort(rnorm(1000))
Normal=dnorm(x,mean(x),sd=sd(x))
Student=dt(x,df=3)
y= c(Normal,Student)
DistBn= rep(c('Normal', 'Student'), each=1000)
# don't call it 'data' that is an R command
df<-data.frame(x=x,y=y, DistBn=DistBn)

head(df)
          x           y DistBn
1 -2.986430 0.005170920 Normal
2 -2.957834 0.005621358 Normal
3 -2.680157 0.012126747 Normal
4 -2.601635 0.014864165 Normal
5 -2.544302 0.017179353 Normal
6 -2.484082 0.019930239 Normal   



ggplot(df,aes(x=x, y=y))+
  geom_histogram(aes(x=x,y=..density..),color="black",fill="darkgrey")+
  geom_line(aes(x=x,y=y,linetype=DistBn))+
  ylab("")+xlab("")+labs(title="Density estimations")+
  scale_color_manual(values=c("Student"="black","Normal"="black"))

绘图

于 2012-11-26T20:47:02.810 回答