OPTIMIZATION FOR STRUCTURAL EQUATION MODELING: APPLICATIONS TO SUBSTANCE USE DISORDERS

Substance abuse is a serious issue in both modern and traditional societies. Besides health complications such as depression, cancer and HIV, social complications such as loss of concentration, loss of job, and legal problems are among the numerous hazards substance use disorder imposes on societies. Understanding the causes of substance abuse and preventing its negative effects continues to be the focus of much research. Substance use behaviors, symptoms and signs are usually measured in form of ordinal data, which are often modeled under threshold models in Structural Equation Modeling (SEM). In this dissertation, we have developed a general nonlinear optimizer for the software package OpenMx, which is a SEM package in widespread use in the fields of psychology and genetics. The optimizer solves nonlinearly constrained optimization problems using a Sequential Quadratic Programming (SQP) algorithm. We have tested the performance of our optimizer on ordinal data and compared the results with two other optimizers (implementing SQP algorithm) available in the OpenMx package. While all three optimizers reach the same minimum, our new optimizer is faster than the other two. We then applied OpenMx with our optimization engine to a very large population-based drug abuse dataset, collected in Sweden from over one million pairs, to investigate the effects of genetic and environmental factors on liability to drug use. Finally, we investigated the reasons behind better performance of our optimizer by profiling all three optimizers as well as analyzing their memory consumption. We found that objective function evaluation is the most expensive task for all three optimizers, and that our optimizer needs fewer number of calls to this function to find the minimum. In terms of memory consumption, the optimizers use the same amount of memory

OpenMx 2.0:Extended Structural Equation and Statistical Modeling

The new software package OpenMx 2.0 for structural equation and other statistical modeling is introduced and its features are described. OpenMx is evolving in a modular direction and now allows a mix-and-match computational approach that separates model expectations from fit functions and optimizers. Major backend architectural improvements include a move to swappable open-source optimizers such as the newly-written CSOLNP. Entire new methodologies such as Item Factor analysis (IRT) and State-space modeling have been implemented. New model expectation functions including support for the expression of models in LISREL syntax and a simplified multigroup expectation function are available. Ease-of-use improvements include helper functions to standardize model parameters and compute their Jacobian-based standard errors, access to model components through standard R $ mechanisms, and improved tab completion from within the R Graphical User Interface

A Bivariate Genetic Analysis of Drug Abuse Ascertained Through Medical and Criminal Registries in Swedish Twins, Siblings and Half-Siblings

Using Swedish nationwide registry data, the authors investigated the correlation of genetic and environmental risk factors in the etiology of drug abuse as ascertained from medical and criminal registries by modeling twin and sibling data. Medical drug abuse was defined using public inpatient and outpatient records, while criminal drug abuse was ascertained through legal records. Twin, full and half sibling pairs were obtained from the national twin and genealogical registers. Information about sibling pair residence within the same household was obtained from Statistics Sweden. Standard bivariate genetic structural equation modeling was applied to the population-based data on drug abuse ascertained through medical and crime registries, using OpenMx. Analyses of all possible pairs of twins (MZ: N = 4482; DZ: N = 9838 pairs), full- (N = 1,278,086) and half-siblings (paternal: N = 7767; maternal N = 70,553) who grew up together suggested that factors explaining familial resemblance for drug abuse as defined through medical or criminal registries were mostly the same. Results showed substantial heritability and moderate contributions of shared environmental factors to drug abuse; both were higher in males versus females, and higher for drug abuse ascertained through criminal than medical records. Because of the low prevalence of both assessments of drug abuse, having access to population data was crucial to obtain stable estimates. Using objective registry data, the authors found that drug abuse—whether ascertained through medical versus criminal records—was highly heritable. Furthermore, shared environmental factors contributed significantly to the liability of drug abuse. Genetic and shared environmental risk factors for these two forms of drug abuse were highly correlated