1. ホーム
  2. r

[解決済み] ガンマー確率密度関数のプロット

2022-02-28 22:04:18

質問

Rで、(k = 1,μ = 1), (k = 2, μ = 1), (k = 2, μ = 2)のy∈(0,10)のガンマ確率密度関数をプロットしようとしているのですが、どうすればいいですか?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

ガンマ分布のプロットは初めてなのですが、何かご教授いただけないでしょうか。

どのように解決するのですか?

以下は、基本的なRグラフィックスを使用して3つの密度をプロットしたものです。

まず、欲しいパラメータ値。私が想定しているのは、あなたの mu で定義されているとおりです。 ガンマ分布のWikipediaページ .

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)