Behavioral Genetics Models

Dear all,

I am trying to add sibling and parents in my heritablity estimation study. But I know little about that. I am wondering do you any example Openmx script of ACE or ADE model which involved not only MZ and DZ but also sibling and parents? If so, could you please send it to me?. Or if you know there is some links or materials which can help me learn by myself, could you also send them to me?
Thank you very much!

I am wondering if someone can answer a nagging question I have about the best way to deal with outliers in heritability modeling.

It has been suggested to me that using z scores, anyone outside the cut off of +/- 2.5 should have their score altered to be the next extreme score (eg if we have scores of 2.5, 2.9, 3.2 etc these would be 2.5, 2.6, 2.7). I am concerned that this will really alter the heritability results. Any thoughts?

Is is better to just remove any outliers from the analyses?

Hello

I am wondering if it is possible to apply the age moderation script from Boulder (attached) to the Common Pathways model? I have tried to do it but end up with multiple problems (one of which is how to add an age moderation to the f component - but that is almost the least of my worries) - my attempt is also attached. Any help would be greatly appreciated.

Thank you in advance

Karen

Hello

My data of cognitive measures is indicating a number of violations of the assumptions of the twin model. I have transformed the data using z scores and still find that MZ and DZ variances are different. Any ideas on how I should proceed, including how to adjust the script in the ACE/ADE model to deal with this?

Many thanks

Karen

This is the first time I use open Mx.
without such covariates, my script runs.

Now, I add age and sex as definition factors. results show :
'Unknown thresholds name 'expThreMZ1' detected in the objective function of model 'MZ'.

I don`t know what is wrong. So if you could give me some help? Thanks.

=================================================
MeanL <-mxMatrix( type="Zero", nrow=1, ncol=ntv, name="M" )

Inc <-mxMatrix( type="Lower", nrow=nth, ncol=nth, free=FALSE, values=1, name="L" )
Tmz <-mxMatrix(type="Full", ........, name="ThMZ" )

Hi,
I am trying to fit a biv model (2 continuous phenotypes) with sex as definition variable to test for mean differences, however I get the same error over and over again:

Error: A definition variable has been declared in model 'Chol' that does not contain a data set

I cannot find where it goes wrong, I suspect it's somewhere in the matrices declared to store linear coefficients for covariate:

grandMean <- mxMatrix(type="Full", nrow=1, ncol=nphen, free = TRUE, values=c(2700, 19), label=c("mean1","mean2"), name="Mean")

Hi,
I am quietly new to Openmx and thankful for any help that I can get.
I have repeated measurement of continuous variable over two time periods to which I would to do similar analysis as in the De Geus et al (2007) : Bivariate Genetic Modeling of Cardiovascular Stress Reactivity: Does Stress Uncover Genetic Variance?.

Hi,

I have fitted an ADE (threshold) model for two separate groups: an adolescent and an adult group.
The estimates of A, D and E differ between the two groups so I would like to test if this difference is significant.
I thought this might be possible by putting both groups together into one model, and then equating
the A, D & E estimates one by one to see if this significantly deteriorated the fit.

Hi,

Hi,

I am trying to carry out heritability analysis with gene expression data.