如何用R中的伽马分布覆盖密度直方图?

时间:2019-05-26 11:40:07

标签: r histogram gamma-distribution

我是R的新手,想对直方图的伽马分布添加拟合。我希望伽玛分布适合叠加我的直方图。

我可以使用WHERE name = ...函数以及mysqli_fetch_array函数来计算伽玛分布。但是,我无法将这种伽马分布叠加到我的直方图上。

这是我尝试的代码:

dgamma

我尝试的代码没有将适合我的直方图的伽马分布覆盖。我只得到不适合的直方图。

1 个答案:

答案 0 :(得分:0)

查看下面的示例是否可以帮助覆盖

  1. 一条黑色的修身线
  2. 带红色虚线的PDF图形

在直方图上。

首先,创建一个数据集。

set.seed(1234)    # Make the example reproducible
mydata <- rgamma(100, shape = 1, rate = 1)

现在使伽马分布适合数据。

param <- MASS::fitdistr(mydata, "gamma")

拟合线需要此向量。

x <- seq(min(mydata), max(mydata), length.out = 100)

并全部绘制。

hist(mydata, breaks = 30, freq = FALSE, col = "grey", ylim = c(0, 1))
curve(dgamma(x, shape = param$estimate[1], rate = param$estimate[2]), add = TRUE)
lines(sort(mydata), dgamma(sort(mydata), shape = 1),
      col = "red", lty = "dotted")

https://www.dropbox.com/s/p1bhacdg8j1qx42/data.csv?dl=0