如何在R中用伽马分布拟合叠加密度直方图?
我是R新手,希望在直方图中添加伽马分布的拟合。我希望伽马分布适合叠加我的直方图 我能够使用如何在R中用伽马分布拟合叠加密度直方图?,r,histogram,gamma-distribution,R,Histogram,Gamma Distribution,我是R新手,希望在直方图中添加伽马分布的拟合。我希望伽马分布适合叠加我的直方图 我能够使用dgamma函数和fitdist函数计算伽马分布。但是,我无法将此gamma分布作为拟合叠加到直方图上 这是我尝试的代码: hist(mydata, breaks = 30, freq = FALSE, col = "grey") lines(dgamma(mydata, shape = 1)) 我尝试的代码没有将gamma分布拟合叠加到我的直方图上。我只得到没有拟合的直方图。看看下面的例子是否有助于叠
dgamma
函数和fitdist
函数计算伽马分布。但是,我无法将此gamma分布作为拟合叠加到直方图上
这是我尝试的代码:
hist(mydata, breaks = 30, freq = FALSE, col = "grey")
lines(dgamma(mydata, shape = 1))
我尝试的代码没有将gamma分布拟合叠加到我的直方图上。我只得到没有拟合的直方图。看看下面的例子是否有助于叠加
set.seed(1234) # Make the example reproducible
mydata <- rgamma(100, shape = 1, rate = 1)
试试行(sort(mydata)、dgamma(sort(mydata)、shape=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")