Behavioral Genetics Models

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

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.
- Read more about Covariates in Saturated Model
- 2 comments
- Log in or register to post comments

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

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
- Read more about Adding covariates to Twin Ace model
- 37 comments
- Log in or register to post comments

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?
- Read more about Simple question re Cholesky decomposition
- 2 comments
- Log in or register to post comments

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.
- Read more about G x E
- 4 comments
- Log in or register to post comments

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.

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",

Converting ACE script into an ADE script
Question from a novice in OpenMx:
I want to convert an ACE script into an ADE script.
- Read more about Converting ACE script into an ADE script
- 4 comments
- Log in or register to post comments

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.
- Read more about Binary Moderation
- 4 comments
- Log in or register to post comments
Pagination
- Previous page
- Page 33
- Next page