Как добавить метки и проценты к печатной форме R sunburst

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

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

Вот мой код.

# libraries
library(dplyr)
library(treemap)
library(sunburstR)
library(readxl)
library(vcd)


## Load Arthritis as example
Data <- data.frame(Arthritis)
Data <- Data %>% select(-ID) %>% 
mutate(Age=ifelse(Age<50,"Young","Old")) %>% group_by(Treatment,Sex,Improved,Age) %>% 
summarise(Count=n()) %>% 
mutate(Path=paste(Treatment,Sex,Improved,Age,sep="-")) %>% 
ungroup() %>% 
select(Path,Count)

sunburst(Data)

Любая помощь была бы замечательной.

Спасибо.


r sunburst-diagram
person Wilze    schedule 17.03.2019    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Я предлагаю пакет ggsunburst https://github.com/didacs/ggsunburst

library(ggsunburst)
library(dplyr)
library(vcd) # just for the Arthritis dataset
Data <- data.frame(Arthritis)

# compute percentage using tally
# add column leaf, with format "name->attribute:value"
# ggsunburst considers everything after "->" as attributes
# the attribute "size" is used as the size of the arc
df <- Data %>% 
  mutate(Age=ifelse(Age<50,"Young","Old")) %>%
  group_by(Treatment,Sex,Improved,Age) %>% 
  tally() %>%
  mutate(percentage = n/nrow(Data)*100, 
         size=paste("->size:",round(percentage,2),sep=""),
         leaf=paste(Improved,size,sep = "")) %>%
  ungroup() %>%
  select(Treatment,Sex,Age,leaf)

# sunburst_data reads from a file so you need to create one
write.table(df, file = 'data.csv', row.names = F, col.names = F, sep = ",")

# specify node_attributes = "size" to add labels with percentages in terminal nodes
sb <- sunburst_data('data.csv', type = "lineage", sep = ',', node_attributes = "size")

# compute percentages for internal nodes
tre <- Data %>%
  group_by(Treatment) %>%
  tally() %>%
  mutate(percent=n/nrow(Data)*100,
         name=Treatment) %>%
  ungroup() %>%
  select(name,percent)

sex <- Data %>%
  group_by(Treatment,Sex) %>%
  tally() %>%
  mutate(percent=n/nrow(Data)*100,
         name=Sex) %>%
  ungroup() %>%
  select(name,percent)

age <- Data %>%
  mutate(Age=ifelse(Age<50,"Young","Old")) %>%
  group_by(Treatment,Sex,Age) %>%
  tally() %>%
  mutate(percent=n/nrow(Data)*100,
         name=Age) %>%
  ungroup() %>%
  select(name,percent)

x <- rbind(tre, sex, age)
# the rows in x are in the same order as sb$node_labels, cbind works here only because of that
x <- cbind(sb$node_labels, round(x[,"percent"],2))
percent <- x %>% mutate(name_percent = paste(label,percent,"%"))

sunburst(sb, node_labels.min = 0) +
  geom_text(data = sb$leaf_labels, aes(x=x, y=0.1, label=paste(size,"%"), angle=angle, hjust=hjust), size = 2) +
  geom_text(data = percent, aes(x=x, y=y, label=name_percent, angle=pangle), size=2)

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

person didac    schedule 23.02.2020