Multivariate Metaanalysis (II)
Reference: Chen, Y, Hong, C and Riley, R. (2015) An alternative pseudolikelihood method for multivariate randomeffects metaanalysis, Statistics in Medicine 34 (3): 361380.
Introduction
Recently, multivariate randomeffects metaanalysis models have received a great deal of attention, despite its greater complexity compared to those for univariate metaanalysis. One of its advantages is the ability to account for withinstudy and betweenstudy correlations. However, the standard inference procedures, such as the maximum likelihood or maximum restricted likelihood inference, require the withinstudy correlations, which are usually unavailable. In addition, the standard inference procedures suffer from the estimation issues such as nonconvergent iterations or a singular estimated covariance matrix.
â€‹
In this paper, we propose a pseudolikelihood method to overcome the aforementioned problem when outcome data are complete or missing completely at random (MCAR), and interest is in functions of the pooled estimates. The pseudolikelihood method does not require withinstudy correlations and is not prone to nonconvergence or a singular covariance matrix problem. In addition, it can properly estimate the covariance between pooled estimates for different outcomes, which enables valid inference on functions of pooled estimates, and can be applied to metaanalysis where some studies have outcomes MCAR.
R Packages and Analysis Functions
In this paper, we analyze the example data using the proposed "Composite Likelihood" method in R.
The proposed method is implemented in the ’mmeta’ function in our R software package ‘xmeta’, and can be installed from CRAN (http://cran.rproject.org/package=xmeta/), the official R package archive. Make sure that you load them before trying to run the examples on this page. If you do not have the package installed, run:
Structure:
Before running the SPLUS/R program code to fit the pseudoREML method, you need to have the following two package in your R package repository:
Then, the following function can be directly used to analyze metaanalysis data:
Input:

mydat: data of n rows and 4 columns, where n is the number of studies and 4 columns are Yi1 (effect size of outcome 1), Si1^2(variance of Yi1), Yi2(effect size of outcome 2), Si2^2(variance of Yi2)
Output:

estim: point estimates of (beta1, tau1^2, beta2, tau2^2)

cov.estim: covariance matrix of the point estimates
Implementation of Working Example
In this working example, we explore the comparison between overall survival and diseasefree survival for prostate cancer Loghazard ratio estimates comparing combined therapy using Goserelin acetate with radiotherapy with respect to the overall survival (Yi1) and that with respect to diseasefree survival (Yi2).
Dataset:
To obtain the point estimates and covariance estimate, you need to run:
The result is:
The model based standard errors of beta1, beta2 and (beta1beta2):
And you will get: