我想创建一个函数，以便针对 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)