[解決済み] ガンマー確率密度関数のプロット
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)
関連
-
[解決済み】ggplotの線幅を変更するには?
-
[解決済み】library(ggplot2)でエラー:'ggplot2'というパッケージは存在しません。
-
[解決済み】RでのMLEエラー:'vmmin'の初期値が有限でない
-
[解決済み】「次のオブジェクトは 'package:xxx' からマスクされています」とはどういう意味ですか?
-
[解決済み】ggplot boxplotでPosition-dodge警告?
-
[解決済み】R ggplot2 で scale_x_discrete を使用する。
-
[解決済み】Rの整数オーバーフローとは何ですか、そしてどのように起こるのですか?
-
[解決済み】r Error dim(X) must have a positive length?
-
[解決済み] [Solved] read.csv warning 'EOF within quoted string' prevents complete reading of file.
-
[解決済み】seq_alongはうまくいくが、seqが意図しない結果を生む例とは?
最新
-
nginxです。[emerg] 0.0.0.0:80 への bind() に失敗しました (98: アドレスは既に使用中です)
-
htmlページでギリシャ文字を使うには
-
ピュアhtml+cssでの要素読み込み効果
-
純粋なhtml + cssで五輪を実現するサンプルコード
-
ナビゲーションバー・ドロップダウンメニューのHTML+CSSサンプルコード
-
タイピング効果を実現するピュアhtml+css
-
htmlの選択ボックスのプレースホルダー作成に関する質問
-
html css3 伸縮しない 画像表示効果
-
トップナビゲーションバーメニュー作成用HTML+CSS
-
html+css 実装 サイバーパンク風ボタン
おすすめ
-
[解決済み】 'height' はベクトルか行列でなければならない barplot のエラー
-
[解決済み】rbind エラー。"名前が以前の名前と一致しない"
-
[解決済み】lis[[i]]のエラー:1つ未満の要素を選択しようとした場合
-
[解決済み] write.tableしようとすると、未実装の型リストが表示される。
-
[解決済み】Rで立方根と対数変換をする
-
[解決済み】Rで文字ベクトルから引用符を削除する
-
[解決済み】R: predict() の数値 'envir' 引数が長さ1でない。)
-
[解決済み】エラー - replacement has [x] rows, data has [y].
-
[解決済み】forループを実行すると「要因のレベルセットが異なる」というエラーが発生する
-
[解決済み】'NULL'型の非(リストまたはベクトル)に適用されるis.na()は何を意味するのか?