A randomised approach for NARX model identification based on a multivariate Bernoulli distribution