A two-stage least squares based iterative parameter estimation algorithm is proposed for identifying a feedback nonlinear system with the open-loop being a controlled autoregressive moving average model from input-output data. The identification model is bilinear on two unknown parameter vectors. By