Nothing
### library(metafor); library(testthat); Sys.setenv(NOT_CRAN="true")
### see also: https://www.metafor-project.org/doku.php/plots:caterpillar_plot
source("settings.r")
context("Checking plots example: caterpillar plot")
test_that("plot can be drawn.", {
expect_equivalent(TRUE, TRUE) # avoid 'Empty test' message
skip_on_cran()
### simulate some data
set.seed(5132)
k <- 250
vi <- rchisq(k, df=1) * .03
yi <- rnorm(k, rnorm(k, 0.5, 0.4), sqrt(vi))
### fit RE model
res <- rma(yi, vi)
png("images/test_plots_caterpillar_plot_test.png", res=200, width=1800, height=1500, type="cairo")
### decrease margins so the full space is used
par(mar=c(5,1,1,1))
### create plot
forest(yi, vi,
xlim=c(-2.5,3.5), ### adjust horizontal plot region limits
order=yi, ### order by size of yi
slab=NA, annotate=FALSE, ### remove study labels and annotations
efac=0, ### remove vertical bars at end of CIs
pch=19, ### changing point symbol to filled circle
col="gray40", ### change color of points/CIs
psize=2, ### increase point size
cex.lab=1, cex.axis=1, ### increase size of x-axis title/labels
lty=c("solid","blank")) ### remove horizontal line at top of plot
### draw points one more time to make them easier to see
points(sort(yi), k:1, pch=19, cex=0.5)
### add summary polygon at bottom and text
addpoly(res, mlab="", cex=1)
text(-2, -2, "RE Model", pos=4, offset=0, cex=1)
dev.off()
expect_true(.vistest("images/test_plots_caterpillar_plot_test.png", "images/test_plots_caterpillar_plot.png"))
})
rm(list=ls())
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.