Reactor neutrino experiments play a crucial role in advancing our knowledge of neutrinos. In this Letter, the evolution of the flux and spectrum as a function of the reactor isotopic content is reported in terms of the inverse-beta-decay yield at Daya Bay with 1958 days of data and improved systematic uncertainties. These measurements are compared with two signature model predictions: the Huber-Mueller model based on the conversion method and the SM2018 model based on the summation method. The measured average flux and spectrum, as well as the flux evolution with the ^{239}Pu isotopic fraction, are inconsistent with the predictions of the Huber-Mueller model. In contrast, the SM2018 model is shown to agree with the average flux and its evolution but fails to describe the energy spectrum. Altering the predicted inverse-beta-decay spectrum from ^{239}Pu fission does not improve the agreement with the measurement for either model. The models can be brought into better agreement with the measurements if either the predicted spectrum due to ^{235}U fission is changed or the predicted ^{235}U, ^{238}U, ^{239}Pu, and ^{241}Pu spectra are changed in equal measure.
