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:
library(lavaan)
Data <- lavaan:::lavData(HS.missing, ov.names=names(HS.missing),
missing="fiml")
# we assume only 1 group
Missing <- lavaan:::getMissingPatternStats(X = Data@X[[1L]], Mp =
Data@Mp[[1L]])
# compute moments using EM algorithm
Moments <- lavaan:::estimate.moments.EM(X=Data@X[[1L]], M=Missing,
verbose=TRUE)
# estimated covariance matrix
Moments$sigma
# estimated mean vector
Moments$mu
Hope this helps,
Yves Rosseel
http://lavaan.org
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
