  Dwarf spheroidal galaxies (dSphs) are promising targets of indirect detection
experiments searching for dark matter (DM) at present universe. Toward robust
prediction for the amount of signal flux originating in DM annihilation inside
dSphs, a precise determination of DM distributions as well as $J$-factors of
the dSphs is particularly important. In this work, we estimate those of Draco,
Sculptor, and Ursa Minor dSphs by an improved statistical method in which both
foreground stars and dSph member stars are simultaneously taken into account.
We define the likelihood function of the method as the so-called conditional
one to remove sampling bias of observed stellar data. This improved method
enables us to estimate DM distributions and $J$-factors of the dSphs directly
from observed stellar data contaminated by foreground stars without imposing
stringent membership criteria on the measured quantities.
