This paper proposes a stochastic modeling of plug-in electric vehicles (PEVs) distribution in power systems, and analyzes the corresponding clustering characteristic. It is essential for power utilities to estimate the PEV charging demand as the penetration level of PEV is expected to increase rapidly in the near future. Although the distribution of PEVs in power systems is the primary factor for estimating the PEV charging demand, the data currently available are statistics related to fuel-driven vehicles and to existing electric demands in power systems. In this paper, we calculate the number of households using electricity at individual ending buses of a power system based on the electric demands. Then, we estimate the number of PEVs per household using the probability density function of PEVs derived from the given statistics about fuel-driven vehicles. Finally, we present the clustering characteristic of the PEV distribution via case studies employing the test systems.