4

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

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)

谢谢！

4

2 回答 2

6

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

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

在此处输入图像描述

于 2013-10-25T21:50:46.743 回答
3

尝试：

curve(dgamma(x, shape=8, scale=0.1), add=TRUE, col='blue')
于 2013-10-25T21:53:43.150 回答