Построение функции плотности гамма-вероятности

Question

Я пытаюсь построить функцию плотности гамма-вероятности в R, где y∈ (0,10) для (k = 1, μ = 1), (k = 2, μ = 1), (k = 2, μ = 2). В R

В R функция pgamma принимает:

pgamma(q, shape, rate = 1, scale = 1/rate, alpha = shape, beta = scale, lower.tail = TRUE, log.p = FALSE)

В R я попробовал:

pgamma(1,1,rate=1,scale = 1/rate, alpha = shape, beta = scale, lower.tail = True, log.p = False)

Но я получаю сообщение

Error in pgamma(1, 1, rate = 1, scale = 1/rate, lower.tail = TRUE, log.p = FALSE) : 
object 'rate' not found

Я впервые строю гамма-распределение и буду признателен за помощь.

Answer 1

На следующем графике показаны три плотности с использованием базовой графики R.

Во-первых, значения параметров, которые вы хотите. Я предполагаю, что ваш mu соответствует определению на странице Википедии гамма-дистрибутива .

k <- c(1, 2, 2)
mu <- c(1, 1, 2)
theta <- mu/k

Теперь участки.

plot(0, 0, xlim = c(0, 10), ylim = c(0, 1), type = "n")
for(i in seq_along(k))
  curve(dgamma(x, shape = k[i], scale = theta[i]), from = 0, to = 10, col = i, add = TRUE)