Higher-order correlations in non-stationary parallel spike trains: statistical modeling and inference