2014-10-02 9 views
0

У меня есть следующий фрагмент кода в R:представляют частоты линии в histagram с помощью FREQ = 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), которая учитывает частоту, а не плотность. Общее количество элементов в rbeta составляет 365 (дни в году), и причина, по которой я умножаю w на 10, состоит в том, что переменная, которую я изучаю, может принимать значение [0,10] каждый день, следуя описанному выше бета-распределению.

Что мне нужно сделать, чтобы представить эту строку?

Подводя итог, гистограмма основана на симулированных значениях, и я хочу показать, как теоретическая бета-функция вели бы по сравнению с симуляцией.

+0

'dbeta' является PDF для бета-распределения. –

+0

Мне это известно, но когда я добавляю его в график, его начальное значение равно 3, а гистограмма с частотой = ИСТИНА начинается с 45, а гистограмма с частотой = FALSE начинается с 0,25 – Pezze

ответ

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

 Смежные вопросы

  • Нет связанных вопросов^_^