如何在直方图上绘制真密度(对于指定的伽玛分布)

时间:2013-10-25 21:42:09

标签: r

我想在直方图上绘制真密度。 这是我的代码,我不确定我的代码是否正确。

set.seed(600)
x <- rgamma(500,shape=8,scale=0.1)
mean(x)
hist(x,prob=T,main='Gamma,scale=0.1')
lines(density(x),col='red',lwd=2)

谢谢!

2 个答案:

答案 0 :(得分:5)

density()给出了密度的局部估计值。 true 密度可以从dgamma导出,并按此绘制......

lines( sort(x) , y = dgamma( sort(x) , shape = 8 , scale = 0.1 ) , col = "blue" , lty = 2 , lwd = 2 )

enter image description here

答案 1 :(得分:3)

尝试:

curve(dgamma(x, shape=8, scale=0.1), add=TRUE, col='blue')