ggplot смесь модель R

Question

У меня есть набор данных с числовыми значениями и категориальной переменной. Распределение числовой переменной отличается для каждой категории. Я хочу построить график плотности для каждой категориальной переменной, чтобы они были визуально ниже всего графика плотности.

Это похоже на компоненты модели смеси без расчета модели смеси (поскольку я уже знаю категориальную переменную, которая разделяет данные).

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

library(ggplot2) 
ggplot(iris, aes(x = Sepal.Width)) + geom_density() + geom_density(aes(x = Sepal.Width, group = Species, colour = 'Species')) 

То, что я хочу, чтобы иметь плотность каждой категории, как к югу от плотности (не интегрирующей 1). Похож на следующий код (который я только реализован для двух из трех видов ирисов)

myIris <- as.data.table(iris) 
# calculate density for entire dataset 
dens_entire <- density(myIris[, Sepal.Width], cut = 0) 
dens_e <- data.table(x = dens_entire[[1]], y = dens_entire[[2]]) 

# calculate density for dataset with setosa 
dens_setosa <- density(myIris[Species == 'setosa', Sepal.Width], cut = 0) 
dens_sa <- data.table(x = dens_setosa[[1]], y = dens_setosa[[2]]) 

# calculate density for dataset with versicolor 
dens_versicolor <- density(myIris[Species == 'versicolor', Sepal.Width], cut = 0) 
dens_v <- data.table(x = dens_versicolor[[1]], y = dens_versicolor[[2]]) 

# plot densities as mixture model 
ggplot(dens_e, aes(x=x, y=y)) + geom_line() + geom_line(data = dens_sa, aes(x = x, y = y/2.5, colour = 'setosa')) + 
    geom_line(data = dens_v, aes(x = x, y = y/1.65, colour = 'versicolor')) 

, в результате чего

Выше я жестко закодировано число для уменьшения значения у. Есть ли способ сделать это с помощью ggplot? Или рассчитать?

Спасибо за ваши идеи.

Answer 1

Вы имеете в виду что-то вроде этого? Однако вам нужно изменить масштаб.

ggplot(iris, aes(x = Sepal.Width)) + 
    geom_density(aes(y = ..count..)) + 
    geom_density(aes(x = Sepal.Width, y = ..count.., 
       group = Species, colour = Species)) 

Другим вариантом может быть

ggplot(iris, aes(x = Sepal.Width)) + 
    geom_density(aes(y = ..density..)) + 
    geom_density(aes(x = Sepal.Width, y = ..density../3, 
        group = Species, colour = Species))