I discuss an analysis to extract the quark and gluon distributions of pion and kaon using meson-induced Drell-Yan and quarkonium production data. Starting from the statistical model first developed for determining the partonic structure of spin-1/2 nucleon, we have extended this approach to perform a global fit to existing meson-induced Drell-Yan and J/ψ production data. These data are well described by the statistical model, allowing a first extraction of the kaon PDFs. We find that both the Drell-Yan and the J/ψ data favor a harder valence distribution for strange quark than for up quark in kaon. The kaon gluon distribution is further constrained by the J/ψ production data. In particular, the momentum fraction carried by gluons is found to be similar for pion and kaon.