Kinetic Monte Carlo (KMC) methods have been a successful technique for accelerating time scales and increasing system sizes beyond those achievable with fully atomistic simulations. However, a requirement for its success is a priori knowledge of all relevant reaction pathways and their rate coefficients. This can be difficult for systems with complex chemistry, such as shock-compressed materials at high temperatures and pressures or phenolic spacecraft heat shields undergoing pyrolysis, which can consist of hundreds of molecular species and thousands of distinct reactions. In this work, we develop a method for first estimating a KMC model composed of elementary reactions and rate coefficients by using large datasets derived from a few molecular dynamics (MD) simulations of shock compressed liquid methane, and then using L1 regularization to reduce the estimated chemical reaction network. We find that the full network of 2613 reactions can be reduced by 89% while incurring approximately 9% error in the dominant species (CH4) population. We find that the degree of sparsity achievable decreases when similar accuracy is required for additional populations of species.