Convolutional spike-triggered covariance analysis for neural subunit models