Flexible joint modeling of mean and dispersion for the directional tuning of neuronal spike counts