testOutliers: Tests for the presence of outliers.

testOutliersR Documentation

Tests for the presence of outliers.


For the t-distribution models this is a test that the degrees of freedom is infinitite, or equivalently that it's inverse is zero. For the mixture-normal model it is a test that the proportion of outliers is zero. As both tests involve a parameter on the boundary of the parameter space, asymptotic theory does not apply, so a parametric bootstrap is performed to determine the empirical distribution of the test statistic under the null hypothesis. The observed likelihood ratio statistic is then compared to this distribution to determine the p value.


## S3 method for class 'metaplus'
testOutliers(object, R = 999, cores = max(detectCores()%/%2, 1))



The meta-analysis for which the presence of outliers is to be tested.


Number of simulations (parametric bootstraps) used in testing the hypothesis. Initially, it may be useful to set this to a smaller value, to allow faster execution time.


Number of cores used to simultaneously perform simulations.



p value obtained from parametric bootstrap


Observed value of the likelihood ratio test statistic


Simulated values of the test statistic under the null hypothesis


Running the default number of bootstrap samples may take considerable time, of the order of hours. Use of the parallel options will improve these times. For a preliminary guide to whether there are outliers, this is indicated by a reduction in AIC or BIC with the robust model, and a change in the estimated mean effect.


Ken Beath <ken@kjbeath.id.au>


cdp3 <- metaplus(yi, sei , plotci = TRUE, slab = study, random = "mixture", cores = 1, data = cdp)
summary(testOutliers(cdp3, cores = 1))

metaplus documentation built on July 4, 2024, 5:08 p.m.