This is a nice question! I've done multigroup analyses with 800 groups. There is no fixed limit on the number of groups. At some point you will start to fill up R's protection pointer stack. We catch this in OpenMx, and tell you to let us know about it on the these forums if it happens.
If I understood it correctly, I could just specify e.g. 100 groups in R, (perhaps with equality constraints), in one analysis, without doing something extraordinary?
No fixed limit
Log in or register to post comments
Operation in R
If I understood it correctly, I could just specify e.g. 100 groups in R, (perhaps with equality constraints), in one analysis, without doing something extraordinary?
Thanks.
Log in or register to post comments
In reply to Operation in R by metavid
Yes
grpmodels <- list()
dataL <- list()
for(k in modelDataIndexes){
dataL[[k]] <- subsetOfYourBigDataBasedOn[k]
grpmodels[[k]] <- mxModel(name=paste('grp', k, sep=''),
model stuff, might depend on k,
paths, matrices, algebras,
expectation, fitfunction,
mxData(dataL[[k]], type='raw')
) #close model
}
groupfits <- paste("grp", modelDataIndexes, ".fitfunction", sep="")
cmod <- mxModel(model='Container', grpmodels, mxFitFunctionMultigroup(groupfits))
cfit <- mxRun(cmod)
Log in or register to post comments
In reply to Yes by mhunter
Thanks
Log in or register to post comments