Как создать круговой график, показывающий ежемесячное присутствие или отсутствие, используя radio.plot/ggplot2

введите здесь описание изображения

Я хочу создать фигуру (похожую на приложенную) в ggplot2, которая является круглой с месяцем года, представленным на оси, поэтому двенадцать делений по окружности, с одной внутренней дугой определенного цвета, соответствующей месяцам, где pa = 1 ( другие варианты: pa = 0) и другая внутренняя дуга (вне первой дуги), соответствующая sampled = 1 (другой вариант: sampled = 0). Данные выглядят так:

 m
   season pa month sampled occupancy
1  spring  1     3       1   present
2  spring  1     4       1   present
3  spring  1     5       1   present
4  summer  1     6       1   present
5  summer  1     7       1   present
6  summer  1     8       1   present
7  winter  0    12       1    absent
8  winter  0     1       0    absent
9  winter  0     2       0    absent
10   fall  1     9       1   present
11   fall  1    10       1   present
12   fall  1    11       1   present

Раньше я использовал ggplot(m, aes(season, fill=season)) с geom_bar(width=1) и coord_polar(), но это просто дает мне круговую диаграмму.

Теперь пытаюсь:

radial.plot(m, radial.pos=c(1,2,3,4,5,6,7,8,9,10,11,12), labels="", rp.type="r",
            line.col=8, lwd=3, add=TRUE)

и получаю ошибку

plot.new еще не вызывался

Я думаю, что неправильно понимаю, какой ввод нужен radial.plot, и потенциально использую неправильную функцию для желаемого результата.


person Community    schedule 10.12.2016    source источник
comment
Не могли бы вы опубликовать пример сюжета, который вы хотите создать (даже если сделать это быстро и вручную в Paint или Gimp)?   -  person zoneparser    schedule 11.12.2016
comment
@zoneparser Я прикрепил фото ... есть идеи? Благодарю вас!   -  person    schedule 11.12.2016
comment
Похоже, вам может понадобиться использовать geom_rect, поскольку все объекты, которые вы пытаетесь построить, являются прямоугольниками (в полярных координатах).   -  person Mist    schedule 12.12.2016


Ответы (2)


Я не был уверен, как связать ваши данные с графиком выше, но поможет ли следующий код?

df <- data.frame(startdate  = as.Date(c("2016-02-20", "2016-02-20", "2016-02-20")),
             finishdate = as.Date(c("2016-04-30", "2016-04-30", "2016-06-10")),
             y          = c(4,5,8))


library(ggplot2)

ggplot(df) +
  geom_rect(aes(xmin = startdate, 
                xmax = finishdate,
                ymin = y-0.2,
                ymax = y + 0.2)) +
  geom_rect(aes(xmin = startdate - 5, 
                xmax = finishdate + 5,
                ymin = y-0.05,
                ymax = y + 0.05)) +
  xlim(as.Date("2016-01-01"), as.Date("2016-12-31")) +
  ylim(0,10) +
  coord_polar()

введите здесь описание изображения

person Mist    schedule 12.12.2016

Вот два разных подхода в ggplot2 для преобразования предоставленных данных в полярный график. Эти два подхода отличаются способом обработки month.

Чтение данных

library(data.table)
m <- fread("id   season pa month sampled occupancy
           1  spring  1     3       1   present
           2  spring  1     4       1   present
           3  spring  1     5       1   present
           4  summer  1     6       1   present
           5  summer  1     7       1   present
           6  summer  1     8       1   present
           7  winter  0    12       1    absent
           8  winter  0     1       0    absent
           9  winter  0     2       0    absent
           10   fall  1     9       1   present
           11   fall  1    10       1   present
           12   fall  1    11       1   present")

Подготовить данные

# reshape from wide to long (as preferred by ggplot)
ml <- melt(m, measure.vars = c("pa", "sampled"))
# create factors to ensure desired order
ml[, variable := factor(variable, levels = c("pa", "sampled"))]

Вариант 1: месяцы как фактор

ml[, fct_month := factor(month, levels = 1:12, labels = month.abb)]
library(ggplot2)
ggplot(ml[value != 0], aes(x = fct_month, y = variable, 
                           group = variable, colour = variable)) + 
  geom_line(size = 5) +
  scale_y_discrete(expand = c(0, 1), breaks = NULL) +
  xlim(month.abb) + 
  coord_polar() +
  theme_bw() + xlab(NULL) + ylab(NULL)

введите описание изображения здесь

Вариант 2: месяцы как период (с датой начала и окончания)

ml[, start := as.Date("2016-01-01") + base::months(month - 1L)]
# last day of month = begin of next month - 1 day
ml[, end := as.Date("2016-01-01") + base::months(month) - 1L]
library(ggplot2)
ggplot(ml, aes(x = start, xend = end, y = variable, yend = variable,
               group = variable, colour = variable,
               size = value)) + 
  geom_segment() +
  scale_y_discrete(expand = c(0, 1), breaks = NULL) +
  scale_x_date(date_breaks = "1 month", date_labels = month.abb,
               limits = range(c(ml$start, ml$end))) +
  scale_size(guide = FALSE) +
  coord_polar() +
  theme_bw() + xlab(NULL) + ylab(NULL)

введите описание изображения здесь

Вариант 3: цветной фон

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

ml[, start_year := as.Date("2016-01-01")]
# last day of month = begin of next month - 1 day
ml[, end_year := as.Date("2016-12-31")]
ml[, start := start_year + base::months(month - 1L)]
ml[, end := start_year + base::months(month) - 1L]

library(ggplot2)
bg_height <- 1.0
ggplot(ml) + 
  geom_rect(aes(xmin = start_year, xmax = end_year, 
                ymin = as.integer(variable) - 0.5 * bg_height, 
                ymax = as.integer(variable) + 0.5 * bg_height,
                group = variable, fill = variable)
            ) + 
  geom_segment(aes(x = start, xend = end, y = variable, yend = variable,
                   group = variable, size = value),
               colour = "gray20") +
  scale_y_discrete(expand = c(0, 1), breaks = NULL) +
  scale_x_date(date_breaks = "1 month", date_labels = month.abb,
               limits = range(c(ml$start, ml$end))) +
  scale_size(guide = FALSE) +
  scale_fill_brewer(palette = "Blues") +
  coord_polar() +
  theme_bw() + xlab(NULL) + ylab(NULL)

введите описание изображения здесь

person Uwe    schedule 12.12.2016
comment
Можно ли добавить круглый цветной фон, например, panel.background = element_rect(fill='grey'), а затем shape = 1 (я знаю, что вы не можете добавить фигуру к прямоугольнику элемента, но, может быть, есть способ сделать это?) - person ; 18.12.2016
comment
Вы можете попытаться нарисовать geom_rect() с правильно выбранными шириной, высотой и цветом в качестве первых слоев, чтобы создать белые, голубые и синие кольца, которые образуют фон на снимке экрана. - person Uwe; 18.12.2016
comment
Я добавил 3-й вариант с цветными фоновыми кольцами, как и просили. - person Uwe; 18.12.2016