R Functions and User Interface

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

Error in single.na(to) : argument is missing, with no default

Hi

Anyone have any idea what this error message means?

Error in single.na(to) : argument is missing, with no default
Calls: mxModel -> mxPath -> generatePath -> single.na
Execution halted

I hesitate to post the script, 308 lines at present...

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

the lowly labeled fixed value

I am using the term "labeled fixed value" to refer to a (row, col) entry in a MxMatrix object where free = FALSE and label is not a NA value. They have not generated much discussion so far. The only restriction we have is that a free parameter in the model cannot have the same name as a labeled fixed value. Here are some additional ideas:

  • Labeled fixed values can be used inside an MxAlgebra expression just as we are going to allow free parameters inside an MxAlgebra expression.
Posted on
Picture of user. mspiegel Joined: 07/31/2009

constant + free parameter substitution proposal

There has been some talk recently of simplifying the MxAlgebra interface in two regards. It would be convenient to specify constant values or matrices in an MxAlgebra statement, without having to explicitly declare a separate MxMatrix that contains constant values. It would also be convenient to specify free parameters in an MxAlgebra statment, without having to explicitly declare a separate 1 x 1 MxMatrix that contains the free parameter. We have been referring to these proposed changes as constant substitution and free parameter substitution, respectively.

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

new 'compute' flag to mxEval()

If 'compute' is FALSE (the default), then MxAlgebra expressions returns their current value as they have been computed by the optimization call (using mxRun). If the 'compute' argument is TRUE, then MxAlgebra expressions will be calculated in R. Any references to an objective function that has not yet been calculated will return a 1 x 1 matrix with a value of NA. Users will like this feature because it allows them to compare the result of algebraic expressions as computed by the back-end against those values computed by R.

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

means vector

The means vector is now always converted into an n x 1 matrix. The argument to mxData continues to accept an R numeric vector (in fact, I don't think the mxData function currently will accept a matrix for the means argument). But the means are now always stored as an n x 1 matrix. We had some sloppiness in the R interface where sometimes it was n x 1 and sometimes it was 1 x n. I updated the test suite in light of these changes, but I did not update the /docs.

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

Ordinal Data Interface

At the last developer's meeting, it was decided that the ideas presented in Ordinal Data Proposal #2 would be implemented. The mxFIMLObjective(), mxMLObjective(), and mxRAMObjective() functions have been updated with a new argument "thresholds" (in the subversion trunk). This new argument is a character string that will be the name of the MxMatrix that will serve as the thresholds matrix. The thresholds matrix is an r x c matrix where r equals the length of the data with the most levels minus one and c must equal the dimensions of the covariance matrix.

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

mxEval

Topic for discussing mxEval

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

loading packages conditionally

The OpenMx repository has a work-in-progress feature that generates graphs from RAM style models. This package requires the Rgraphviz library. I don't want to add Rgraphviz as a required library. Not all users will need to make graphs, and it turns out that the Rgraphviz library can be a pain to install on Windows (you need a specific version of the Graphviz library that is hidden on their website). So this feature conditionally works if the Rgraphviz package is installed.

Posted on
Picture of user. pdeboeck Joined: 08/04/2009

Bootstrapping Confidence Intervals

I couldn't find info on bootstrapping confidence intervals on the site, so I wanted to put in a plug for such an option. Kris Preacher and I were talking about this yesterday.