This method function performs a Cochran Q test of (residual) heterogeneity on fitted meta-analytical models represented in objects of class "mixmeta".
# S3 method for mixmeta
qtest(object, …)# S3 method for qtest.mixmeta
print(x, digits=3, …)
objects of classes "mixmeta" and "qtest.mixmeta", respectively.
an integer specifying the number of digits to which printed results must be rounded.
further arguments passed to or from other methods.
A list object of class "qtest.mixmeta" with the following components:
the vector of test statistics for overall and outcome-specific tests, distributed under the null hypothesis as a Chi-square with degrees of freedom df.
the vector of degrees of freedom of the null distribution for overall and outcome-specific tests. For the overall test, equal to the number of observations used for estimation minus the number of coefficients in the fixed part of the model. For outcome-specific test, equal to number of observed values minus the number of coefficients.
the vector of p-values for overall and outcome-specific tests.
logical switch indicating if a meta-regression model is assessed, meaning that the tested heterogeneity is residual.
dimensionality of the overall test, that is the number of outcome parameters in the model.
As usual, the print method function for class "qtest.mixmeta" does not return any value.
The test assesses the null hypothesis that the variability in the distribution of the outcomes is explained only in terms of estimation error in each unit, measured by the within-unit (co)variance matrices stored in the component S of mixmeta objects. This is equal to test the hypothesis that the random-effects (co)variance matrix (or all matrices in multilevel models) is a zero matrix, and there is no random deviation in unit-specific estimates. For multivariate models, tests for single outcome parameters, comparable to estimates from multiple univariate meta-analysis, are also reported. This test reduces to the standard Q test in univariate single-level models.
The function compute the statistics by actually fitting the related fixed-effects model, re-evaluating the call of the model with method changed to "fixed".
Sera F, Armstrong B, Blangiardo M, Gasparrini A (2019). An extended mixed-effects framework for meta-analysis.Statistics in Medicine. 2019;38(29):5429-5444. [Freely available here].
Cochran WG (1950). The comparison of percentages in matched samples". Biometrika. 37(3/4):256--266.
See qtest for the generic method function. See mixmeta-package for an overview of the package and modelling framework.
# NOT RUN {
# RUN THE MODEL
model <- mixmeta(cbind(PD,AL) ~ 1, S=berkey98[5:7], data=berkey98)
# MULTIVARIATE COCHRAN Q TEST FOR HETEROGENEITY
test <- qtest(model)
print(test, digits=2)
unclass(test)
# }
Run the code above in your browser using DataLab