The increased production and price of rare earth elements(REEs) are indicative of their importance and of growing global attention. More accurate and practical exploration procedures are needed for REEs, and for other geochemical resources. One such procedure is a multivariate approach. In this study, five classifiers, including multilayer perceptron(MLP), Bayesian, k-Nearest Neighbors(KNN), Parzen, and support vector machine(SVM),were applied in supervised pattern classification of bulk geochemical samples based on REEs, P, and Fe in the Kiruna type magnetite-apatite deposit of Se-Chahun,Central Iran. This deposit is composed of four rock types:(1) High anomaly(phosphorus iron ore),(2) Low anomaly(metasomatized tuff),(3) Low anomaly(iron ore), and(4)Background(iron ore and others). The proposed methods help to predict the proper classes for new samples from the study area without the need for costly and time-consuming additional studies. In addition, this paper provides a performance comparison of the five models. Results show that all five classifiers have appropriate and acceptable performance. Therefore, pattern classification can be used for evaluation of REE distribution. However, MLP and KNN classifiers show the same results and have the highest CCRs in comparison to Bayesian, Parzen, and SVM classifiers. MLP is more generalizable than KNN and seems to be an applicable approach for classification and predictionof the classes. We hope the predictability of the proposed methods will encourage geochemists to expand the use of numerical models in future work. The increased production and price of rare earth elements (REEs) are indicative of their importance and of growing global attention. More this and their exploration of resources are needed for REEs, and for other geochemical resources. One such procedure is a multivariate approach. study, five classifiers, including multilayer perceptron (MLP), Bayesian, k-Nearest Neighbors (KNN), Parzen, and support vector machine (SVM), were applied in supervised pattern classification of bulk geochemical samples based on REEs, P, and Fe in the Kiruna type magnetite-apatite deposit of Se-Chahun, Central Iran. This deposit is composed of four rock types: (1) High anomaly (phosphorus iron ore), (2) Low anomaly (metasomatized tuff), (3) Low anomaly (iron ore), and (4) Background (iron ore and others). The proposed methods help to predict the proper classes for new samples from the study area without the need for costly and time-consuming additional studies. In addition, this paper provides a performance Comparison of the five models. Therefore, pattern classification can be used for evaluation of REE distribution. However, MLP and KNN classifiers show the same results and have the highest CCRs in comparison to Bayesian , Parzen, and SVM classifiers. MLP is more generalizable than KNN and seems to be applicable applicable for classification and prediction of the classes. We hope the predictability of the proposed methods will encourage geochemists to expand the use of numerical models in future work.