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

Posted on Fri, 09/02/2011 - 01:04

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

Ryne 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.

Log in or register to post comments

Replied on Fri, 09/02/2011 - 20:37

rgore 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.