Incorporating neural response variability into models for neural coding