对数正态拟合的密度直方图

时间:2019-08-30 14:37:27

标签: data-fitting

亲爱的 我使用以下代码将对数正态分布拟合到我的一组正数据的直方图中。当我运行它时,Rstudio不会显示任何错误消息,但是同时,没有结果,也没有生成图形。

怎么回事?

data("BoreholesR", package = "MASS")
my_data <- BoreholesR$UCS

x <- seq(0,max(my_data), length=100)
hst <- hist(my_data, breaks=x)
fit_params <- fitdistr(my_data,"lognormal")
fit <- dlnorm(x, fit_params$estimate['meanlog'], 
fit_params$estimate['sdlog'])
plot(x, fit, type="l", ylab="Density",
xlab="UCS", ylim=c(0,max(hst$density)), xlim=c(0,100))
title(main = "Density histogram with lognormal fit")
lines(hst$mid, hst$density, type="l", col="red")

0 个答案:

没有答案