我想在直方图上绘制真密度。 这是我的代码,我不确定我的代码是否正确。
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)
谢谢!
答案 0 :(得分:5)
density()
给出了密度的局部估计值。 true 密度可以从dgamma
导出,并按此绘制......
lines( sort(x) , y = dgamma( sort(x) , shape = 8 , scale = 0.1 ) , col = "blue" , lty = 2 , lwd = 2 )
答案 1 :(得分:3)
尝试:
curve(dgamma(x, shape=8, scale=0.1), add=TRUE, col='blue')