  We present an analysis to extract kaon parton distribution functions (PDFs)
for the first time 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 and later applied to the spin-0 pion, we
have extended this approach to perform a global fit to existing kaon-induced
Drell-Yan and $J/\psi$ 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/\psi$ 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/\psi$ production data. In particular, the
momentum fraction carried by gluons is found to be similar for pion and kaon.
