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)

到目前为止,一切都很好。我现在要做的是添加一条线,表示具有参数的 beta 函数alpha=1, beta=3(就像rbeta我使用的函数一样),它考虑了频率而不是密度。中的元素总数rbeta为 365(一年中的天数),我乘以w10 的原因是因为我正在研究的变量可以[0,10]按照上述 beta 分布每天取值。

我该怎么做才能代表这条线?

总而言之,直方图是基于模拟值的,我想展示理论 beta 函数与模拟相比的表现。

4

1 回答 1

0

如果您希望它们匹配,您将需要匹配直方图和密度图曲线下的区域。这应该使它们处于相同的规模。一种方法是

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 回答