我在 R 中有一个图,我试图在 ggplot2 中复制它。我有以下代码:
theta = seq(0,1,length=500)
post <- dgamma(theta,0.5, 1)
plot(theta, post, type = "l", xlab = expression(theta), ylab="density", lty=1, lwd=3)
我试图在 ggplot2 中复制这个情节,这是我能得到的最接近的情节。
df=data_frame(post,theta)
ggplot(data=df,aes(x=theta))+
stat_function(fun=dgamma, args=list(shape=1, scale=.5))