Question

我在使用quantile在我正在编写的Shiny应用程序中获取用户指定变量的百分位数时遇到问题，并使用这些百分位数值来设置直方图输出的x限制。

我在相关问题的答案中尝试了一系列修正，但不断得到与反应性，变量类型和NA相关的不同错误（即使我有na.rm = TRUE）。

这个例子是从一个更复杂的应用程序中减少的，所以我无法显示我已经尝试过的所有内容，但我很确定问题是在quantile命令中发生的，并且可能有几个我的代码有不同的问题。我也很难通过group_by_获取百分位数，而且我不确定我是否已成功解释NSE。

这是我的服务器文件（这里有问题）：

  {if (!require("devtools"))
  install.packages("devtools")
  if (!require("ggplot2")) 
    install.packages("ggplot2")
  if (!require("dplyr")) 
    install.packages("dplyr")
  if (!require("lazyeval")) 
    install.packages("lazyeval") }


#load libraries 
library(ggplot2) 
library(dplyr) 
library(lazyeval) 
library(shiny) 

#now the server function
shinyServer(function(input, output) {

  #output for histogram 
  output$histplot <- renderPlot({
    #make objects based on percentiles to trim the histogram you see
    lower.cut <- quantile(as.numeric(input$hist), probs = as.numeric(input$bottom), na.rm = TRUE) #error: missing value where TRUE/FALSE needed
    upper.cut <- quantile(as.numeric(input$hist), probs = as.numeric(input$top), na.rm = TRUE)
    q <- ggplot(iris, aes_string(input$hist)) + geom_histogram(binwidth = (as.numeric(input$n_breaks))) + 
      coord_cartesian(xlim = c(lower.cut, upper.cut)) #I'm not sure whether or not this is working because I can't get past the quantile commands
    q    
  })

  #output for summary stats - among other things, I would like to output the 95th, 98th, and 99th percentile values of the selected histogram variable
  output$summary <- renderPrint({  
  w <- iris %>% group_by_(as.numeric(iris$Species)) %>% summarize_(p95 = quantile(as.numeric(input$hist)), .95, na.rm=TRUE) #error: missing values and NaN's not allowed if 'na.rm' is FALSE
  w 
  })
})

这是我的UI文件（这似乎工作正常）：

{if (!require("devtools"))
  install.packages("devtools")
  if (!require("ggplot2")) 
    install.packages("ggplot2")
  if (!require("dplyr")) 
    install.packages("dplyr")
  if (!require("lazyeval")) 
    install.packages("lazyeval")
  if (!require("psych")) 
    install.packages("psych") }

library(ggplot2) 
library(dplyr) 
library(lazyeval) 
library(shiny) 

shinyUI(fluidPage(                      
  sidebarPanel( #this stuff will be in the left side panel - all user inputs and notes        
    h3("Histogram"),
    selectInput('hist', 'Histogram Variable', names(iris)), 
    numericInput('n_breaks', 'Histogram Bin Width (type a number, based on axis scale)', value = 10),
    numericInput('bottom', 'Histogram lower limit percentile (type a number between 0 and 1, .01 = exclude bottom 1%)', value = 0), #default no bottom trim
    numericInput('top', 'Histogram upper limit percentile', value = 1), #default no top trim

    hr() 
  ),

  mainPanel( #this stuff will be in the main body of the page - make placeholder spots for the charts
    fluidRow(
      plotOutput("histplot"), #simple histogram

      h4("Summary Statistics for Histogram Variable (grouped by species)"),
      verbatimTextOutput("summary") #output summary statistics below
  ))
))

以下是我收到的错误：

Warning in quantile(as.numeric(input$hist), probs = as.numeric(input$bottom),  :
  NAs introduced by coercion
Warning in quantile(as.numeric(input$hist), probs = as.numeric(input$top),  :
  NAs introduced by coercion
Warning: Error in if: missing value where TRUE/FALSE needed
Stack trace (innermost first):
    68: output$histplot
     1: runApp
Warning: Truncating vector to length 1
Warning in quantile(as.numeric(input$hist)) :
  NAs introduced by coercion
Warning: Error in quantile.default: missing values and NaN's not allowed if 'na.rm' is FALSE
Stack trace (innermost first):
    89: quantile.default
    88: quantile
    87: as.lazy_dots
    86: lazyeval::all_dots
    85: summarise_.tbl_df
    84: summarize_
    83: function_list[[k]]
    82: withVisible
    81: freduce
    80: _fseq
    79: eval
    78: eval
    77: withVisible
    76: %>%
    75: renderPrint [R:\Project\EPAR\Working Files\326 - Changes in Yield and Crop Allocation\R analysis\toy_app/server.R#32]
    74: func
    73: eval
    72: eval
    71: withVisible
    70: evalVis
    69: utils::capture.output
    68: paste
    67: output$summary
     1: runApp

这是sessionInfo（）输出：

> sessionInfo()
R version 3.2.4 Revised (2016-03-16 r70336)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252    LC_MONETARY=English_United States.1252 LC_NUMERIC=C                          
[5] LC_TIME=English_United States.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] shinythemes_1.0.1 psych_1.6.6       lazyeval_0.2.0    dplyr_0.5.0       foreign_0.8-66    ggplot2_2.1.0     devtools_1.12.0   shiny_0.13.2     

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.5      magrittr_1.5     mnormt_1.5-4     munsell_0.4.3    colorspace_1.2-6 xtable_1.8-2     R6_2.1.2         stringr_1.0.0    plyr_1.8.4      
[10] tools_3.2.4      parallel_3.2.4   grid_3.2.4       gtable_0.2.0     DBI_0.4-1        withr_1.0.2      htmltools_0.3.5  digest_0.6.9     assertthat_0.1  
[19] tibble_1.0       reshape2_1.4.1   memoise_1.0.0    mime_0.4         labeling_0.3     stringi_1.1.1    scales_0.4.0     jsonlite_0.9.22  httpuv_1.3.3

Answer 1

names(iris)在您的ui脚本中生成一个字符向量。然后，您的服务器脚本会尝试通过将此字符向量转换为数字来查找分位数，该数字会将您的所有名称转换为NA。

闪亮：“分位数”在反应性环境中不起作用

1 个答案: