This study is an attempt to investigate the applicability of the nonextensive statistics involving the effects of magnetic fields on the thermodynamics of quantum chromodynamics (QCD). The non-extensive statistics are controlled by the entropic index, q ̸= 1. In the case of q = 1, the Boltzmann–Gibbs statistics (extensive) are recovered. The thermodynamics such as pressure, entropy, and magnetization are determined for zero and non-zero magnetic fields. Therefore, the magnetic field is divided into strong, eB = 0.2, 0.3 GeV2, and weak eB = 0.002, 0.003, 0.005 GeV2 magnetic fields. The magnetic field effect is caused by adding a vacuum contribution to the free energy alongside the thermal contribution. The theoretical results are confronted with the lattice results which show overestimation especially at high temperatures and with higher entropic index q. It is concluded that QCD matter is considered to be a para-magnetic matter. Nevertheless, the non-extensive statistics might not be a favorable tool for describing the strongly coupled media responses to the external magnetic fields.