r中直方图中的垂直线

r中直方图中的垂直线,r,ggplot2,R,Ggplot2,我正在用R中的和平代码进行一些尝试。我试图在同一个图形中创建6个不同的直方图。它工作正常,但我需要在6个直方图中的每个直方图中放置1条垂直线。我正在使用的代码可能如下所示: require(ggplot2) require(reshape2) require(gdata) MC_beta=rbind(rnorm(1000,-2,0.1),rnorm(1000,-1,0.1),rnorm(1000,0,0.1),rnorm(1000,0.5,0.1),rnorm(1000,1,0.1),rnor

我正在用R中的和平代码进行一些尝试。我试图在同一个图形中创建6个不同的直方图。它工作正常,但我需要在6个直方图中的每个直方图中放置1条垂直线。我正在使用的代码可能如下所示:

require(ggplot2)
require(reshape2)
require(gdata)

MC_beta=rbind(rnorm(1000,-2,0.1),rnorm(1000,-1,0.1),rnorm(1000,0,0.1),rnorm(1000,0.5,0.1),rnorm(1000,1,0.1),rnorm(1000,2,0.1))


df <- data.frame(MC_beta[1,], MC_beta[2,], MC_beta[3,], MC_beta[4,],MC_beta[5,],MC_beta[6,])
names(df)[1:6]<-c("1", "2", "3", "4","5","6")

df2 = melt(df)

z=c(-2,-1,0,0.5,1,2)

ggplot(df2, aes(x=value, fill = variable)) + geom_vline(xintercept = z, colour="black") +geom_histogram(binwidth=0.03,colour = "black") + scale_fill_manual(name = "",values = c('red','blue',"red","blue","red","blue")) +
  facet_wrap(~variable,nrow=6, ncol=1) + scale_x_continuous(breaks=seq(-2.5,2.5,0.5)) + guides(fill=FALSE) +
  theme_bw() + theme(strip.background = element_blank(),axis.text=element_text(size=14.5),strip.text.x = element_text(size = 14.5)) + stat_function(fun = dnorm)
require(ggplot2)
要求(2)
要求(gdata)
MC_β=rbind(rnorm(1000,-2,0.1),rnorm(1000,-1,0.1),rnorm(1000,0,0.1),rnorm(1000,0.5,0.1),rnorm(1000,1,0.1),rnorm(1000,2,0.1))

df您的
z
需要是一个
data.frame
,对于定义刻面的变量的每个值,都有相应的
xintercept
。尝试以下更改:

z <- data.frame(variable=levels(df2$variable), 
                mean=c(-2,-1,0,0.5,1,2))

ggplot(df2, aes(x=value, fill = variable))+ 
  geom_vline(data=z, aes(xintercept = mean), colour="black") +
  geom_histogram(binwidth=0.03,colour = "black") + 
  scale_fill_manual(name = "",values = c('red','blue',"red","blue","red","blue")) +
  facet_wrap(~variable,nrow=6, ncol=1) + 
  scale_x_continuous(breaks=seq(-2.5,2.5,0.5))+ guides(fill=FALSE) +
  theme_bw() + 
  theme(strip.background = element_blank(), axis.text=element_text(size=14.5), strip.text.x = element_text(size = 14.5)) + 
  stat_function(fun = dnorm)

z您在数据外部有
z
,因此您将在每个方面绘制一条垂直线。使用

df2 <- (merge(df2, cbind.data.frame(variable=names(df), z)))

谢谢,看起来很好。很好,+1。但是,这将多次绘制垂直线。它可能无法区分,但如果是文本,它看起来会很奇怪。
geom_vline(aes(xintercept = z), colour="black")