You are here

Errors when including covariates in metaSEM

4 posts / 0 new
Last post
forscher's picture
Offline
Joined: 03/05/2013 - 09:31
Errors when including covariates in metaSEM
AttachmentSize
File Data for metaSEM forum.csv134.09 KB
Binary Data Code for metaSEM forum.R2.06 KB

I have been attempting to add covariates to a network meta-analytic model that I'm fitting in metaSEM. As a brief bit of background, the network meta-analytic model is designed to model comparisons between a reference group and a set of other groups. Each comparison between the reference and other groups is modeled as a separate outcome. The S matrix in these models is often quite sparse because it is often the case that only ~half the studies contain multi-group designs, so one often needs to place constraints (e.g., with RE.constraints in metaSEM) on the between-studies covariance matrix for the model to be identifiable.

I am able to fit the model with a covariate using the mvmeta package in R, and I have verified that there is variance in my covariate for each of the outcomes (i.e., comparisons between reference and other groups) in the meta-analysis. However, when I attempt to fit the model with the single covariate in metaSEM, I receive the following error:

Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
0 (non-NA) cases

Does anybody have any idea what's happening here? For reference, I have tried fitting the covariate model on subsets of the data (i.e., using only one of the 11 outcomes from the meta-analysis) without generating these errors. My data and a script are attached.

Mike Cheung's picture
Offline
Joined: 10/08/2009 - 22:37
Hi Patrick, The meta()

Hi Patrick,

The meta() function uses the regression coefficients of a regression analysis as starting values. In your case, there is no enough data for the regression analysis after excluding the missing data with listwise deletion. You may skip this step by manually providing starting values to meta(). Please see the attached analysis.

Cheers,
Mike

forscher's picture
Offline
Joined: 03/05/2013 - 09:31
Thanks for your response! I

Thanks for your response!

I noticed that the "Explained variances" section of the output has a bunch of NAs. I assume that this section of the output relies on these same functions that use listwise deletion?

Mike Cheung's picture
Offline
Joined: 10/08/2009 - 22:37
No. The so-called "explained

No. The so-called "explained variances" are based on comparing the heterogeneity variances with and without the predictors. Since there are constraints on your variance components, I don't think that the meta() function is smart enough to calculate the correct "explained variances."