0

Rに次のコードがあります：

w=rbeta(365,1,3,ncp=0)
hist(10*w,breaks=25,freq=TRUE,xlim=c(0,10), ylim=c(0,60))
h=seq(0,1,0.05)

ここまでは順調ですね。ここでやりたいことは、alpha=1, beta=3(rbeta使用した関数のように) パラメーターを持つベータ関数を表す行を追加することです。これは、密度ではなく周波数を考慮に入れます。の要素の総数rbetaは 365 (1 年の日数) であり、w10 を掛ける理由は、私が調査している変数が、[0,10]上記のベータ分布に従って毎日値を想定できるためです。

この線を表すにはどうすればよいですか?

要約すると、ヒストグラムはシミュレートされた値に基づいており、シミュレーションと比較して理論上のベータ関数がどのように動作したかを示したいと思います。

4

1 に答える 1

0

それらを一致させたい場合は、ヒストグラムと密度プロットの曲線の下の領域を一致させたいと思うでしょう. それはそれらを同じスケールに置くべきです。それを行う1つの方法は

set.seed(15) #to make it reproducible
w <- rbeta(365, 1, 3, ncp=0)
hh <- hist(w*10, breaks=25, freq=TRUE, xlim=c(0,10), ylim=c(0,60))
ss <- sum(diff(hh$breaks)*hh$counts)
curve(dbeta(x/10, 1, 3, ncp=0)*ss/10, add=T)

これは与える

ここに画像の説明を入力

于 2014-10-02T20:09:37.763 に答える