R: Диаграмма с накоплением не складывается - PullRequest
0 голосов
/ 18 мая 2018

У меня есть данные, которые я хочу построить как график с областями сложения.На оси X у меня есть данные, которые являются непрерывными, а на оси Y у меня есть непрерывные данные, которые я готовлю для накопления.Это код, который я использую с некоторыми фиктивными данными:

library(data.table)
library(ggplot2)

set.seed(1)
dt <- data.table(var=sample(1:6,1000,replace=TRUE),xdata=runif(1000),ydata=runif(1000))
setorder(dt, var, xdata)

dt$cumydata <- dt[,
                  cumsum(ydata),
                  by = .(var)]$V1/sum(dt$ydata)

ggplot(dt, aes(x = xdata, y = cumydata, fill = as.factor(var))) +
  geom_area(position = "stack")

Вот график вывода: enter image description here

Моя проблема в том, что данныенеправильно составлятьЯ думаю, это может быть из-за непрерывности данных?

Ответы [ 2 ]

0 голосов
/ 04 июня 2018

Так вот, наконец, я так решил, основываясь на информации Джимбоу.Это всего лишь предварительная обработка.Я также сделал все это логарифмическим.

enter image description here

library(data.table)
library(ggplot2)

set.seed(1)
dtt <- data.table(var=sample(1:6,1000,replace=TRUE),xdata=runif(1000),ydata=runif(1000))

setorder(dtt, var, xdata)

log.min.xdata <- log(min(dtt$xdata))
log.max.xdata <- log(max(dtt$xdata))

nbreaks <- 101

temp <- hist(log(dtt$xdata[dtt$var==1]),
             breaks = seq(log.min.xdata, log.max.xdata, length=nbreaks),
             plot = FALSE)


dt <- data.table(var = unlist(lapply(sort(unique(dtt$var)),
                                     function(x){rep(x,nbreaks-1)})),
                 bin = rep(1:(nbreaks-1),length(unique(dtt$var))),
                 mid = rep(temp$mids))

dt$count <- dt[,
               hist(log(dtt$xdata[dtt$var==var]), 
                    breaks = seq(log.min.xdata, log.max.xdata, length=nbreaks),
                    plot = FALSE)$counts,
               by = .(var)]$V1

dt$cumcount <- dt[,
                  cumsum(count),
                  by = .(var)]$V1



pp <- ggplot(dt, aes(x = exp(mid), y = cumcount, fill = as.factor(var))) +
  geom_area(position = "stack") +
  scale_x_log10() +
  theme_bw() +
  theme(legend.position = c(0.1, 0.70),
        legend.background = element_rect(fill="lightgrey", 
                                         size=0.5, linetype="solid")) +
  labs(title = "y",
       fill = " var",
       x = "xdata",
       y = "cumcount") +
  theme(title = element_text(face = "bold"),
        axis.title = element_text(face = "bold"),
        legend.title = element_text(face = "bold"),
        legend.text = element_text(face = "bold"))
0 голосов
/ 18 мая 2018

Для диаграммы с областями с накоплением значения х, а также количество вхождений должны быть одинаковыми.Таким образом, при изменении данных выборки вы получите ожидаемый результат:

set.seed(1)
dt <- data.table(ydata=runif(1002))
dt$var <- rep(1:6, each=167)
dt$xdata <- rep(runif(167), 6)
setorder(dt, var, xdata)

dt$cumydata <- dt[,
                  cumsum(ydata),
                  by = .(var)]$V1/sum(dt$ydata)

ggplot(dt,aes(x = xdata, y = cumydata, fill = as.factor(var))) +
  geom_area(position = "stack")

enter image description here

...