use isSymmetric(vcov, tol = sqrt(.Machine$double.eps) to check for symmetry of covariance matrices in `parm' (suggested by Fritz) better error messages suggested by Rich