亲爱的 我使用以下代码将对数正态分布拟合到我的一组正数据的直方图中。当我运行它时,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")