For an extended sample of Seyfert galaxies we compile from the literature the infrared fluxes in the four IRAS bands, the ground-based small-beam (~ 5-10") fluxes in the standard Q, N, M, L (or L') bands, and the nuclear (non-stellar) estimated fluxes in the JHK bands. We estimate nuclear fluxes in the L band by applying a correction for stellar light. From the statistical study of the infrared colors and luminosities, we derive the typical SEDs of Seyfert 1 and 2 nuclei and the typical differences in luminosities between the two types of objects in the mid- and near-infrared spectral ranges. The observational data appear to severely challenge many models of dusty tori, which hardly account for the shapes of the SEDs and the degree of anisotropy observed in Seyfert galaxies. In particular, at variance with some earlier claims, very thick and compact tori are basically inconsistent with these observational constraints. The most successful models, though having problems in accounting for several details, can fit the major infrared observational data of both Seyfert 1 and Seyfert 2 nuclei with tori which extend up to several hundreds pc and have fairly low optical thickness.