algebra names of MxFIMLObjective

6 replies [Last post]
brandmaier's picture
Offline
Joined: 02/04/2010

Just a minor issue:
Looking into the output of "objective" of an MxFIMLObjective, I find that the covariance matrix is names "covariance" but the means vector only has a generic title like "untitled4391". For consistency, I suggest to name the means vector "means".

compare this output:

> model@objective
MxFIMLObjective 'objective'
@covariance : 'covariance'
@means : 'untitled4391'
@vector : FALSE
@dims : NA
@thresholds : NA

best,
Andreas

mspiegel's picture
Offline
Joined: 07/31/2009
Ah, this is not a bug. There

Ah, this is not a bug. There is already a free parameter in the script that has the name 'means'. This script is one of our test cases, it used to throw an error complaining that a free parameter and a matrix were assigned the same name. There is now a check that looks for anything with the name 'means'. If that is found, then an anonymous name is used. I could change the anonymous name from:

 # formerly known as omxUntitledName()
meansName <- imxUntitledName()

to

meansName <- paste('means_', imxUntitledName(), sep = '')

Is that OK?

brandmaier's picture
Offline
Joined: 02/04/2010
Thanks again for the fast

Thanks again for the fast response!

My feeling is that "means_untitledXXX" is a more intuitive solution.

mspiegel's picture
Offline
Joined: 07/31/2009
Can I confirm that you were

Can I confirm that you were using the RAM objective function with raw data and that was transformed into the FIML objective function? That is the first workflow that comes to mind that could auto-generate a means vector for you.

mspiegel's picture
Offline
Joined: 07/31/2009
Wait, I just tested RAM + raw

Wait, I just tested RAM + raw data and that generates a means vector that is named "means". Could you let us know what is the workflow of your script?

brandmaier's picture
Offline
Joined: 02/04/2010
I am using $mxVersion

I am using $mxVersion "1.0.4-1540".

I attached a script which builds a linear LGCM via RAM notation and uses some fake raw data. For me, I get still an "untitledXXX" as means vector name.

AttachmentSize
meansVectorNames.R 1.67 KB
mspiegel's picture
Offline
Joined: 07/31/2009
Aha. I have a better idea of

Aha. I have a better idea of where the problem lies. I was testing a matrix type model that uses the RAM objective function. I should have been testing a RAM type model. I'll take a look at this over the next couple of days. The fix will be in the next OpenMx pre-release as well as in the next major release series (1.1.0 or whatever numbering system is adopted).