如何在 Shiny 中创建可点击的直方图？

Question

我想在 shiny 中创建一个可点击的直方图，但我不知道是否可行。

几个月前我看到了一个可点击的火山图，它为您提供了点击内容的 table。

来源：https://2-bitbio.com/2017/12/clickable-volcano-plots-in-shiny.html

我发现关于创建可点击直方图的最接近 post 的是这个

但是，我不想获取坐标。我想要数据框的行名。

有了这个数据框，我可以在每次点击直方图中的条形时获取行名吗？

mtcars <- mtcars %>% 
  select("hp")
mtcars <- as.matrix(mtcars)

闪亮的一个例子（但不可点击）：

library(shiny)
library(ggplot2)
library(scales)
library(dplyr)

ui <- fluidPage(
  
  titlePanel("Histogram"),
  
  sidebarLayout(
    sidebarPanel(
    ),
    
    mainPanel(
      plotOutput("hist"),
    )
  )
)

mtcars <- mtcars %>% 
  select("hp")
mtcars <- as.matrix(mtcars)

server <- function(input, output) {
  
  output$hist <- renderPlot({
    
    pp <- qplot(mtcars, geom = "histogram", bins = 10, xlab="values", 
                ylab="Frequency", main="Histogram",
                fill=I("red"), col=I("black"), alpha=I(0.4))
    
   pp + scale_x_continuous(breaks=pretty(mtcars, n=10))
  })
  
  
}

shinyApp(ui = ui, server = server)

有人知道怎么做吗？

非常感谢！

此致

Answer 1

这是一个很好的问题，而具有挑战性的是 qplot/ggplot 图表是静态图像。下面的 app.r 是我将如何做的一个例子。我很乐意看到其他方法。

本质上：

创建一个数字序列，既可以用作直方图中的间隔，也可以用作数据框中的间隔。我将这些基于用户输入，但您可以对它们进行硬编码。
根据值所在的区间为数据框中的每一行分配一个“bin”值。
记录用户点击事件的 x 坐标，并根据同一组间隔为其分配一个“bin”值。
子集化您的数据框并仅保留数据的“bin”值与用户点击事件的 x 坐标的“bin”值相匹配的记录。

否则，如果您愿意走 d3 路线，您可以探索 something like this 发布的 R Views。

#Load libraries ----------------------------------------------------
library(shiny)
library(ggplot2)
library(scales)
library(dplyr)


# Prepare data -----------------------------------------------------
df <- mtcars
df <- cbind(model = rownames(df), data.frame(df, row.names = NULL)) # setting the rownames as the first column
dm <- df$hp %>% as.matrix()


# UI function ------------------------------------------------------
ui <- fluidPage(

  titlePanel("Histogram"),

  sidebarLayout(
    sidebarPanel(

      tags$h5("I added the below text output only to demonstrate shiny's way for tracking user interaction on static plots. You can click, double-click, or click & drag (i.e. brushing). These functions are AWESOME when exploring scatterplots."),

      tags$h3("Chart click and brushing"),

      verbatimTextOutput("info"),

      tags$h5("Now I'm applying the below UI inputs to the `vec` and `breaks` arguments in `findInterval()` and `qplot()` respectively; I'm using `findInterval()` to bin the values in the dataframe AND to bin the x-value of the user's click event input on the chart. Then we can return the dataframe rows with the same bin values as the x-value of the click input."),

      sliderInput("seq_from_to"
                  , label = h3("Sequence 'From' and 'To'")
                  , min = 0
                  , max = 500
                  , value = c(50, 350)
                  ),

      sliderInput("seq_by"
                  , label = h3("Sequence 'By'")
                  , min = 25
                  , max = 200
                  , value = 50
                  , step = 5)

    ),

    mainPanel(

      plotOutput("hist",
                 click = "plot_click",
                 dblclick = "plot_dblclick",
                 hover = "plot_hover",
                 brush = "plot_brush"),

      dataTableOutput("table")

    )
  )
)


# Server function --------------------------------------------------
server <- function(input, output) {

  # Render Histogram Plot
  output$hist <- renderPlot({

    # Using the same `qplot` function but inserting the user inputs to set the breaks values in the plot
    pp <- qplot(dm
                , geom = "histogram"
                , breaks = seq(from = input$seq_from_to[1], to = input$seq_from_to[2], by = input$seq_by)
                , xlab = "values"
                , ylab = "Frequency"
                , main = "Histogram"
                , fill = I("red")
                , col = I("black")
                , alpha = I(0.4)
                )

    # Also using the user inputs to set the breaks values for the x-axis
    pp + scale_x_continuous(breaks = seq(from = input$seq_from_to[1], to = input$seq_from_to[2], by = input$seq_by))
  })

  # This is purely explanatory to help show how shiny can read user interaction on qplot/ggplot objects
  # It's taken from the Shiny docs here: https://shiny.rstudio.com/articles/plot-interaction.html
  output$info <- renderText({

    # Retain the x and y coords of the user click event data
    xy_str <- function(e) {
      if(is.null(e)) return("NULL\n")
      paste0("x=", round(e$x, 1), " y=", round(e$y, 1), "\n")
    }

    # Retain the x and y range coords of click & drag (brush) data
    xy_range_str <- function(e) {
      if(is.null(e)) return("NULL\n")
      paste0("xmin=", round(e$xmin, 1), " xmax=", round(e$xmax, 1),
             " ymin=", round(e$ymin, 1), " ymax=", round(e$ymax, 1))
    }

    # Paste this together so we can read it in the UI function for demo purposes
    paste0(
      "click: ", xy_str(input$plot_click),
      "dblclick: ", xy_str(input$plot_dblclick),
      "hover: ", xy_str(input$plot_hover),
      "brush: ", xy_range_str(input$plot_brush)
    )
  })

  # Back to the story. Set a listener to trigger when one of the following is updated:
  toListen <- reactive({list(
    input$plot_click    # user clicks on the plot
    , input$seq_from_to # user updates the range slider
    , input$seq_by      # user updates the number input
    )
  })

  # When one of those events are triggered, update the datatable output
  observeEvent(toListen(), {

    # Save the user click event data
    click_data <- input$plot_click
    print(click_data) # during your app preview, you can watch the R Console to see what click data is accessible

    # Assign bin values to each row using the intervals that are set by the user input
    df$bin <- findInterval(dm, vec = seq(from = input$seq_from_to[1], to = input$seq_from_to[2], by = input$seq_by))

    # Similarly assign a bin value to the click event based on what interval the x values falls within
    click_data$x_bin <- findInterval(click_data$x, vec = seq(from = input$seq_from_to[1], to = input$seq_from_to[2], by = input$seq_by))

    # Lastly, subset the df to only those records within the same interval as the click event x-value
    df_results <- subset(df, bin == click_data$x_bin)

    # Select what values to view in the table
    df_results <- df_results %>% select(model, hp)

    # And push these back out to the UI
    output$table <- renderDataTable(df_results,
                                     options = list(
                                       pageLength = 5
                                     )
    )

  })


}

shinyApp(ui = ui, server = server)

Answer 2

嗯，有人回答了。由于我花时间将它放在一起，这里是另一个可能的解决方案。

library(shiny)
library(ggplot2)
library(scales)
library(dplyr)
library(DescTools)                             # added for Closest()

ui <- fluidPage(                    
  
  titlePanel("Histogram"),
  sidebarLayout(
    sidebarPanel(
    ),
    
    mainPanel(
      plotOutput("hist", click = 'plot_click'),  # added plot_click
      verbatimTextOutput("x_value"),             # added queues for interactivity
      verbatimTextOutput("selected_rows")        # added table for bin values
    )
  )
)

# this can be a dataframe or matrix for qplot or ggplot 
           # (not sure if there was another reason you had this code?)
# mtcars <- mtcars %>% 
#   select("hp")                      # if you only want hp
# mtcars <- as.matrix(mtcars)         # I suggest making row names a column
                                      # to keep 2 columns 

pp <- ggplot(mtcars) +
  geom_histogram(aes(x = hp),
                 bins = 10,
                 fill = "red",
                 color = "black",
                 alpha = .4) +
  labs(x = "values",
       y = "Frequency",
       title = "Histogram")

# extract data from plot to find where each value falls within the histogram bins
        # I kept the pkg name, function in more than one library
bd <- ggplot_build(ggplot2::last_plot())$data[[1]]   

# add the assigned bin number to the mtcars frame; used for filtering matches
mtcars$bins <- lapply(mtcars$hp,
                      function(y) {
                        which(bd$x == Closest(bd$x, y))
                      }) %>% unlist()

server <- function(input, output) { 

  output$hist <- renderPlot({
    # moved the plot outside of server, so that global variables could be created
    
    # pp <- qplot(mtcars[,"hp"], geom = "histogram", bins = 10, xlab="values", 
    #             ylab = "Frequency", main = "Histogram",
    #             fill = I("red"), col = I("black"), alpha = I(0.4))
    # scale_x_continuous(breaks=pretty(mtcars, n=10)) # can't use this

    pp
  })
  # # Print the name of the x value                 # added all that's below with server()
  output$x_value <- renderPrint({
    if (is.null(input$plot_click$x)) return()

    # find the closest bin center to show where the user clicked on the histogram
    cBin <- which(bd$x == Closest(bd$x, input$plot_click$x))
    paste0("You selected bin ", cBin)   # print out selected value based on bin center
  })
  # Print the rows of the data frame which match the x value
  output$selected_rows <- renderPrint({
    if (is.null(input$plot_click$x)) return()
    
    # find the closest bin center to show where the user clicked on the histogram
    cBin <- which(bd$x == Closest(bd$x, input$plot_click$x))
    mtcars %>% filter(bins == cBin)
    # mtcars
  })
}

shinyApp(ui = ui, server = server)

Answer 3

以防万一有人在这个 post 中寻找 方法来包含 brushedPoints... 受到这个 post 的启发，我找到方法了！

代码：

#Load libraries ----------------------------------------------------
library(shiny)
library(ggplot2)
library(scales)
library(dplyr)


# Prepare data -----------------------------------------------------
df <- mtcars
df <- cbind(model = rownames(df), data.frame(df, row.names = NULL)) # setting the rownames as the first column
breaks_data = pretty(mtcars$hp, n=10)
my_breaks = seq(min(breaks_data), to=max(breaks_data), by=30)


# UI function ------------------------------------------------------
ui <- fluidPage(
  
  titlePanel("Histogram"),
  
  sidebarLayout(
    sidebarPanel(
      actionButton("draw_plot", "Draw the plot")
    ),
    
    mainPanel(
      
      plotOutput("hist",
                 brush = brushOpts("plot_brush", resetOnNew = T, direction = "x")),
      
      dataTableOutput("table"),
    )
  )
)


# Server function --------------------------------------------------
server <- function(input, output) {
  

  observeEvent(input$plot_brush, {
   info_plot <- brushedPoints(df, input$plot_brush)
   output$table <- renderDataTable(info_plot)
   
  })
  
  # If the user didn't choose to see the plot, it won't appear.
  output$hist <- renderPlot({
    df %>% ggplot(aes(hp)) +
      geom_histogram(alpha=I(0.4), col = I("black"), fill = I("red"), bins=10) +
      labs(x = "values",
           y = "Frequency",
           title = "Histogram") +
      scale_x_continuous(breaks = my_breaks)
    
    
  })
  
}

shinyApp(ui = ui, server = server)

如何在 Shiny 中创建可点击的直方图？

How to create a clickable histogram in Shiny?

r

histogram

clickable

shiny