我尝试使用fitdistrplus::fitdistr
将截断的正态分布拟合到数据并指定上限和下限。然而,当将MLE拟合参数与没有界限的MLE拟合参数进行比较时,它们看起来是相同的。
library(fitdistrplus)
library(MASS)
dt <- rnorm(100, 1, 0.5)
cat("truncated:", fitdistr(dt, "normal", lower = 0, upper = 1.5, method = "mle")$estimate,
"original:", fitdist(dt, "norm", method = "mle")$estimate, sep = "\n")
truncated:
1.034495
0.4112629
original:
1.034495
0.4112629
我不是统计天才,但我很确定参数应该是不同的,因为截断分布,均值和sd都会改变(因为分布重新调整)。这是对的吗?
感谢您的建议
干杯,
西蒙