A method is presented for the estimation of the parameters of a noncausal nonminimum phase ARMA model for non-Gaussian random processes. Using certain higher-order cepstra slices, the Fourier phases of two intermediate sequences, hmin(n) and hmax(n), can be computed, where hmin(n) is composed of the minimum phase parts of the AR and MA models, and hmax(n) of the corresponding maximum phase parts. Under the condition that the AR and MA models do not have common zeros, these two sequences can be estimated from their phases only, and lead to the reconstruction of the AR and MA parameters, within a scalar and a time shift. The AR and MA orders do not have to be estimated separately, but they are a byproduct of the parameter estimation procedure. Unlike existing methods, the estimation procedure is fairly robust if a small order mismatch occurs. Since the robustness of the method in the presence of additive noise depends on the accuracy of the estimated phases of hmin(n) and hmax(n), the phase errors occurring due to finite length data are studied.