2015-02-23 4 views
0

Просто быстрый вопрос: Я пытаюсь добавить строку плотности моей гистограммыкоманда «LINES» в R

lines(x,dbeta(x,6,3)) 

Я определено ранее:

x = runif(100000,0,1) 

Однако довольно чем «присоединение к точкам» бета-дистрибутива, оно «заполняется цветом», все под ним. Поэтому я в основном остался с цветным бета-распределением, наложенным на мою гистограмму.

Я попытался изменить типы линий, но это, похоже, не помогает. Как я могу просто превратить его в однострочный оверлейный сюжет?

ответ

2

вы замышляете пары случайных чисел с его соответствующей плотностью ... если вы просто хотите, чтобы перекрывать линию, не используйте runif() , но seq():

x <- seq(from=0, to=1, length.out=10000) 
lines(x, dbeta(x, 6,3)) 
1

Вот решение с использованием ggplot2

x = runif(100000,0,1) 
x = data.frame(beta=dbeta(x,6,3)) 
library(ggplot2) 

ggplot(x, aes(beta)) + 
    geom_histogram(aes(y=..density..), 
       binwidth=.1, 
       colour="red", fill="white") + 
    geom_density(alpha=.2) 

enter image description here