Pǿlya–Aeppli regression model for overdispersed count data