Parallel standard OpenMx scripts & matching umx versions

This site offers OpenMx and matching umx scripts to fit standard biometrical models to data collected in MZ and DZ twins. The models can be used to estimate the role of genetic (A: additive genetic factors; D: dominance genetic factors) and environmental factors (C: common/shared environmental factors, E: unique environmental factors) to the variance of phenotypes of interest and covariances between phenotypes of interest.

Scripts are organized in pages by number of phenotypes, addition of covariates etc.  Within each page, scripts are organized by type of model in rows (Saturated, ACE estimating variance components, ADE estimating variance components, ACE estimating path coefficients, ADE estimating path coefficients) and by type of data in columns (continuous, binary, ordinal (estimating all thresholds), ordinal (fixing two thresholds and estimating means/variances) using standard code, and using umx (for different data types) in the last column. All scripts source the R code attached here that includes a number of functions that automate various aspects of the models such as labels, starting values, output generated etc.

Note that each of the scripts is represented by a path diagram.  If you click on the diagram, a PDF of the associated script will be displayed in a separate window.  If you click on the filename.R below the diagram, the R script will be downloaded.

Comments, suggestions, corrections welcome!!  Email hmaes@vcu.edu.

New pages/scripts will be added as they are ready to go!  If you have scripts that you'd like to add, send them my way!

Last updated: 02/28/2020