The non-linear dynamics of a chain of stationary, saturated magnetic islands is studied by solving a four-field system of equations that include non-ideal effects, lowest order finite Larmor radius corrections and neoclassical terms. The magnetic island rotation velocity is calculated self-consistently with the fields profiles. The solutions for the island rotation velocity and for the ion polarization current are determined as a function of the characteristic parameters of the system and the results are discussed. The results of the calculations show that island rotation velocity and the ion polarization current depend in a non-trivial way on the parameters characterizing the system, and some clear patterns emerge only in particular cases. An analysis of magnetic island rotation velocity is performed on experiments in COMPASS tokamak. Measured island rotation velocity is compared with the calculated ion and electron flow velocities, for different hypotheses on the toroidal rotation of the plasma. The comparison shows that the island rotation velocity is consistent with the ion flow velocity, under the hypothesis of slow toroidal rotation and low collisionality. Theoretical calculation of the island rotation velocity according to the model here developed suggests that the islands rotate weakly in the ion direction, in the hypothesis of slow toroidal rotation and high collisionality. The impossibility of directly measuring the plasma rotation velocity makes it difficult to distinguish between these different regimes.