ctmaBiG: ctmaBiG

In CoTiMA: Continuous Time Meta-Analysis ('CoTiMA')

ctmaBiG

Description

Analysis of publication bias and generalizability. The function takes a CoTiMA fit object (created with ctmaInit) and estimates fixed and random effects of single drift coefficients, heterogeneity (Q, I square, H square, tau square), PET-PEESE corrections, Egger's tests, and z-curve analysis yielding expected replication and detection rates (ERR, EDR).

Usage

ctmaBiG(
  ctmaInitFit = NULL,
  activeDirectory = NULL,
  PETPEESEalpha = 0.1,
  activateRPB = FALSE,
  digits = 4,
  zcurve = FALSE,
  undoTimeScaling = TRUE,
  dt = NULL
)

Arguments

ctmaInitFit	fit object created with ctmaInit containing the fitted ctsem model of each primary study
activeDirectory	the directory where to save results (if not specified, it is taken from ctmaInitFit)
PETPEESEalpha	probability level (condition) below which to switch from PET to PEESE (cf. Stanley, 2017, p. 582, below Eq. 2; default p = .10)
activateRPB	if TRUE, messages (warning, finished) could be send to smart phone (default = FALSE)
digits	rounding (default = 4)
zcurve	performs z-curve analysis. Could fail if too few studies (e.g. around 10) are supplied. default=FALSE
undoTimeScaling	if TRUE, the original time scale is used (timeScale argument possibly used in ctmaInit is undone )
dt	A scalar indicating a time interval across which discrete time effects should be estimated and then used for ctmaBiG.

Value

ctmaBiG returns a list containing some arguments supplied, the results of analyses of publication bias and generalizability, model type, and the type of plot that could be performed with the returned object. The arguments in the returned object are activeDirectory, and coresToUse. Further arguments, which are just copied from the init-fit object supplied, are, n.studies, n.latent, studyList, statisticsList, modelResults (all parameter estimates and their standard error), and parameter names. All new results are returned as the list element "summary", which is printed if the summary function is applied to the returned object. The summary list element comprises a title (model='Analysis of Publication Bias & Generalizability') and "estimates", which is another list comprising "Fixed Effects of Drift Coefficients", "Heterogeneity", "Random Effects of Drift Coefficients", "PET-PEESE corrections", "Egger's tests" (constant of the WLS regression of drift coefficients on their standard errors (SE) with 1/SE^2 as weights), "Egger's tests Alt. Version" (constant of the OLS regression of the standard normal deviates of the drift coefficients on their precision), and "Z-Curve 2.0 Results". Plot type is plot.type=c("funnel", "forest") and model.type="BiG".

Examples

## Not run: 
# perform analyses of publication bias and generalizability
CoTiMAInitFit_D_BO$activeDirectory <- "/Users/tmp/" # adapt!
CoTiMABiG_D_BO <- ctmaBiG(ctmaInitFit=CoTiMAInitFit_D_BO, zcurve=FALSE)

## End(Not run)

# display results
summary(CoTiMABiG_D_BO)

## Not run: 
# get funnel & forest plots
CoTiMABiG_D_BO$activeDirectory <- "/Users/tmp/" # adapt!
plot(CoTiMABiG_D_BO)

## End(Not run)

CoTiMA documentation built on May 29, 2024, 11:39 a.m.