Создать функцию для интегрирования 2 распределений по x
Я хотел бы создать функцию, чтобы интегрировать произведение двух распределений (PDF и CDF) по x. Это делается для того, чтобы выполнить вероятностную оценку риска и определить "ожидаемый общий риск".
Мой код у меня до сих пор, который не работает, это...
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)