Гистограмма с весами в R

Question

Мне нужно построить взвешенную гистограмму плотности, а не частоту. Я знаю, что freq = FALSE доступен в hist(), но вы не можете указать вес. В ggplot2 я могу это сделать:

library(ggplot2) 
w <- seq(1,1000) 
w <-w/sum(w) 
v <- sort(runif(1000)) 

foo <- data.frame(v, w) 

ggplot(foo, aes(v, weight = w)) + geom_histogram() 

Но где эквивалент freq = FALSE?

Answer 1

По умолчанию geom_histogram() будет использовать частоту, а не плотность по оси y. Тем не менее, вы можете изменить это, установив y эстетику в ..density.. так:

ggplot(foo, aes(x = v, y = ..density.., weight = w)) + geom_histogram() 

Это даст взвешенную гистограмму v с плотностью на оси у.

Вы также можете сделать это с freq аргументом в weighted.hist() из plotrix пакета:

library(plotrix) 
with(foo, weighted.hist(v, w, freq = FALSE))