# ---------------------------------------------------------------------
# Program: TwoFactorOrthoRaw-OpenMx100221.R
#  Author: Steven M. Boker
#    Date: Sun Feb 21 14:39:50 EST 2010
#
# This program fits a covariance two factor orthogonal FIML model to the 
#     factorExample1.csv simulated data.
#
#
# ---------------------------------------------------------------------
# Revision History
#    -- Sun Feb 21 14:39:47 EST 2010
#      Created TwoFactorOrthoRaw-OpenMx100221.R.
#
# ---------------------------------------------------------------------

# ----------------------------------
# Read libraries and set options.

options(width=80)

require(psych)
require(OpenMx)

# ----------------------------------
# Read the data and print descriptive statistics.

factorExample1 <- read.csv("factorExample1.csv")

describe(factorExample1)

# ----------------------------------
# Build an OpenMx two factor orthogonal FIML model with fixed variance

indicators <- names(factorExample1)
latents <- c("F1", "F2")
loadingLabels <- c(paste("b_F1", indicators, sep=""), paste("b_F2", indicators, sep=""))
uniqueLabels <- paste("U_", indicators, sep="")
meanLabels <- paste("M_", indicators, sep="")
factorVarLabels <- paste("Var_", latents, sep="")
factorCovLabels <- c("Cov_F1F1", "Cov_F2F2", "Cov_F1F2")


twoFactorOrthoRaw1 <- mxModel("Two Factor Orthogonal FIML Model with Fixed Variance",
    type="RAM",
    manifestVars=indicators,
    latentVars=latents,
    mxPath(from=latents, to=indicators, 
           arrows=1, all=TRUE, 
           free=TRUE, values=.2, 
           labels=loadingLabels),
    mxPath(from=indicators, 
           arrows=2, 
           free=TRUE, values=.8, 
           labels=uniqueLabels),
    mxPath(from=latents,
           arrows=2, 
           free=FALSE, values=1, 
           labels=factorVarLabels),
    mxPath(from="one", to=indicators, 
           arrows=1, free=TRUE, values=.1, 
           labels=meanLabels),
    mxData(observed=factorExample1, type="raw")
    )

twoFactorOrthoRaw1Out <- mxRun(twoFactorOrthoRaw1)

summary(twoFactorOrthoRaw1Out)


# ----------------------------------
# Build an OpenMx two factor orthogonal FIML model with fixed loading


twoFactorOrthoRaw2 <- mxModel("Two Factor Orthogonal FIML Model with Fixed Loading",
    type="RAM",
    manifestVars=indicators,
    latentVars=latents,
    mxPath(from=latents, to=indicators, 
           arrows=1, all=TRUE,
           free=TRUE, values=.2, 
           labels=loadingLabels),
    mxPath(from=indicators, 
           arrows=2, 
           free=TRUE, values=.8, 
           labels=uniqueLabels),
    mxPath(from=latents,
           arrows=2,
           free=TRUE, values=.8, 
           labels=factorVarLabels),
    mxPath(from=latents, to=c("x9", "x7"),
           arrows=1, 
           free=FALSE, values=1),
    mxPath(from="one", to=indicators, 
           arrows=1, free=TRUE, values=.1, 
           labels=meanLabels),
    mxData(observed=factorExample1, type="raw")
    )

twoFactorOrthoRaw2Out <- mxRun(twoFactorOrthoRaw2)

summary(twoFactorOrthoRaw2Out)