I've changed the functions 'analyze_covariate_lists', 'checkCovariateNAs', 'check_covariates', and 'printNACovariates' so that analyze calls check_covariates. After check_covariates verifies the presence/absence of a phenotype, it calls printNACovariates to check for NAs within the present phenotype.