Как нарисовать истинную плотность (для указанного гамма-распределения) на гистограмме

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) 

Answer 2

Try:

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