revise ttestGC() so that formula version will not fail when factor variable has NAs