Плотность в R участок

Question

Я пытаюсь построить график плотности гамма-распределения.

x<-seq(0,10000,length.out = 1000)
plot(density(rgamma(1000,shape = 7,scale = 120)))
plot(dgamma(x,shape=7,scale=120),col="red")

Но я не понимаю, почему оба графика совершенно разные.

Answer 1

Поскольку вы не указали x в последнем вызове, координаты x по умолчанию соответствуют индексам 1,2,3, ... 1000 вектора dgamma(x,shape=7,scale=120), а не предполагаемым 0,10,20, ....

Если вы сделаете:

x<-seq(0,10000,length.out = 1000)
plot(density(rgamma(1000,shape = 7,scale = 120)))
points(x,dgamma(x,shape=7,scale=120),type = "l", col="red")

Тогда график будет: