FIML for multilevel modeling in R?

[Bonnie Dixon]

Does anyone know of a way of doing full information maximum likelihood (FIML) estimation in R in the context of a mixed effects model?  This method is also sometimes referred to as "direct maximum likelihood" or "raw maximum likelihood".  It is a method that can handle missing data without the need for imputation and is commonly used in structural equation modeling (such as in the R package, lavaan).  But, after much searching, I cannot find an implementation of this method for multilevel modeling in R.  If anyone knows of a way of doing this, please let me know.  Thanks.

Bonnie