# Groups analysis

5 posts / 0 new
Offline
Joined: 06/23/2014 - 16:15
Groups analysis
Offline
Joined: 07/31/2009 - 15:26
No fixed limit

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.

Offline
Joined: 06/23/2014 - 16:15
Operation in R

Hi,

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.

Offline
Joined: 07/31/2009 - 15:26
Yes

Yes, that's correct. The way I've done this before is like the following:

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)
Offline
Joined: 06/23/2014 - 16:15
Thanks

Thank you for providing an example of the analysis.