GLM for an ANCOVA with count data and multiple Covariates