[en] A global analysis of nuclear medium modifications of parton distributions is presented using deeply inelastic scattering data of various nuclear targets. Two obtained data sets are provided for quark and gluon nuclear modification factors, referred as nIMParton16. One is from the global fit only to the experimental data of isospin-scalar nuclei (Set A), and the other is from the fit to all the measured nuclear data (Set B). The scale-dependence is described by DGLAP equations with nonlinear corrections in this work. The Fermi motion and off-shell effect, nucleon swelling, and parton–parton recombination are taken into account together for modeling the complicated x-dependence of nuclear modification. The nuclear gluon shadowing in this paper is dynamically generated by the QCD evolution of parton splitting and recombination processes with zero gluon density at the input scale. Sophisticated nuclear dependence of nuclear medium effects is studied with only two free parameters. With the obtained free parameters from the global analysis, the nuclear modifications of parton distribution functions of unmeasured nuclei can be predicted in our model. Nuclear modification of deuteron is also predicted and shown with recent measurement at JLab.