Two recent papers, P. H. Yoon and G. Choe, Phys. Plasmas 28, 082306 (2021) and Yoon et al., Phys. Plasmas 29, 112303 (2022), utilized in the derivation of the kinetic equation for the intensity of turbulent fluctuations the assumption that the wave spectra are isotropic, that is, the ensemble-averaged magnetic field tensorial fluctuation intensity is given by the isotropic diagonal form, <delta B-i delta B-j >(k) = <delta B-2 >(k)delta(ij). However, it is more appropriate to describe the incompressible magnetohydrodynamic turbulence involving shear Alfvenic waves by modeling the turbulence spectrum as being anisotropic. That is, the tensorial fluctuation intensity should be different in diagonal elements across and along the direction of the wave vector, <delta B-i delta B-j >(k) = 1/2 <delta B-perpendicular to(2)>(k) (delta(ij) - k(i)k(j)/k(2)) + <delta B-parallel to(2)>(k) (k(i)k(j)/k(2)). In the present paper, we thus reformulate the weak magnetohydrodynamic turbulence theory under the assumption of anisotropy and work out the form of nonlinear wave kinetic equation. (C) 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/)