The problem of state estimation for noisy discrete-time nonlinear dynamic systems is discussed by proposing the multi-stage iteration filter-smoother in conjunction with the algorithms for prediction, filtering and smoothing which are based upon the viewpoint of marginal maximum likelihood estimation. Results of the numerical example indicate greatly improved preformance over the filtering only.