Statistics of Model reduction and frequency residuals for a robust estimation of nonlinearities in subspace identification