powered by
estimate the variance and covariance parameters
estVarCovarCoef( data, g, eqns, pt, ordered, vc.coef.estimator = "ML", vc.coef.missing = "FIML" )