使用freq = TRUE表示直方图中的频率线

时间:2014-10-02 17:54:44

标签: r simulation beta-distribution

我在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的beta函数的线(如我使用的rbeta函数),它考虑了频率而不是密度。 rbeta中的元素总数为365(一年中的天数)以及我将w乘以10的原因是因为我正在研究的变量可以假设每个值[0,10]一天,遵循上述beta分布。

我需要做些什么才能代表这条线?

总结一下,直方图是基于模拟值的,我想说明理论beta函数与模拟相比会如何表现。

1 个答案:

答案 0 :(得分: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)

这给出了

enter image description here