[en] The neutron capture gamma-ray spectra for 8 nuclides, ⁸⁹Y, ⁹³Nb, ¹²⁷I, ¹³³Cs, ¹⁴¹Pr, ¹⁹⁷Au, ^natTl, and ²⁰⁹Bi, were calculated by using the Hauser-Feshbach statistical model, and their results were compared with the available experimental data. Two dominant ingredients to perform the statistical calculation were the level densities described by the Gilbert-Cameron approach with an improved systematics, and the gamma-ray transmission coefficients described by gamma-ray strength functions. Although various gamma-ray strength functions with a Lorentzian formula have been developed by using the photonuclear data or a microscopic analysis, they have failed to reproduce the occasional anomalous bumps observed near or below a neutron binding. In this work, we could reproduce the bumps well by adding a Lorentzian with an energy-temperature dependent width into a giant dipole resonance with an enhanced generalized Lorentzian. In addition, we introduced a correction function so as to compensate for the shortcomings due to missing levels or level-cuts. (author)