Errors in OpenMx 1.1 "Example: Common Factor Model for Ordinal Data"

Posted on
Picture of user. Mike Cheung Joined: 10/08/2009
Forums

Thanks for implementing the analysis of ordinal data in OpenMx 1.1. After running the code, it seems that there are some errors in the example.

Original S matrix:
mxMatrix(
type="Symm",
nrow=4,
ncol=4,
values=c(1,0,0,0,0,0,0,
0,1,0,0,0,0,0,
0,0,1,0,0,0,0,
0,0,0,1,0,0,0,
0,0,0,0,1,0,0,
0,0,0,0,0,1,0,
0,0,0,0,0,0,1),
free=c(T, F, F, F,
F, T, F, F,
F, F, T, F,
F, F, F, F),
labels=c("e1", NA, NA, NA,
NA, "e2", NA, NA,
NA, NA, "e3", NA,
NA, NA, NA, "varF1"),
byrow=TRUE,
name="S"
)

Correct S matrix:
mxMatrix(
type="Symm",
nrow=4,
ncol=4,
values=c(1,0,0,0,
0,1,0,0,
0,0,1,0,
0,0,0,1),
free=c(F, F, F, F,
F, F, F, F,
F, F, F, F,
F, F, F, F),
labels=c("e1", NA, NA, NA,
NA, "e2", NA, NA,
NA, NA, "e3", NA,
NA, NA, NA, "varF1"),
byrow=TRUE,
name="S"
)

Original mxRAMObjective:
mxRAMObjective("A", "S", "F", "M",
dimnames=c("z1", "z2", "z3", "F1"),
thresholds="thresh",
threshnames=c("z1, "z2", "z3"))

Correct mxRAMObjective:
mxRAMObjective("A", "S", "F", "M",
dimnames=c("z1", "z2", "z3", "F1"),
thresholds="thresh",
threshnames=c("z1", "z2", "z3"))

Replied on Fri, 09/02/2011 - 09:00
Picture of user. Ryne Joined: Jul 31, 2009

Thanks, Mike. We spotted those the other day, and we should have new docs up on the site in the next day or so.

Replied on Fri, 09/02/2011 - 20:37
Picture of user. rgore Joined: Jan 27, 2011

In reply to by Ryne

The PDF of the User Guide (http://openmx.psyc.virginia.edu/docs/OpenMx/latest/OpenMxUserGuide.pdf) is now updated. The remaining documentation will be updated when the first bug fixes to OpenMx 1.1 are released.