Non-linear Auto-Regressive Models for Cross-Frequency Coupling in Neural Time Series