Methylamine clusters consisting of up to four molecules were studied by employing Hartree-Fock, density functional theory, and Moller-Plesset calculations with the 6-31+G* basis set. The dimer was found to exhibit two minima with similar interaction energies (-13 kJ/mol) and involving a hydrogen bond. The dipole moment for the dimer differs by up to 20% from the vector addition of the dipole moments for the individual molecules by effect of the interaction; also, the N-H bond distance in the group involved in the hydrogen bond is lengthened by up to 0.006 A as a result. The minima identified for the trimer and tetramer possess cyclic structures that differ in the position of the methyl groups with respect to the plane containing the hydrogen bonds. The contribution of nonadditivity to the interaction in these structures is quite significant (12%-18% of the overall interaction energy). N-H distances in the donor molecule are longer than in the dimer and increase with increasing cluster size. Likewise, the hydrogen bonding energy increases with cluster size. These results expose the significance of cooperative phenomena in the interactions between methylamine molecules. The computations predict slight shifts in the C-N stretching frequencies, which are consistent with their experimental counterparts. The N-H stretching and the NH2 wagging modes undergo large shifts, with a magnitude that increases as the size of the cluster grows.