r - 如何在直方图上绘制真实密度(对于指定的 Gamma 分布)

Question

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

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)

谢谢!

Answer 1

density()给出密度的局部估计。真密度可以从 dgamma 推导出来，并像这样绘制......

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