我在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函数与模拟相比会如何表现。
答案 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)
这给出了