Running OpenMx with multicore

openmc.R

require(OpenMx)
library(parallel)

data(demoOneFactor)
manifests <- names(demoOneFactor)
latents <- c("G")

factorModel <- mxModel(name="One Factor",
    type="RAM",
    manifestVars = manifests,
    latentVars = latents,
    mxPath(from=latents, to=manifests),
    mxPath(from=manifests, arrows=2),
    mxPath(from=latents, arrows=2, free=FALSE, values=1.0),
    mxData(observed=cov(demoOneFactor), type="cov", numObs=500)
)

# limit to 1 thread for multicore
mxOption(NULL, 'Number of Threads', 1)

# vary number of observations for data
factorFits <- mclapply(seq(500, 1000, by=100),
                      function(obs, model) mxRun(mxModel(model, mxData(observed=cov(demoOneFactor), type="cov", numObs=obs))),
                      model=factorModel, mc.cores=6)