Vibrational energy levels, wave functions, and ir absorption intensities are computed for (H2O)(n) clusters with n = 2, 3, 4, and 5. The calculations were carried out by the vibrational self-consistent field (VSCF) approximation, with corrections for correlations between the modes by perturbation theory. This correlation corrected VSCF (CC-VSCF) is analogous to the familiar Moller-Plesset method in electronic structure theory. Test calculations indicate that this method is of very good accuracy also for very anharmonic systems. While the method is of highest relative accuracy for the stiffest modes, it works very well also for the soft ones. Some of the main results are (1) the frequencies calculated are in good but incomplete agreement with experimental data available for some of the intramolecular mode excitations. The deviations are attributed to the inaccuracy of the coupling between intramolecular and intermolecular modes for the potential function used. (2) Insight is gained into the pattern of blue- or redshifts from the corresponding harmonic excitation energies for the various modes. (3) Anharmonic coupling between the modes dominates in general over the intrinsic anharmonicity of individual modes in determining the spectrum. (4) The anharmonic corrections to the frequencies of some intermolecular modes (shearing, torsional) are extremely large, and exceed 100% or more in many cases. (5) An approximation of quartic potential field in the normal mode displacement is tested for the clusters. It works well for the high and intermediate frequency modes, but is in error for very soft shearing and torsional modes. (6) The relative errors of the VSCF approximation are found to decrease with the cluster size. This is extremely encouraging for calculations of large clusters, since the VSCF level is computationally simple.