Clustered factor analysis of multineuronal spike data