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坐标默认为向量dgamma(x,shape=7,scale=120)的索引1,2,3，... 1000而非预期的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")

然后图表是：