Behavioral Genetics Models

Hello everyone!

During the process of getting the result of my twin submodel, I failed to get the CI. For example, in the AE model, I get the a2, e2, but how to show the CI? Of course, I get the CI in the full model ACE. The scripts of submodel are below, should I add something?

# Run AE model
AeModel <- mxModel( AdceFit, name="AE" )
AeModel <- omxSetParameters( AeModel, labels="d11", free=FALSE, values=0 )
AeModel <- omxSetParameters( AeModel, labels="c11", free=FALSE, values=0 )
AeFit <- mxRun(AeModel)
mxCompare(AdceFit, AeFit)
round(AeFit@output$estimate,4)

Hello

Hi all,

Hello,

I was wondering if p-values for estimated genetic correlations can be obtained directly from the model output, in addition to confidence intervals?

I know that a genetic correlation is significant when zero is not included in the CI, but I need a p-value for multiple comparison correction as I will be running many models.

Any help would be much appreciated,

Marc

Hello

I downloaded from the OpenMX website the correlated factors sex limitation model, and amended it for our data (removing the opposite sex components as our twins are all same sex).

The script runs and converges but there is no output (seeming that it is not fitting). There are no warning messages and no error messages. The output we get is:

observed statistics: 4500
estimated parameters: 0
degrees of freedom: 4500
-2 log likelihood: NA
saturated -2 log likelihood: NA
number of observations: 750
chi-square: NA
p: NA
AIC (Mx): NA
BIC (Mx): NA

Hello,
I wanted to work the following as an example of how to apply SEM in genes scenario
http://www.biomedcentral.com/1753-6561/1/S1/S76

or

https://www.google.mu/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&cad=rja&ved=0CC8QFjAA&url=http%3A%2F%2Fwww8.sas.com%2Fscholars%2FProceedings%2F2006%2FSTAT%2FST06_06.PDF&ei=1xl9UvehHYnwrQec34Ew&usg=AFQjCNFEKotV1T6cZ_sJVS8N1oHYrXwDAg&bvm=bv.56146854,d.bmk

Hi all,

I would like to covariate for age & gender in a multivariate analysis. I found scripts to covariate but in ordinal data and in univarite analysis.

It would be really helpful for me it any of you could show me an example of a script on how to covariate (by age and gender) in a MV analysis, because I now from my previous univariate analysis that gender could be causing some genetic sex differences in 4 of the 6 traits that I’m using.

Thank you very much!!!

Hi all,
I’m trying to run a saturated model with restrictions in a multivarite analysis with 6 continuous variables.

In this saturated model with restrictions (see below) I had the following problem (that I did not have in the saturated model without restrictions):
Error: The job for model 'Con' exited abnormally with the error message: Expected covariance matrix is not positive-definite in data row 431 at major iteration 0 (minor iteration 1).

purcell_Transf&DF.R

Dear OpenMx users,

For those interested in DF analysis, we have written a script which
in part 1 prepares / transforms data as in Purcell's means.gawk and prepare.gawk scripts
and in part 2 runs DF-extremes analysis.

Ref: Maximum likelihood DF-extremes analysis based on Purcell S. & Sham PC, Behav Genet. 2003 May;33(3):271-8.
"A model-fitting implementation of the DeFries-Fulker model for selected twin data".

Please let us know if you have any problems or questions.