# --------------------------------------------------------------------- # Program: BivariateStd-OpenMx100214.R # Author: Steven M. Boker # Date: Sun Feb 14 12:12:17 EST 2010 # # This program fits a bivariate model to the multiData simulated data. # # # --------------------------------------------------------------------- # Revision History # -- Sun Feb 14 12:12:15 EST 2010 # Created BivariateStd-OpenMx100214.R. # # --------------------------------------------------------------------- # ---------------------------------- # Read libraries and set options. options(width=80) require(psych) require(OpenMx) # ---------------------------------- # Read the data and print descriptive statistics. multiData1 <- read.csv("multiData.csv") describe(multiData1) # ---------------------------------- # Build an OpenMx bivariate regression model using y and x1 manifests <- c("x1", "x2", "y") multiData1Cov <- cov(multiData1[,c(1,2,5)]) biRegModel <- mxModel("Bivariate Regression of y on x1 and x2", type="RAM", manifestVars=manifests, mxPath(from=c("x1","x2"), to="y", arrows=1, free=TRUE, values=.2, labels=c("b1", "b2")), mxPath(from=manifests, arrows=2, free=TRUE, values=.8, labels=c("VarX1", "VarX2", "VarE")), mxPath(from="x1", to="x2", arrows=2, free=TRUE, values=.2, labels=c("CovX1X2")), mxData(observed=multiData1Cov, type="cov", numObs=500) ) biRegModelOut <- mxRun(biRegModel) summary(biRegModelOut)