We report global analysis results of experimental data for nuclear structure-function ratios F_2^A/F_2^{A'} and proton-nucleus Drell-Yan cross-section ratios sigma_{DY}^{pA}/sigma_{DY}^{pA'} in order to determine optimum parton distribution functions (PDFs) in nuclei. An important point of this analysis is to show uncertainties of the distributions by the Hessian method. The results indicate that the uncertainties are large for gluon distributions in the whole x region and for antiquark distributions at x>0.2. We provide a code for calculating any nuclear PDFs at given x and Q^2 for general users. They can be used for calculating high-energy nuclear reactions including neutrino-nucleus interactions, which are discussed at this workshop.