Suggestion Box

Hi

I just discovered now OpenMx, whihc seems to be really nice! So thanks in advance:-)

Is there a reason the package is not on CRAN and not listed in the CRAn task-views?

I think this would be a significant improvement, especially if it can be on the task-views... I had done quite a lot of searches on R packages for sem and did not find OpenMx on usual R search tools (mailing list, task-view), just found it randomly. So could be of advantage to be listed on task-views!!

Thanks

is there internal documentation on the effects of mxOptions and resultant quantities in the @output slots of an MxModel object?

specific problem involves the elements of ...@output$hessianCholesky and ...@output$hessian. for the simple problem i outlined in "summary(mxRun(myMxModel)) gives strange results" i get identical results using all mxOptions defaults and using
this <- mxOption(this, "Hessian", "No"). that is, both ...@output$hessianCholesky and ...@output$hessian give identical results.

(1) run the attached. (NOTE: data sets at the beginning should be into R matrix objects ObsCovMZ and ObsCovDZ). then examine the result of "summary(test2)." i get (after the parameters):

Observed statistics: 0
Estimated parameters: 45
Degrees of freedom: -45
-2 log likelihood: 29434.03
Saturated -2 log likelihood:
Chi-Square:
p:
AIC (Mx):
BIC (Mx):
adjusted BIC:
RMSEA:

(2) suggest placing the chi square and contribution to chi square for each group in the MxModel object, so users can access these quantities directly.

there are inconsistencies in the documentation and execution of a model involving dimnames for matrices

i read in an observed cov matrix using read.table(..., header=TRUE, ...) and got this error in mxRun:
>Running thisOMxModel
>Error: The dataset associated with the ML objective in model 'MZ' does not contain identical row and column non-NULL dimnames.

fine. fixed that up, and the next mxRun gave this error:
>Running thisOMxModel
>Error: The expected covariance matrix associated with the ML objective in model 'MZ' does not contain dimnames.

again, pardon intrusions from an R neophyte.

# this works
thisOMxModel <- mxModel("thisOMxModel", mxMatrix(name='VA', type='Symm', free=TRUE, nrow=8, ncol=8) )

# but this does not
thatOMxModel <- new("MxModel")
thatOMxModel <- mxModel(thatOMxModel, mxMatrix(name='VA', type='Symm', free=TRUE, nrow=8, ncol=8) )

greg

(1) i know that forums are for discussion, but i cannot find a place on the website for potential bug reports. would a forum of this name be useful?

(2) here is a place to start. submit the following:
library(OpenMx)
thisOMxModel <- mxModel(name="testModel")
thisOMxModel <- mxModel(thisOMxModel, mxMatrix(name="Afac", type="Full",nrow=8,ncol=5))
Efac <- mxMatrix(name="Efac", type="Full",nrow=8,ncol=5)
thisOMxModel <- mxModel(thisOMxModel, Efac)
testit <- mxModel(thisOMxModel, mxMatrix(name=Afac), remove=TRUE)
testit <- mxModel(thisOMxModel, Efac, remove=TRUE)

(1) am inexperienced in R, so please ignore if this is stupid
(2) the scope of some OpenMx objects allows an object to be replaced by an identically-named object in a user's R session, function, whatever. run the following and note the error "Type must be one of ... "

library(OpenMx)
thisOMxModel <- mxModel(name="testModel")
thisOMxModel <- mxModel(thisOMxModel, mxMatrix(name="Afac", type="Full",nrow=8,ncol=5))
matrixTypes <- c("Full", "Symm: Symmetric", "Diag: Diagonal",
"Stand: Correlation", "Lower: Diagonal and below", "Sdiag: Sub-diagonal",

Can we add an "Administration" topic to the OpenMx Developer Forums area? I wanted to ask about the revision number for Friday's release, but there is no good place to ask that question.

This forum is a place for you to give us suggestions about any part of the OpenMx and OpenSEM projects. If you have suggestions about making the web site better, or changing the organization of the forum topics, or services we could provide that we aren't providing, this is the place. We also welcome general comments about the projects. We want to do our best to provide the most helpful community area we can, and your ideas will make it better!