OpenMx General Help

Posted on
Picture of user. carey Joined: 10/19/2009

* and %x% operators in mxAlgebra

anyone explain this?

> model2 <- mxModel(model1, mxAlgebra(.5*CTf, name="CTDZf"))
> model2 <- mxRun(model2)
Running SameMZDZf
Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
BLAS/LAPACK routine 'DPOTRF' gave error code -4

> model2 <- mxModel(model1, mxAlgebra(.5 %x% CTf, name="CTDZf"))
> model2 <- mxRun(model2)
Running SameMZDZf
Warning message:
In model 'SameMZDZf' NPSOL returned a non-zero status code 1. blah blah ...

yet
> mat <- matrix(rnorm(16), 4)
> .5 * mat == .5 %x% mat
[,1] [,2] [,3] [,4]
[1,] TRUE TRUE TRUE TRUE

Posted on
No user picture. xlzqz Joined: 02/08/2011

Where can I download the source codes.

Does anybody know how can we download the source codes of OpenMx. Looks like no links to the source codes from the OpenMx Website.
Posted on
Picture of user. Lawrence L Lo Joined: 02/24/2010

Appending mxAlgebra expressions

Hello,

I am trying to generate mxAlgebra expressions from a size variable. I have a variable amount of mx submodels and am trying to add the objectives from these together for a full model objective. I have tried things like
mxAlgebra(sum(get(paste("submodel",1:x,".objective",sep=""))),
name="TotObj")
and other various methods. I also tried the method described on
http://openmx.psyc.virginia.edu/thread/537
but don't think this will help in my particular situation.
Any suggestions?

Best,
-LLL

Posted on
Picture of user. trinewa Joined: 11/30/2009

How do I add a demo script to an OpenMx directory?

Hi
I woulld like to share a demo script to the following directory:
http://openmx.psyc.virginia.edu/repoview/1/trunk/demo

The Wiki says where to put the file, but not technically how to add the file to the directory.
Any help?

Posted on
Picture of user. carey Joined: 10/19/2009

@submodels slot

run:

MZf <- mxModel(name="MZf", mxData(ObsCovMZf, type="cov", numObs=nMZf),
mxMLObjective(covariance='myModel.PreCovMZf'))
DZf <- mxModel(name="DZf", mxData(ObsCovDZf, type="cov", numObs=nDZf),
mxMLObjective(covariance='myModel.PreCovDZf'))
myModel1 <- mxModel(name="myModel", MZf, DZf)
myModel1

then run and compare:

MZf <- mxModel(name="MZf", mxData(ObsCovMZf, type="cov", numObs=nMZf),
mxMLObjective(covariance='myModel.PreCovMZf'))
DZf <- mxModel(name="DZf", mxData(ObsCovDZf, type="cov", numObs=nDZf),

Posted on
Picture of user. carey Joined: 10/19/2009

remove = TRUE issue

try the following code

VA <- mxMatrix(name="VA", type="Symm", nrow=5, ncol=5)
DZ <- mxAlgebra(.5 %x% VA, name="covDZA")
test <- mxModel(model="test", VA, DZ )
test
names(test)
rm1 <- mxModel(test, test@matrices$VA, remove=TRUE)
rm2 <- mxModel(test, test@algebras$covDZA, remove=TRUE)
rm3 <- mxModel(test, test@algebras$DZ, remove=TRUE)

why doesn't remove=TRUE work with rm1 and rm2 while it fails to throw an error or warning with rm3?

Posted on
Picture of user. carey Joined: 10/19/2009

MxModel@intervals

some old (over 1 year ago) MxModel objects that i saved are giving me the following error in R:
Error: no slot of name "intervals" for this object of class "MxModel"

if i create a new MxModel, it shows a slot for intervals of type list. but i cannot find documentation for this slot in http://openmx.psyc.virginia.edu/docs/OpenMx/latest/_static/Rdoc/MxModel-class.html (or in the doc for mxModel either).

any help?
best,
greg

Posted on
No user picture. deb193 Joined: 12/01/2010

quiet failure without error message

I am not getting an error message, but the model seems to be not running. I run the ful model and get a GREEN status, then I set the paths moderatign the loadings to 0 and run the nested model ... and nothing. It even reports 0 observed statistics.

How can I troubleshoot when there is no eror message?

------------------------------------------------------------------------------------

> summary(FactorConfigFIout)
data:
$`Moderated FI Model.data`

Posted on
No user picture. Hermine Joined: 07/31/2009

precision of integration

Is it possible to change the precision of the numerical integration for binary FIML? I don't get the same likelihood in classic versus OpenMx.

Thanks, Hermine.

Posted on
No user picture. pgseye Joined: 10/13/2009

Set probability precision

Hi,

Is there a way to set the precision of the reported p values to more than 2 decimal places. It might be something simple, but I can't find it in the documentation.

Thanks,

Paul