Login

News

path style ACE model

Posted on
Picture of user. tbates Joined: 07/31/2009
Forums
Behavioral Genetics Models
hi,
just working on the 'univACEP.R' script in /failing

It gives this error.. which is obscure to me

> twinACEFit <- mxRun(twinACEModel)
Running twinACE
Error in mxRun(twinACEModel) :
REAL() can only be applied to a 'numeric', not a 'character'

Replied on Sat, 09/26/2009 - 12:31
Picture of user. tbates Joined: 07/31/2009

ok, so that error is from a numeric converted to a factor... probably clear enough

if a supermodel doesn't recognise paths, but its submodels do, is there a way to inject new paths into the submodels? Tried this without sucess:

twinAEModel <- mxModel(twinACEModel, name="twinAE",
+ mxPath(from=c("MZ.C1"), to="MZ.bmi1", arrows=1, free=FALSE, values=0, label="c"),
+ mxPath(from=c("DZ.C1"), to="DZ.bmi1", arrows=1, free=FALSE, values=0, label="c")
+ )
Error: The model type of model 'twinACE' does not recognize paths.
+

Log in or register to post comments
Replied on Sat, 09/26/2009 - 14:32
Picture of user. sboker Joined: 07/28/2009

Interesting problem. I agree that manipulating paths within RAM-type submodels would be good to implement. Michael might be able to tell us how difficult that might be. I'm not sure if he is keeping the submodel's RAM type or converting based on the parent.
Log in or register to post comments
Replied on Sat, 09/26/2009 - 14:54
Picture of user. mspiegel Joined: 07/31/2009

Use the following function. Otherwise the implication is that matrix-style models need to be able to process paths. Which really implies that any-style model needs to be able to process weird entities that are going to be used by any other style model. Which breaks all our efforts at encapsulation. Plus there's the problem with misusing the "." notation.

updateSubmodel <- function(model, submodelname, ...) {
   if (model@name == submodelname) {
       return(mxModel(model, ...))
   }
   model@submodels <- lapply(model@submodels, updateSubmodel, submodelname, ...)
   return(model)
}
Log in or register to post comments
Replied on Sat, 09/26/2009 - 14:40
Picture of user. mspiegel Joined: 07/31/2009

No, C1 is a latent variable and bmi1 is a manifest variable, both of which are not named entities. The "." notation currently only works with named entities, which is why it can't be used in the 'from' or 'to' arguments of mxPath, because these argument do not accept named entities.
Log in or register to post comments
Replied on Sat, 09/26/2009 - 14:42
Picture of user. mspiegel Joined: 07/31/2009

"REAL() can only be applied to a 'numeric', not a 'character'" is a bug in the backend. Can you check in a version of the script that is failing with this error to models/failing? I can't reproduce this error with the current version of the script.
Log in or register to post comments
Replied on Sat, 09/26/2009 - 21:27
Picture of user. mspiegel Joined: 07/31/2009

Thanks that was useful. Now we have error checking if the data is a matrix that is not of type 'double'. I will add error checking for data.frames after I consult with Tim Brick to confirm my suspicions on what to do about ordinal data.
Log in or register to post comments
  Login

News