私はRを初めて使用し、ヒストグラムのガンマ分布に適合を追加したいと思います。ガンマ分布をフィットさせてヒストグラムをオーバーレイしたいと思います。
dgamma
関数と関数を使ってガンマ分布を計算することができfitdist
ます。ただし、このガンマ分布をヒストグラムに合わせることができません。
これは私が試したコードです:
hist(mydata, breaks = 30, freq = FALSE, col = "grey")
lines(dgamma(mydata, shape = 1))
私が試したコードは、ガンマ分布の適合をヒストグラムにオーバーレイしていません。フィットなしのヒストグラムしか取得できません。
次の例がオーバーレイに役立つかどうかを確認します
ヒストグラム上。
まず、データセットを作成します。
set.seed(1234) # Make the example reproducible
mydata <- rgamma(100, shape = 1, rate = 1)
次に、ガンマ分布をデータに適合させます。
param <- MASS::fitdistr(mydata, "gamma")
このベクトルは、近似直線に必要です。
x <- seq(min(mydata), max(mydata), length.out = 100)
そして、それらすべてをプロットします。
hist(mydata, breaks = 30, freq = FALSE, col = "grey", ylim = c(0, 1))
curve(dgamma(x, shape = param$estimate[1], rate = param$estimate[2]), add = TRUE)
lines(sort(mydata), dgamma(sort(mydata), shape = 1),
col = "red", lty = "dotted")
この記事はインターネットから収集されたものであり、転載の際にはソースを示してください。
侵害の場合は、連絡してください[email protected]
コメントを追加