OpenMx Developer Forums

Posted on
No user picture. CharlesD Joined: 04/30/2013

Change the 'running xxxmodel' output from mxRun to a message

instead of cat or print or however it may be output to the console at the moment, the message function interacts nicer with other R / knitr functionality, I think.
Posted on
Picture of user. tbates Joined: 07/31/2009

allow mxOption to get as well as set

currently

mxOption(m1, key="Number of Threads")

and

mxOption(m1)

Both fail (no value, no key errors, respectively)

It would be nice if asking for a key returned its value, and asking about a model returned the currently active options

Posted on
No user picture. Anonymous (not verified) Joined: 09/16/2024

Using analyticaly derived gradient and hessian in OpenMX

I know NPSOL will accept first derivatives (gradient) and second derivatives (Hessian) of the fitting function with respect to the parameters.

If I have determined these analytically for my model, can i make OpenMX parse them to NPSOL?

The only slot I can find for the gradient and the hessian are in a fitted model object not in an model hat still has to be fitted.

I did find the mxOption slot where i can tell NPSOL the gradient and hessian will be provided and need not be estimated

Michel

Posted on
Picture of user. tbates Joined: 07/31/2009

diag()

I'd really like diag implemented in Openmx - just so that the same code works inside and outside of OpenMx and to ease the learning transition. But, If we aren’t going to support the standard r implementation of diag(), can we implement


diag <- function(x=1, nrow, ncol) {
message("We’re a bit worried that diag can be ambiguous, so we’ve broken it into multiple new functions...")
# then either say
stop("You need to choose among list up our replacement functions to handle everything that diag() does...")

Posted on
Picture of user. neale Joined: 07/31/2009

Backward compatibility issue with extracting or inserting certain results of optimization into mxModel objects

The attached file runs fine under 1.3 but fails unpleasantly under the current svn trunk (2484).

> LCA2 <- LCAfun(data[4:8], 2, 5, 3)
Error in as.character.default() :
no method for coercing this S4 class to a vector
In addition: Warning messages:
1: In mxFIMLObjective(covariance = "R", means = "M", dimnames = nameList[1], :
Objective functions have been deprecated. Please use mxExpectationNormal() and mxFitFunctionML() instead.
2: In mxFIMLObjective(covariance = "R", means = "M", dimnames = tempVar, :

Posted on
Picture of user. mhunter Joined: 07/31/2009

Defunct omxRAMExpectation.h file?

I noticed that the typedef struct for the RAM expectation was moved its .h file to the .c file. The comment for the svn checkin says "Improve encapsulation and remove deadcode". Could anyone explain to me how this improves things? Speed ups? Bug removals? New features? Single point of truth?

Posted on
Picture of user. Ryne Joined: 07/31/2009

Many many submodels (how S4 builds objects and a request to run lists of models)

So I'm following up on some of my factor score work and happened across and old question I had regarding holder or parent models for embarassingly many submodels.

I'm trying to run the same model on 500 (or some other large n) datasets with the independence flag set to TRUE. I'm doing this by creating one model per dataset, building a list of mxModels, then putting this list into a single mxModel for optimization. However, I'm spending most of my time the fourth line of this code:


singScore <- transformFactorScores(spRes, 1, "mu", "sigma", "epsilon")

Posted on
Picture of user. mhunter Joined: 07/31/2009

Updated Copyright

I updated the copyright on the OpenMx website, and the source code. The source code update was done with a few lines of sed at the Unix terminal.

find . -name "*.h" -exec sed -i "s/Copyright 2007-2012 The OpenMx/Copyright 2007-2013 The OpenMx/g" '{}' \;

This finds files in the current directory (.) with the name ending in ".h" and executes sed inplace (-i) editing the copyright (replace "Copyright 2007-2012 The OpenMx" with "Copyright 2007-2013 The OpenMx"). Analogous commands where done in models/passing, demo/, man/, and src/.

Posted on
Picture of user. mhunter Joined: 07/31/2009

Developers Meeting Notes 2/1/13

At the developers meeting on February 1 we discussed the following.

  • The meeting had more than typical technical difficulties.
  • We were introduced to Mahsa Zahery, a computer science graduate student at VCU. Her two main tasks are to allow GPUs to work on numerical integration routines, and to help make OpenMx more finely grained parallel.
  • We discussed the difficulty of Mahsa Zahery's second task.
Posted on
Picture of user. mhunter Joined: 07/31/2009

Developers Meeting Notes 1/25/13

The OpenMx development team discussed the following on January 25:

  • Andreas Brandmaier, Tim Brick, and Timo von Oertzen presented a demo version of an mxPlot command that calls Onyx to draw a diagram corresponding to an MxModel object. I believe the demo was largely made by Andreas. Several aspects of the user interface were discussed at length.