我有一个数据集，数据集归一化为最大值（值介于 [0,1] 之间），我尝试拟合伽马分布。我正在使用 fitdistrplus，我在获得对数似然值以及 AIC 和 BIC 的同时估计分布的参数。使用我的数...

我在数值积分方面遇到了一个问题，无法对后验分布进行归一化。我的代码是： a1 <- 0.012 sigma_a1 <- 0.001 likelihood <- function(a_true) ...