Third Party Software
benruk
Joined: 12/09/2019
Multiple Moderators
Hello,
some days ago I read the paper by Shaun Purcell (2002) "Variance Components Models for Gene-Environment Interaction in Twin Analysis" and I'm especially interested in modelling a three-way interaction (GxExE). So, for example, the genetic effect may be moderated by the socioeconomic status and by age, so that the additive genetic variance component is:
$$
(a+\beta_{Age}M_{Age}+\beta_{SES}M_{SES}+\beta_{AgexSES}M_{Age}M_{SES})^{2}
$$
- Read more about Multiple Moderators
- 3 comments
- Log in or register to post comments
vildan Bağcı
Joined: 12/04/2019
Complex data anaysis in REM doesn't result
Dear Mike & colleagues,
First of all, I am grateful for asking my question to this forum where gurus of meta-SEM are helping people like me.
I'm trying to make analysis in your metasem package. I have 41 studies and so 41 correlation matrices. and the correlation matrices consist of 20x20.
- Read more about Complex data anaysis in REM doesn't result
- 17 comments
- Log in or register to post comments
qingwen
Joined: 08/15/2019
How to add CIs in umx?
I made a common pathway model in umx. I add CIs through adding umxCI(model, run="yes"). However, the output are NAs and there is a warning that "3: The nonlinear constraints and bounds could not be satisfied. The problem may have no feasible solution." I wonder what it means and what I can do to deal with it.
- Read more about How to add CIs in umx?
- 5 comments
- Log in or register to post comments
aqzhou
Joined: 10/21/2019
OpenMx status1: 6
Hi All,
Thanks in advance for your help and time! It is much appreciated.
I have 11 datasets (15 variables: 3 latent/3 indicators each, 2 observed DVs, 4 control variables) and am interested in using TSSEM. However, I'm running into errors and am getting stuck at stage 1. I have uploaded my data, syntax, and output.
- Read more about OpenMx status1: 6
- 4 comments
- Log in or register to post comments
Guoquan
Joined: 09/30/2019
Problem with TSSEM
Dear Mike & colleagues,
I am currently working on a TSSEM paper. However, some errors occur when I run the codes in R software. Please see data, my model and R codes in the attachments.
First, when I run TSSEM1, OpenMx status 1 is 6 and Warning message appears: "In print.summary.meta(x) : OpenMx status1 is neither 0 or 1. You are advised to 'rerun' it again." What should I do to solve it?
- Read more about Problem with TSSEM
- 1 comment
- Log in or register to post comments
raklein
Joined: 09/07/2019
OpenMX status1: 5 -- Problematic in this case?
Hi all,
I'm attaching some data and code to demonstrate my problem, but long story short I'm wondering if there is a fix for the 'OpenMX status1: 5' problem I'm having, AND/OR whether it is OK to report these values as-is? (perhaps with a grain of salt?) We're hoping to submit this paper for publication shortly.
Pruealexis
Joined: 08/18/2019
Outliers/Influential cases and Publication Bias
Hi all,
I am currently using metaSEM to create a pooled correlation matrix from 20 studies in my meta-analysis.
I previously used metafor to do this, conducting a metaanalysis via Hunter-Schmidt estimation/method on uncorrected data (maModel), and identified any outliers/influential cases. Using this method, I also completed a leave one out analysis. Metafor also enabled creating Forest Plots and Funnel Plots (to check for publication bias).
nastjuscha
Joined: 06/28/2016
TSSEM2 with different combinations of variables
Dear Mike,
I have a question about possibilities for model specification at stage 2 of MetaSEM.
So I have a dataset with 10 variables - at the end of stage 1 I get a pooled correlation matrix. In stage 2 I would like to test some models with different constellations of the variables, e.g. with variables V1, V2, V3, V4, V5, V6, V7 and V8 in the first model and with V1, V2, V3, V4, V5, V9, V10 in the second. Is that possible? And if it is possible, how should I specify this in the syntax?
Thanks for your help!
forscher
Joined: 03/05/2013
Best way to test constraints on parameters in metaSEM models
I have metaSEM model (fit using metaSEM::meta) for which I wish to test constraints on the estimated intercepts. I've been using car::linearHypothesis because it's extremely flexible -- it works on any model that provides a variance-covariance matrix (via vcov()) and a vector of estimated coefficients (via coef()). Given a metaSEM model, car::linearHypothesis seems to allow the testing of joint hypotheses (ex, the hypothesis that Intercept1 = 0 AND Intercept2 = 0) and complex hypotheses (ex, the hypothesis that Intercept1 = Intercept2).
Xin Lin
Joined: 06/10/2019
MetaSEM Error Message-Missing value where TRUE/FALSE needed
Hi Dr. Cheung,
I am trying to use my dataset to mimic your code “fixed1 <- tssem1(mat,n, method = "FEM")”
I got the error message “Error in if (!all(isPD)) warning(paste("Group ", (1:no.groups)[!isPD], :
missing value where TRUE/FALSE needed”
I have attached the first five correlation matrices and the sample size (n) for each study.