2014-10-02 83 views
0

我有R中的下面的代碼段:表示使用頻率在histagram頻率線= TRUE

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(如我所用的rbeta函數)的beta函數的線,該線考慮了頻率而不是密度。 rbeta中的元素總數爲365(一年中的天數),並且我將w乘以10的原因是因爲我正在研究的變量可以假設值爲每天[0,10],遵循上述的beta分佈。

我需要做些什麼來表示這條線?

總結,直方圖是基於模擬值,我想展示理論beta函數與模擬相比的表現。

+0

'dbeta'是beta分佈的PDF文件。 – 2014-10-02 18:10:52

+0

我意識到這一點,但是當我將它添加到圖中時,其初始值等於3,而freq = TRUE的直方圖從45開始,並且freq = FALSE的直方圖從0.25開始 – Pezze 2014-10-02 18:13:04

回答

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