Behavioral Genetics Models

Is anyone aware of a script to do a Cholesky decomposition between a binary and a continuous variable?

I attempted to adapt a script I found here (http://www.genepi.qimr.edu.au/staff/sarahMe/workshop14/Day3.zip) and listed it below, but it's not running.

require(OpenMx)
require(psych)
source( "http://www.vipbg.vcu.edu/~vipbg/Tc24/GenEpiHelperFunctions.R")

nl <- read.delim("Cholesky.dat", header=T, sep="\t", na.strings=".")

nl$family1 <- nl$fam
nl$family2 <- nl$famb
nl$happy1 <- nl$max
nl$happy2 <- nl$maxb

# Select Variables for Analysis

Hi there,
I want to test for equal means and variances across twins and zygosities in a trivariate model, separately for the measures.

In the trivariate script I am using, http://ibg.colorado.edu/cdrom2014/bartels/Multivariate/Trivariate.R , the means and variances are tested together and I am not sure how to separate them.

Can I use my results from univariate saturated sub models or is there a way to test for it in a trivariate model?

Many thanks!

Hi,
I am trying to get the CI for correlations in trivariate model. I did get them in my univariate models and tried to readjust the script to my multivariate models:

# Algebra for expected Means, Covariances and Correlation Matrices in MZ & DZ twins
Saturated_Model <- mxModel("Saturated",
mxModel("MZM", mxMatrix( type="Full", nrow=1, ncol=ntv, free=TRUE,
values=svMe, labels=c("MZM1_1","MZM1_2","MZM2_1","MZM2_2","MZM3_1","MZM3_2"),
name="expMeanMZM" ),

mxMatrix( type="Symm", nrow=ntv, ncol=ntv, free=TRUE,
values=c(68,6,1,41,4,1,
6,8,3,5,2.7,2.5,
1,3,12,1,2,3.5,

Hi,
I am conducting a trivariate twin analysis and after running ACE model my estimates of C for males are significant for phenotypes 1 and 3 but not significant for phenotype 2. AE model results in significant worsening of the goodness of fit.
I have run ACE model with c21, c22 and c32 fixed to 0 using the attached script.
Is that correct way of doing it or is there another way?
Any comments will be appreciated
Many thanks!

Hi everyone,

Greetings, fairly new student of OpenMx here. I've successfully run a univariate twin model with confidence intervals thanks to mxCI, but in the quadvariate script below (obtained from http://ibg.colorado.edu/cdrom2014/bartels/Multivariate/Trivariate.R), I don't know where to put the mxCI command to request CIs for A/V, C/V, and E/V. The output is in the form of a matrix too, so I am unsure how to proceed. Any assistance would be greatly appreciated.

# Matrices declared to store a, c, and e Path Coefficients
pathA <- mxMatrix( type="Lower", nrow=nv, ncol=nv, free=TRUE,

Hi,
I am new to OpenMx and trying to use it for twin study (monozygotic vs dizygotic twins). I have some questions and hope someone may help me with the following questions.

I tried the following two scripts:
http://openmx.psyc.virginia.edu/svn/tags/stable-1.2/demo/UnivariateTwinAnalysis.R
http://openmx.psyc.virginia.edu/svn/tags/stable-1.2/demo/UnivariateTwinAnalysis_MatrixRaw.R

I knew their input data were different, but is there any difference for data processing steps between these two?

I'm trying to pull out estimates in a matrix using genEpi helper functions in a trivariate model with a sibling but I'm running into some trouble with the genEpi_FormatOutputMatrices function.

Error:
Error in print(genEpi_FormatOutputMatrix(matrix = genEpi_EvalQuote(expstring = matricesList[[k]], :
error in evaluating the argument 'x' in selecting a method for function 'print': Error: The following error occurred while evaluating the expression 'iSD %*% a' in model 'CholACE' : requires numeric/complex matrix/vector arguments

Full Script is below

require(OpenMx)

Hi all,

Hi everyone!
Sorry to ask another question again. I'd like to analyse the trait correlation, for example BMI and WC, which function should I try to use or is there any scripts which I can study. Thanks a lot.