When multiple outliers occur in linear regression model or the distribution of residuals is not normal, we can use residuals rank as weight function to get some resist estimator.

The ARMA model was used to describe the prior distribution of observed discharge and the ar model was adopted to simulate the likelihood function of forecasting error.