Monte Carlo simulations of some typical order-disorder ferroelectrics such as TGS, NaNO2 and DKDP nanocrystals were studied using a Transverse Ising Model Hamiltonian with four-spins interactions. The microscopic parameters corresponding to this Hamiltonian were adjusted to fit the experimental polarization-temperature curves for each one of the materials in the bulk phase. Then the dependences of the ferroelectric-paraelectric phase transition temperatures, Tc, on the sizes of those crystals were studied with Monte Carlo simulations of the order-disorder system. We report a weak dependence of Tc on the size of the crystal (d) for these materials above d∼6nm. The addition of surface effects showed that the expected lowtemperature shift of Tc due to size effects, can be reverted.