Packaging during the passaging of viruses in cell cultures yields various phenotypes and is regulated by viral protein expression in infected cells. Although such a packaging mechanism has a profound effect in controlling the virus yield, little is known about the underlying statistical models followed by virus packaging and protein expression among cells infected with the virus. A predictive framework combining identification of the probability density function (PDF) based on log-likelihood and using the PDF for Monte-Carlo simulations is developed. The Birnbaum–Saunders distribution was found to be consistent with all three-virus packaging levels, including nucleocapsids/occlusion-derived virus (ODV), ODVs/polyhedra, and polyhedra/cell for both wild-type and genetically modified AcMNPV. Next, it was demonstrated that PDF fitting could be used to compare two viruses having distinctly different genetic configurations. Finally, the identified PDF can be incorporated in RNA synthesis parameters for baculovirus infection to predict the cell-to-cell variability in protein expression using Monte-Carlo simulations. The proposed tool can be used for the estimation of uncertainty in the kinetic parameter and prediction of cell-to-cell variability for other biological systems. © 2020 Wiley Periodicals, Inc.