The ongoing study aims to establish a direct probabilistic load flow (PLF) for the analysis of wind integrated radial distribution systems. Because of the stochastic output power of wind farms, it is very important to find a method which can reduce the calculation burden significantly, without having compromising the accuracy of results. In the proposed approach, a K-means based data clustering algorithm is employed, in which all data points are bunched into desired clusters. In this regard, probable agents are selected to run the PLF algorithm. The clustered data are used to employ the Monte Carlo simulation (MCS) method. In this paper, the analysis is performed in terms of simulation run-time. Also, this research follows a two-fold aim. In the first stage, the superiority of data clustering-based MCS over the unsorted data MCS is demonstrated properly. Moreover, the impact of data clustering-based MCS and unsorted data-based MCS is investigated using an indirect probabilistic forward/backward sweep (PFBS) method. Thus, in the second stage, the simulation run-time comparison is carried out rigorously between the proposed direct PLF and the indirect PFBS method to examine the computational burden effects. Simulation results are exhibited on the IEEE 33-bus and 69-bus radial distribution systems.