EFA with binary indicators

OpenMx doesn't have any built-in methods for exploratory factor analysis, in part because we defer to the existing EFA methods like factanal() that are already a part of R.

There are several ways to supply user-defined models in OpenMx using things like the algebra objective function, so its possible for an interested user to fit just about whatever they want. It may be simpler to use OpenMx or something else to estimate a polychoric correlation matrix and just run that through factanal(), but that method cannot account for missing data.

The front page of the website has an EFA example.

Path-based RAM models can also use binary, ordinal, and continuous data in any mixture.

So, you should be able to run the model you wish using OpenMx. if you need to do rotation, or factor scores (including with missing data) then there are some extra steps needed.