metavcov (version 2.1.5)

plotCI: Plot Confidence Intervals for a Meta-Analysis

Description

The function plot.CI generates confidence interval figures for effect sizes from each study and the estimated effect sizes across studies.

Usage

plotCI(y, v,
       name.y = NULL,
       name.study = NULL,
       y.all, y.all.se,
       hline = 0,
       up.bound = Inf, low.bound = -Inf,
       return.data = FALSE)

Author

Min Lu

Arguments

y

A \(N \times p\) matrix or data frame that stores effect sizes from \(N\) primary studies. Usually the output value ef produced by r.vcov for correlation coefficients or mix.vcov for other types of effect sizes.

v

A \(N\)-dimensional list of \(p \times p\) matrices that stores within-study (co)variance matrices of the estimated effect sizes for each one of the \(N\) studies. Usually the output value list.vcov produced by r.vcov for correlation coefficients or mix.vcov for other types of effect sizes.

name.y

A \(p\)-dimensional vector that stores names for the effect sizes in y. By default, it is the columname of y.

name.study

A \(N\)-dimensional vector that stores names for the primary in y. By default, it is the row name of y.

y.all

A \(p\)-dimensional vector that stores the estimated effect sizes across studies.

y.all.se

A \(p\)-dimensional vector that stores the standard errors for the estimated effect sizes in y.all.

hline

A \(p\)-dimensional vector that specifies the position of the dash line in the figures to compare the coefficients across studies. If its length is one instead of \(p\), this number will be adopted for all the \(p\) effect sizes.

up.bound

A \(p\)-dimensional vector that specifies the upper bound in the figures. If its length is one instead of \(p\), this number will be adopted for all the \(p\) effect sizes.

low.bound

A \(p\)-dimensional vector that specifies the lower bound in the figures. If its length is one instead of \(p\), this number will be adopted for all the \(p\) effect sizes.

return.data

Should the data for the confidence interval plots be returned?

Details

The difference between a forest plot and a confidence interval plot is that a forest plot requires a symbol on each confidence interval that is proportional to the weight for each study. Because the weighting mechanism in multivariate meta-analysis is too complex to be visualized, such a propositional symbol is omitted for multivariate meta-analysis.

References

Lu, M. (2023). Computing within-study covariances, data visualization, and missing data solutions for multivariate meta-analysis with metavcov. Frontiers in Psychology, 14:1185012.

Examples

######################################################
# Example: Craft2003 data
######################################################
data(Craft2003)
computvcov <- r.vcov(n = Craft2003$N,
                     corflat = subset(Craft2003, select = C1:C6),
                     method = "average")
y <- computvcov$ef
Slist <- computvcov$list.vcov
MMA_FE <- summary(metafixed(y = y, Slist = Slist))
obj <- MMA_FE
# pdf("CI.pdf", width = 4, height = 7)
plotCI(y = computvcov$ef, v = computvcov$list.vcov,
        name.y = NULL, name.study = Craft2003$ID,
        y.all = obj$coefficients[,1],
        y.all.se = obj$coefficients[,2])
# dev.off()

######################################################
# Substitute obj for Random-effect model
######################################################
# library(mvmeta)
# S <- computvcov$matrix.vcov
# MMA_RE <- summary(mvmeta(cbind(C1, C2, C3, C4, C5, C6),
#                            S = S, data = y, method = "reml"))
# obj <- MMA_RE