包装 R 的绘图函数（或 ggplot2）以防止绘制大型数据集

Question

我不想问如何绘制大数据集，而是想包装plot，以便生成大量绘图的代码在绘制大对象时不会受到打击。如何以非常简单的方式包装plot，以便保留其所有功能，但首先要测试以确定传递的对象是否太大？

此代码适用于对plot 的非常普通的调用，但它缺少与plot 相同的通用性（见下文）。

myPlot <- function(x, ...){
    isBad <- any( (length(x) > 10^6) || (object.size(x) > 8*10^6) || (nrow(x) > 10^6) )
    if(is.na(isBad)){isBad = FALSE}
    if(isBad){
        stop("No plots for you!")
    }
    return(plot(x, ...))
}

x = rnorm(1000)
x = rnorm(10^6 + 1)

myPlot(x)

一个失败的例子：

x = rnorm(1000)
y = rnorm(1000)
plot(y ~ x)
myPlot(y ~ x)

是否有一些简单的方法来包装plot 以启用对要绘制的数据的检查，同时仍然传递所有参数？如果没有，那么ggplot2 怎么样？我是一个机会均等的非策划者。 （在数据集比较大的情况下，我会使用hexbin、下采样、密度图等，但这不是这里的重点。）

---

注1：在测试想法时，我建议测试大小> 100（或设置变量，例如myThreshold <- 1000），而不是与> 1M的大小相比 - 否则会遇到很多痛苦绘图。 :)

Answer 1

您遇到的问题是，按照当前的编码，myplot() 假定 x 是一个数据对象，但是您尝试向它传递一个公式。 R 的 plot() 通过方法实现这一点 - 当 x 是一个公式时，plot.formula() 方法被调度到而不是基本的 plot.default() 方法。

你也需要这样做：

myplot <- function(x, ...)
    UseMethod("myplot")

myplot.default <- function(x, ....) {
    isBad <- any((length(x) > 10^6) || (object.size(x) > 8*10^6) || 
                    (nrow(x) > 10^6))
    if(is.na(isBad)){isBad = FALSE}
    if(isBad){
        stop("No plots for you!")
    }
    invisible(plot(x, ...))
}

myplot.formula <- function(x, ...) {
    ## code here to process the formula into a data object for plotting
    ....
    myplot.default(processed_x, ...)
}

您可以从plot.formula() 窃取代码以用于将x 处理为对象所需的代码。或者，您可以在standard non-standard evaluation rules (PDF) 之后自行滚动。