new function lavCor() fits unrestricted model