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

2 replies [Last post]
Mike Cheung's picture
Offline
Joined: 10/08/2009

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"))

Ryne's picture
Offline
Joined: 07/31/2009
Thanks, Mike. We spotted

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

rgore's picture
Offline
Joined: 01/27/2011
PDF User Guide updated

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.