Я хочу вычислить интегральную функцию распределения в R для данных, следующих за гамма-распределением. Я понял, как сделать это с логнормальным распределением, используя уравнение из Википедии; однако уравнение гаммы кажется более сложным, и я решил использовать функцию pgamma ().
Я новичок в этом и не понимаю следующее:
- Почему я получаю три разных значения из pgamma, и как имеет смысл, что они отрицательные?
- Должен ли я брать журнал всех квантилей, так же как я использовал log (среднее значение) и log (стандартное отклонение) при выполнении расчетов с распределением lognorm?
- Как я концептуально понимаю CDF, рассчитанный по pgamma? Для lognorm имело смысл, что я вычислял вероятность того, что X примет значение <= x, но в этой функции pgamma нет "x". </li>
Действительно ценю помощь в понимании этого.
shape <- 1.35721347
scale <- 1/0.01395087
quantiles <- c(3.376354, 3.929347, 4.462594)
pgamma(quantiles, shape = shape, scale = scale, log.p = TRUE)