在R中绘制引导结果的图形

Question

假设我有以下代码：

library(boot)
samplemean <- function(x, d) {
  return(mean(x[d]))
}
results_qsec <- boot(data=mtcars$qsec, statistic = samplemean, R=1000)
results_wt <- boot(data=mtcars$wt, statistic = samplemean, R=1000)
plot(results_qsec)
plot(results_wt)

plot(results_wt)我得到以下内容：

我可以编辑绘制的内容吗？例如，我想摆脱右边的图形，将标题从“直方图t”更改为“香蕉直方图”，并在同一图表上显示results_qsec和results_wt的直方图。 可以这样做吗？我查看了启动文档，但找不到任何有用的信息。 谢谢 贝

Answer 1

这个怎么样？当然，您可以更改纸槽宽度等，使其看起来更像原始纸张。

wt_t <- plot(results_wt)$t
sec_t <- plot(results_qsec)$t

par(mfrow = c(1, 2))
hist(wt_t, main = "banana")
hist(sec_t, main = "apple")

Answer 2

谢谢，看起来不错。 我找到了另一种方法，使用ggplot：

library(ggplot2)

boot_qsec <- as.data.frame(results_qsec$t)
boot_wt <- as.data.frame(results_wt$t)

ggplot() + geom_histogram(data=boot_qsec,aes(V1)) + 
      geom_histogram(data=boot_wt ,aes(V1))