Behavioral Genetics Models

Posted on
No user picture. diana Joined: 04/17/2019

Why mxTryHardOrdinal()dosen't report start values?

Hi, everyone!
When i useSatFit <- mxTryHardOrdinal(SatModel, intervals=TRUE,extraTries = 15,finetuneGradient=FALSE,bestInitsOutput=TRUE) ,the result is Solution found! Final fit=22505.957 (started at 22819.562) (16 attempt(s): 16 valid, 0 errors). I don't know why it dosen't report start values from the best fit.

Posted on
No user picture. diana Joined: 04/17/2019

saturated univariate SEM model of a binary variable

Dear all,
I have some problems of fitting saturated univariate SEM model of a binary variable "whtr_g". I would appreciate it if anyone could help me.
1. I always get warning message when fitting saturated SEM model(see warning message.png), whether using SLSQP, NPSOL or CSOLNP optimizer. What’ more, the result of three optimizers are different. But I don’t know if there are errors in my code. Please help me to check it.
Posted on
No user picture. yaning Joined: 04/12/2019

How to calculate bivariate h2(or c2, e2) in bivaraite SEM model with two binary traits?

Hi,
I have encountered several problems in doing bivariate SEM models, which is an important part of my graduation project.The main problems are as follows:
Posted on
No user picture. Qiuzhi Xie Joined: 12/27/2018

G*E interaction with a moderator

I am currently doing G*E interaction with SES as a moderator. The plot (attached here) suggests that SES likely moderates with the A factor. To analyze whether the moderation effect is significant on A factor, I drop the moderator on the Path A only and compare the model fit (diffLL). I find that the p value is insignificant. I want to ask whether this is the correct way to see the significance of moderator effect.

Also, is it appropriate to see the unstandardized variance components by SES moderation, rather than the standardized variance components?

Posted on
Picture of user. Julia Joined: 03/29/2012

Nonlinear transformation and estimates

Hi. I have a question about analysis of non-linear transformed variables. I have a continuous variable which is skewed (1.2). Box-Cox transformation suggested 1/y transformation. If I run ACE models (and saturated) on this new variables, how do I interpret and report the results? Is it common to report the estimates for the transformed variable? And is it possible to get the estimates for the original raw variable then?
Thank you in advance!
Julia
Posted on
No user picture. Qiuzhi Xie Joined: 12/27/2018

Question about Confidence Interval for Standardized path coefficients (multivariate model)

Hi everyone, I am new here and learning OpenMX. I have a problem of getting confidence interval for standardized path coefficients. Could anyone please provide me the OpenMx language about this? I also provide my scripts below and hope someone can let me know what the problem is (the script about 'ciACE <- mxCI("VC[1:5, 16:30]")' ) or how to read the output about standardized CI.

# Load Data
nl <- read.csv(file="C:/R/TwinPairACE.csv", header=TRUE, sep=",")

Posted on
Picture of user. chancey Joined: 09/09/2017

Fit tests for ACE models

Good morning,

I am running an ACE model for an ordered categorical outcome using script adapted from [Maes](http://ibg.colorado.edu/cdrom2016/maes/UnivariateAnalysis/one/oneACEo.R). I seem to be able to generate most indicators of model fit, including AIC, BIC, -2LL and p-values associated with likelihood-ratio tests for nested models. However, I cannot seem to generate an associated chi-square value. Is there a simple way to generate chi-square value for these tests? Thanks,

Chance