Behavioral Genetics Models

Posted on
No user picture. svrieze Joined: 01/03/2010

Means for growth curve with ACE decomposition

Hello,

How do I specify the intercept and slope means for a growth curve, using an ACE parameterization.

I can follow the RAM notation for growth curves in the documentation, but I'm not sure how to specify the intercept and slope means outside of the RAM framework.

Maybe another way to ask the question is as such: In the ACE factor model examples (common pathway model) the factor means are fixed at zero. How can I change the code to estimate the factor mean?

Thanks,

Scott

Posted on
No user picture. pgseye Joined: 10/13/2009

Covariates in Saturated Model

Hi All,

I'm having similar difficulties adding covariates as in the ACE thread - but I thought this deserved a separate thread as my problem is related to equating means in a saturated model.

I have the model broken down by zygosity - the following code is for MZF. Real values for MZF data are shown after the code, followed by expected values and then estimated (free) parameters after that.

Posted on
No user picture. MariBalt Joined: 08/18/2010

Multigroup Path Analysis with (very) numerous groups

Hello
I am trying to fit multigroup models with 15 groups. When I want to assess the fit of my models, I encounter problems :
(1) on one hand, the fit measures like Chi or RMSEA are never good for any model different from the saturated model

Posted on
No user picture. Sabha Joined: 05/18/2010

Adding covariates to Twin Ace model

Hi ,

I am new to Twin data analysis and would like to get help on how to adjust for covariates ( Age and Gender) in TWIN ACE model to estimate the genetic and environmental variance from the total molecular phenotypic variance . What modifications are needed in the script UnivariateTwinAnalysis_MatrixRaw.R for age and gender adjusted estimation of A, C and E .

Any insight would be of great help.

Sanchita

Posted on
Picture of user. Dorothy Bishop Joined: 02/04/2010

Simple question re Cholesky decomposition

I'm running a multivariate ACE model with raw data and matrices input, based on a version by Hermine Maes.
Changing the start value for the paths from .6 to .5 makes the c path estimates flip from positive to negative.
Am I right in assuming I can just ignore the sign, because the true estimate should be the (positive) square root of the squared path?

Posted on
Picture of user. Dorothy Bishop Joined: 02/04/2010

G x E

I'm trying to understand the analysis used in a paper that has recently created a splash in Science:
Taylor, J., Roehrig, A. D., Hensler, B. S., Connor, C. M., & Schatschneider, C. (2010). Teacher quality moderates the genetic effects on early reading. Science, 328, 512-514.
They have used a moderator variable analysis described by:Purcell, S. (2002). Variance components models for gene-environment interaction in twin analysis. Twin Research, 6, 377-382.
This is in old Mx.

Posted on
No user picture. svrieze Joined: 01/03/2010

Modeling factor means to test SNP associations

I have a single common factor and several manifest phenotypic variables. I'd like to determine the relationship between the factor and several SNPs of interest.

I'm not sure exactly how to model the factor means with respect to SNP status (coded 0,1,2 for homozygotes, heterozygotes, and the other homozygotes). A preliminary thought was to fit a multi-group CFA with one group for each genotype, and allow the means to vary across groups. There must be a better way than this that involves a single-group CFA means model.

Posted on
Picture of user. trinewa Joined: 11/30/2009

Error in MxConstraint after updating to OpenMx 0.3.0

I have just installed the OpenMx version 3.0.3, and corrected the MxConstraint in my multivariate script.
I keep getting error messages when running the MxConstraint in the new script.
Anyone who can help me getting this part of the script right?

# Fit Multivariate Model with Equal Means & Variances across Twin Order and Zygosity
# -----------------------------------------------------------------------
# Constrain expected means and variances to be equal across twin order
equateMeansVarsTwinModel <- mxModel(multiTwinSatFit, name="equateMeansVarsTwin",

Posted on
Picture of user. trinewa Joined: 11/30/2009

Converting ACE script into an ADE script

Question from a novice in OpenMx:

I want to convert an ACE script into an ADE script.

Posted on
No user picture. kpm Joined: 04/16/2010

Binary Moderation

Hello, I'm trying to adapt the moderation script from the twin workshop in Boulder. At the Boulder workshop the binary moderator was divorce (0/1, and identical for twin 1 and twin 2). The moderator I'm using also is scored 0/1 but may differ for twin 1 and twin 2. Everything else is the same, and I changed very little, but instead of estimating moderation effects (aM, cM, eM), those parameters always come out as whatever I put in for the starting values.