Question

上下文：我有一个包含50多个特征的数据集，我想为每个特征生成一个箱线图，直方图和摘要统计信息，以便进行演示。这使得150多个情节。我上面提到的代码是这样的：

library(ggplot2)
library(dplyr)
library(ggpubr)
library(ggthemes)
library(Rmisc)
library(gridExtra)    

myplots <- list()  # new empty list

for (i in seq(2,5,3)){
  local({
    i <- i
    p1 <- ggplot(data=dataset,aes(x=dataset[ ,i], colour=label))+ 
      geom_histogram(alpha=.01, position="identity",bins = 33, fill = "white") +
      xlab(colnames(dataset)[ i]) + scale_y_log10()  + theme_few()
    p2<- ggplot(data=dataset, aes( x=label, y=dataset[ ,i], colour=label)) +
      geom_boxplot()+ylab(colnames(dataset)[ i]) +theme_few()
    p3<- summary(dataset[ ,i])
    print(i)
    print(p1)
    print(p2)
    print(p3)
    myplots[[i]] <<- p1  # histogram
    myplots[[i+1]] <<- p2 # boxplot
    myplots[[i+2]] <<- p3 # summary
  })
}

myplots[[2]]
length(myplots)

n <- length(myplots)
nCol <- floor(sqrt(n))
do.call("grid.arrange", c(myplots, ncol=nCol)) # PROBLEM: cant print summary as  grob

我创建了一个绘图列表，每3个元素代表每个特征的直方图，箱图和摘要的结果。我遍历了50多个功能中的每一个，将每个结果附加到我的列表中（不是我知道的最好的方法）。当我尝试通过网格排列打印列表时，我遇到了以下问题：

Error in gList(list(grobs = list(list(x = 0.5, y = 0.5, width = 1, height = 1,  : 
  only 'grobs' allowed in "gList"

可以理解，因为摘要函数不会产生图形对象。关于如何克服这一挫折，除了不包括摘要统计数据之外的任何想法？

Answer 1

在结合了几个建议后，我设法弄清楚如何在循环浏览数据集的不同功能后，将每个要素的摘要统计数据绘制为grob对象。

library(skimr)
library(GridExtra)
library(ggplot2)
library(dplyr)
mysumplots <- list() # new empty list

for (i in seq(2,ncol(dataset))){
  local({
    i <-         
    sampletable <- data.frame(skim((dataset[ ,i]))) #creates a skim data frame 
    summarystats<-select(sampletable, stat, formatted) #select relevant df columns
    summarystats<-slice(summarystats , 4:10) #select relevant stats
    p3<-tableGrob(summarystats, rows=NULL) #converts df into a tableGrob

    mysumplots[[i]] <<- p3 # summary #appends the grob of to a list of summary table grobs
  })
}

do.call("grid.arrange", c(mysumplots, ncol=3)) # use grid arrange to plot all my grobs

这样做是为每个列（特征）创建一个skim数据帧，然后我选择了相关的统计数据，并将该grob分配给变量p3，然后将变量p3迭代地附加到每个特征的tablegrobs列表中。然后我使用gridarrange来打印所有tableGrobs！

如何将plot（）的结果显示为plot / grob？

1 个答案: