Refer to rq by name

Author: hadley

R/stat-quantile.r

@@ -16,14 +16,16 @@ StatQuantile <- proto(Stat, {
   default_aes <- function(.) aes(group = ..quantile..)
   required_aes <- c("x", "y")
 
-  calculate <- function(., data, scales, quantiles=c(0.25, 0.5, 0.75), formula=y ~ x, xseq = NULL, method=rq, na.rm = FALSE, ...) {
+  calculate <- function(., data, scales, quantiles=c(0.25, 0.5, 0.75), formula=y ~ x, xseq = NULL, method="rq", na.rm = FALSE, ...) {
     try_require("quantreg")
     if (is.null(data$weight)) data$weight <- 1 
 
     if (is.null(xseq)) xseq <- seq(min(data$x, na.rm=TRUE), max(data$x, na.rm=TRUE), length=100)
 
     data <- as.data.frame(data)
     data <- remove_missing(data, na.rm, c("x", "y"), name = "stat_quantile")
+    
+    method <- match.fun(method)
     model <- method(formula, data=data, tau=quantiles, weight=weight, ...)
 
     yhats <- stats::predict(model, data.frame(x=xseq), type="matrix")