A microscopic (Hamiltonian-based) method for the quantum statistics of bosonic excitations in a two-mode magnon system is developed. Both the exchange and the dipole–dipole interactions, as well as the Zeeman term for an external applied field, are included in the spin Hamiltonian, and the model also contains the nonlinear effects due to parallel pumping and four-magnon interactions. The quantization of spin operators is achieved through the Holstein–Primakoff formalism, and then a coherent magnon state representation is used to study the occupation magnon number and the quantum statistical behaviour of the system. Particular attention is given to the cross correlation between the two coupled magnon modes in a ferromagnetic nanowire geometry formed by two lines of spins. Manipulation of the collapse-and-revival phenomena for the temporal evolution of the magnon number as well as the control of the cross correlation between the two magnon modes is demonstrated by tuning the parallel pumping field amplitude. The role of the four-magnon interactions is particularly interesting and leads to anti-correlation in some cases with coherent states.