Abstarct: The current thesis focused on SHAMKAN stream sediment lithogeochemistry investigation. The study area is located on the north-east of Iran in Khorasan-Razavi Province. During the exploration program, 398 stream sediments (30 duplicate samples) and 131 heavy mineral samples were taken from the study area. Size fraction of -230 mesh of stream sediment and -20 mesh of heavy minerals were separated and sent for analysis. The stream sediment samples were analyzed by ICP-OES for 54 elements in the geology society of Iran-Tehran lab. Data preparation was the first step to confirm the data set. In the next step, single, two and multi – variant analyses were done on data set. Elements such as Cr, B, Pb, Cu, and Au were introduced as first priority baxsed on the single – variant analysis and secondary dispersion haloes. According to the statistical analyses, Sb, Al, V, Tl, Zn, Ni, Mg, Mn were considered as second priority and the elements with the third priority were As and S. According to the correlation coefficient, Au was shown a good correlation with V, Bi , Na , Cu , Ge , Cd , Hf , Sc , U. The geochemical dispersion hole maps of the corresponding elements were illustrated a good overlap in the center , north and NW of the study area. Therefore, these elements were considered as first priority. There were a very good correlation coefficients between U- Hf-Ge with B, Al- V- Sc- Ti- Na- Ga- Bi with Cu, and V- Cu- Al with Na. Cluster analysis and principal component analysis were applied and identified probable mineralization of Au, Pb, Zn, and Cu. The relationship of Cu, V, Mg, Co, Cr, Ni and Au in a cluster was interpreted as of hydrothermal fluids effects on ultramafic, mafic to intermediate rocks. The gold anomaly area close to the outcrops of granitoid rocks might be assumed as relationship of these intrusive rocks and mineralization. The cluster analysis was grouped Pb togather with effects of sedimentary rocks. is result was confirmed by the heavy mineral geochemistry investigations. Gold was used as high priority element for the next statistical analysis. The probability plot modeling method was used to separate the background from transition and anomaly data. Three codes of 0, 1 and 2 were used for background, transition and anomaly data and discriminant function modeling, DFM. The general trend of DFM was used to evaluate the correlation of all elements with their sub population. The stepwise method was applied to determine the most important elements and their correlation with their sub-population. The general trend could evaluate 92% of relation baxsed on all elements and stepwise could recognize 85.4% of variations baxsed on seven elements. The seven elements were introduced for the follow up exploration program use. In the current research, Au was introduced as target for DFM. In the case of polymextal mineralization, DFM technique needs to use for all of the polymextals.