0

我开始说我对 R 真的很陌生,因此这可能是一个非常简单的问题。

我想用 beta 分布做一个蒙特卡罗模拟。情况是我有一个网站,该网站将在一年内每天收到 0 到 10 个可视化,我想对此进行模拟。

我论文的第一步是制作一个显示 beta 分布形状的图表。

为此,我愿意

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]:

第二张图

谁能帮我?提前谢谢你的帮助!

回答部分:通过使用 plot(density(10*x), xlim = c(0, 10)) 我获得 答案 1 了它再次具有域偏移量的 ptblem 并且不仅看起来函数被移动到对,因为我希望给定 x1

4

1 回答 1

0

你真的只是想要一个带有参数 1、4 的 Beta 分布图吗?:

# 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")

在此处输入图像描述

如果你真的想要模拟的东西,你也可以做一个这样的直方图:

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

在此处输入图像描述

于 2014-10-02T15:44:05.980 回答