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.