Multivariate statistical methods are widely used in several disciplines of fundamental sciences. In the present study, the data analysis of the chemical analysis of the sands of Moonlight Beach in the Kemer region was examined using multivariate statistical methods. This study consists of three parts. The multivariate statistical analysis tests were described in the first part, then the pollution indexes were studied in the second part. Finally, the distribution maps of the chemical analyses and pollution indexes were generated using the obtained data. The heavy metals were mostly observed in location K1, while they were sorted as follows based on their concentrations: Mg > Fe > Al > Ti > Sr > Mn > Cr > Ni > Zn > Zr > Cu > Rb. Also, strong positive correlations were found between Si, Fe, Al, K, Ti, P. According to the results of factor analysis, it was found that four factors explained 83.5% of the total variance. On the other hand, the coefficient of determination (R-2) was calculated as 63.6% in the regression model. Each unit increase in the value of Ti leads to an increase of 0.022 units in the value of Si. Potential Ecological Risk Index analysis results (RI < 150) revealed that the study area had no risk. However, the locations around Moonlight Beach are under risk in terms of Enrichment Factor and Contamination Factor values. The index values of heavy metals in the anomaly maps and their densities were found to be successful; and higher densities were observed based on heavy metal anomalies.