Но чтобы добавить значения счетчика, вычисленные в geom_histogram, сверху столбцов в R - PullRequest
Новая
       51

Но чтобы добавить значения счетчика, вычисленные в geom_histogram, сверху столбцов в R

0 голосов
/ 10 октября 2019

Я хотел бы добавить значения счетчика, вычисленные в функции geom_histogram на ggplot2. Я поместил ggplot2 в цикл, чтобы я мог создавать несколько графиков, в моем случае 30, но для простоты, вот фиктивный набор только для четырех графиков. Обтекание фасет не сработало, поскольку плотность геом объединяла данные по всем факторам до вычисления пропорций, а не внутри фактора / переменной. Чтобы создать этот сюжет, я, по сути, смешал целую кучу кода из разных источников, так что поверьте им. 

              library(dplyr)
              library(ggplot2)
              library(ggridges)
              library(reshape2)
              library(gridExtra)
       #Make the data#       
        df.fact <- data.frame("A"=rnorm(400, mean = 350, sd=160),"B"=rnorm(400, mean = 300, sd=100), "C"=rnorm(400, mean = 200, sd=80), names=rep(factor(LETTERS[23:26]), 100))
        df.test<-melt(df.fact, id.vars = "names", value.name = "Length2")
        names(df.test)[names(df.test) =="variable"] <- "TSM.FACT"

    #Create the plotlist##      
        myplots <- list()

    #Loop for plots##
        for(i in 1:(length(unique(df.test$names)))){ 
          p1 <- eval(substitute(
            ggplot(data=df.test[df.test$names == levels(df.test$names)[i],], aes(x=Length2, group=TSM.FACT, colour = TSM.FACT, fill=TSM.FACT)) +
              geom_histogram(aes( y = stat(width*density)), position = "dodge", binwidth = 50, alpha =0.4, show.legend=T)+
              ggtitle(paste0(levels(df.test$names)[i]))+
              geom_density_line(stat="density", aes(y=(..count..)/sum(..count..)*50), alpha=0.3, size=0.5, show.legend=F) +
              geom_vline(data=ddply(df.test[df.test$names == levels(df.test$names)[i],], ~ TSM.FACT,  numcolwise(mean)), mapping=aes(xintercept = Length2, group=TSM.FACT, colour=TSM.FACT),  linetype=2, size=1, show.legend=F) +
              scale_y_continuous(labels = percent_format()) +
              ylab("relative frequency") +
              scale_color_manual(values= c("#00B2EE",  "#1E90FF",  "#104E8B")) +
              scale_fill_manual(values= c("#00B2EE",  "#1E90FF", "#104E8B")) +
            theme_bw() + theme(
              plot.title = element_text(lineheight=0.5, hjust= 0.5, size=10),
              strip.text.y = element_text(hjust = 1, angle = 0),
              strip.text.x = element_text(size=10, vjust = 0.9), 
              strip.text=element_text(margin = margin(t=0.3,r=1,b=0.3,l=1), size=8, debug = F, vjust=0.2),
              strip.background = element_blank(),
              axis.text.x = element_text(size=8, angle=0, vjust=0.2, margin = margin(t=0.3,r=0.1,b=0.3,l=0.1)),
              axis.title.x=element_blank(),
              axis.title.y=element_blank(),
              axis.line.x=element_line(colour="black"),
              axis.line.y=element_line(colour="black"),
              panel.grid.minor = element_blank(), 
              panel.border=element_blank(),
              panel.background=element_blank(),  
              legend.position=(c(0.9,0.9)),
              legend.title = element_blank(), 
              legend.key = element_blank()),
            list(i = i)))
          print(i)
          print(p1)
          myplots[[i]] <- p1 
          plot(p1)
        }

#Join the plots
panelplot=grid.arrange(plotlist = myplots, grobs = myplots, shared.legend=T)

enter image description here

1 Ответ

0 голосов
/ 22 октября 2019

К сожалению, я не могу воспроизвести ваш пример. Я могу порекомендовать добавить столбец, содержащий сумму каждого столбца (назовем его «Бар»)

Требуемое дополнение к коду ggplot включает: 

geom_text(aes(label = Bar), position = position_stack(vjust = 1)) +

Высота текста вышебар можно отрегулировать с помощью vjust

...