You are here

Is it possible to estimate an AE model using the ACE model functions?

2 posts / 0 new
Last post
tbates's picture
Offline
Joined: 07/31/2009 - 14:25
Is it possible to estimate an AE model using the ACE model functions?

Simple question: How do I estimate an AE model using the umxACE function?

tbates's picture
Offline
Joined: 07/31/2009 - 14:25
umxModify

Yes, you can estimate an "AE" umxACE model!

This is done via umxModify and labels - the labels umx applies to every matrix cell/parameter in the model.

in twin models, these labels take the form matrix name + "_" + "r" + row-number + "c" + col-number. So the link from the first latent variable to the second measured variable in the c matrix is labeled c_r2c1.

Just run umxACE to get your saturated model, then drop the c parameters using umxModify

There is an example of this in ?umxACE

This is a powerful tool. You can provide a single label to update. e.g. update = "c_r1c1", this can be a list of labels.
By default they are set to zero, but you can set them to anything via values. By default they are fixed, but you can free variables with free = T.

You can even use regular expressions: umxModify(m1, regex = "c_r.c.”) which would drop all free parameters from the c matrix (as long as there are fewer than 10 rows and columns), which is what you want to create an AE model!

then umxCompare to compare them

Because I am lazy, you can update the model , rename it, and do the comparison all in 1 line!

m2 = umxModify(m1, regex = "c_r.c.”, name = "AE", comp= TRUE)