Is the $IndependenceLikelihood also the -2 * likelihood from the independence model?

In my case it is very different from the $minus2Likelihood of fitting the independence model itself.

-2loglik of independence model : 3096.288

IndependenceLik of target model: 126040060

For the independence model I used a diagonal model.

See http://openmx.psyc.virginia.edu/thread/78

Per Hao's post http://openmx.psyc.virginia.edu/thread/1061 there is definitely odd behavior with omxMnor() with singular and non-positive definite covariance matrices. I seem to remember that Genz's integration routines were extended to handle positive semi-definite matrices (one or more eigenvalues of zero) by reducing the dimensionality of the problem. However, the answer is clearly wrong in this case:

> omxMnor(array(1,dim=c(2,2)),cbind(0,0),cbind(-Inf,-Inf),cbind(0,0))
[,1]
[1,] 0.375
(should be .5)

1.) Bug in type='RAM'. There's no slot (and no argument) in mxRAMObjective for threshold names. The attached model yields the following error message during the mxRun call for the first model (oneFactorModel1):

Running Common Factor Model Path Specification
Error in updateThresholdDimnames(.Object, job, flatJob, model@name) :
no slot of name "threshnames" for this object of class "MxRAMObjective"

Make sure that the help page for mxAlgebra() (?mxAlgebra) contains all the functions that are listed in data/omxSymbolTable.tab. This should be done for both the /trunk and /tags/stable-1.1. I don't believe there are any functions in the trunk that are missing from the 1.1 branch, but you should check.

I had thought that there would be a default order of operations of mxAlgebra() operators. Indeed, consultation of the wiki http://openmx.psyc.virginia.edu/wiki/operator-precedence seems to confirm that binary + should be done after %*%

I am trying to fill either a standardized or a subdiagonal (m x m) matrix with values read from an R matrix. The matrix has 26 rows & 26 columns, so there should be (m*(m-1))/2 = 325 elements in the subdiagonal part. I note that the documentation says:

When ‘type’ is either ‘Lower’, ‘Sdiag’, ‘Symm’, or ‘Stand’, and the arguments to ‘free’, ‘values’, ‘labels’, ‘lbound’, or ‘ubound’ are vectors with enough elements to populate exactly one half of the matrix, then mxMatrix() populates the lower triangle of the matrix (and transposes the lower triangle if the matrix is symmetric).

When I run the attached script I get:

Running model
Error in strsplit(flatName, imxSeparatorChar, fixed = TRUE) :
non-character argument

This is because the algebra for threshold matrix is mis-spelled Thresholds instead of thresholds. It would be much nicer to get an error message that rudely informs me that I can't spell - or at least that I've referred to an object that doesn't exist :). Data are here: http://openmx.psyc.virginia.edu/repoview/1/trunk/models/passing/data/mddndzf.dat