[R] easy way to fit saturated model in sem package?

yrosseel yrosseel at gmail.com
Mon Jul 23 10:41:51 CEST 2012

> I will check out the lavaan package.

Dear Joshua,

The lavaan package may help you. The FIML estimator typically starts 
with the EM algorithm to estimate the moments of the unrestricted model. 
There is no 'one-shot' function for it, at the moment, but if you only 
need those moments, you can do something like this:

Suppose your data is a data.frame called 'HS.missing', then the 
following commands can be used to get the estimated moments:


Data <- lavaan:::lavData(HS.missing, ov.names=names(HS.missing), 

# we assume only 1 group
Missing <- lavaan:::getMissingPatternStats(X = Data at X[[1L]], Mp = 
Data at Mp[[1L]])

# compute moments using EM algorithm
Moments <- lavaan:::estimate.moments.EM(X=Data at X[[1L]], M=Missing, 

# estimated covariance matrix

# estimated mean vector

Hope this helps,

Yves Rosseel

More information about the R-help mailing list