0

私はRに本当に慣れていないので、これは非常に単純な問題である可能性があると言い始めます。

ベータ分布でモンテカルロ シミュレーションを実行したいと考えています。場合は、1 年間毎日 0 から 10 のビジュアライゼーションを受け取る Web サイトがあり、これをシミュレートしたいと考えています。

私の論文の最初のステップは、ベータ分布の形状を示す図を作成することです。

そうするために、私はします

x=rbeta(365,1,4,ncp=0)
y=dbeta(x,1,4)
plot(10*x,y,type="l")

次のような醜い図が表示されます。 最初の図

それから私はやろうとした

x=rbeta(365,1,4,ncp=0)
plot(density(10*x))

この場合、この欠陥のある右の図を取得します。0 側でドメイン [0,10] を超えています。

2番目の図

誰でも私を助けることができますか?事前に助けてくれてありがとう!

回答部分: plot(density(10*x), xlim = c(0, 10)) 答え 1 を使用して、ドメインのオフセットを持つ ptoblem を再度取得するだけでなく、関数が正しい、与えられた x1 を期待しているので

4

1 に答える 1

0

Do you literally just want a plot of the Beta distribution with parameters 1, 4?:

# Points over which we will evaluate Beta PDF
x <- seq(0, 1, 0.05)

# PDF values
y <- dbeta(x, 1, 4, ncp=0)

plot(x*10, y, type="l")

enter image description here

If you really want something simulated you could also do a histogram like this:

hist(rbeta(2000, 1, 4, ncp=0), breaks=25)

enter image description here

于 2014-10-02T15:44:05.980 に答える