Negative binomial regression for spike trains using Polya-gamma augmentation