We have developed a method for analyzing the (hyper)polarizabilities of open-shell molecular systems. This method employs the (hyper)polarizability densities based on the natural orbitals and occupation numbers, which enables us to analyze the contributions of odd electrons having various open-shell (diradical) characters. Within broken-symmetry, i. e., spin-unrestricted, single-determinant molecular orbital and density functional theory approaches, we can also remove the spin contamination effects on these quantities through spin projection. To do that, an approximate spin projected method has been elaborated and applied to the analysis of the (hyper)polarizability of multi-radical systems. As examples, typical open-shell singlet systems, 1,3-dipoles and rectangular graphene nanoflakes, are examined.