Categorical Outcomes

Hello everyone!
In the process of learning a bit more about how OpenMX works and getting up to speed on twin study methodology, I have been running some of the presentation scripts from the Boulder conference's intro seminar. I'm having a bit of a problem with one of the scripts for binary outcome twin models and was hoping that someone recognized the issue. The script I am using is here: http://ibg.colorado.edu/cdrom2010/rijsdijk/OrdinalTwinAnalysis/Binary/UnivACE_MatrRawOrd.R

Hi
I'm trying to get an ordinal ACE model working but it keeps crashing R.
I'm using OpenMx .5 R 2.10.1 under Windows 7 64 bit

Any suggestions would be greatly appreciated
thanks
s

Hi,

Why is the mxConstraint command not working anymore?
How else is one to constrain the variances to unity?

Thanks.

Fruhling

# Algebra to compute total variances and standard deviations (diagonal only)
mxAlgebra( expression=A+C+E, name="V" ),
mxMatrix( type="Iden", nrow=nv, ncol=nv, name="I"),
mxAlgebra( expression=solve(sqrt(I*V)), name="sd"),
# Constraint on variance of ordinal variables
mxConstraint( alg1="V", "=", alg2="I", name="Var1"),

Error in mxConstraint(alg1 = "V", "=", alg2 = "I", name = "Var1") :

The mxFactor() function was introduced in the 0.3.0 release. In the source code repository, it has been updated to accept multiple columns of a data frame. Here is an example of the new functionality:

nthresh1 <- 1
nthresh2 <- 12
data <- read.table("data/mddndzf.dat", na.string=".",
col.names=c("t1neur1", "t1mddd4l", "t2neur1", "t2mddd4l"))
data[,c(1,3)] <- mxFactor(data[,c(1,3)], c(0 : nthresh2))
data[,c(2,4)] <- mxFactor(data[,c(2,4)], c(0 : nthresh1))