我想创建一个函数，以便针对x集成2个分布的乘积（PDF和CDF）。这是为了进行概率风险评估并确定“预期总风险”。
我到目前为止无法使用的代码是...
F <- function(x) {dgamma(x, shape=0.4259325, rate=8.1490741)*plnorm(x, meanlog=3.4906738, sdlog=0.2938556)}
c <- integrate( F, lower = 0, upper = inf)
我在哪里哪里做错了什么？
此外，我已经使用'fitdist'函数拟合了我的发行版，现在我很难对发行版参数进行编码，以尝试使集成工作正常进行。但是，我如何直接从拟合的分布中获取形状/速率/均值等-下面的示例代码中的“拟合”对象？
x_vect <- seq(0.1,10, length.out = 1e3)
example <- rgamma(x_vect, shape=0.4259325, rate=8.1490741)
plot(x_vect, example, type = "l", log="x")
Fitted <- fitdist(example, "gamma")
plot(Fitted)