There are various treatment modalities for prostate cancer, which has a high incidence. In this study, it is aimed to make predictions with machine learning in order to determine the optimal treatment option for prostate cancer patients. The study included 88 male patients diagnosed with prostate cancer. Independent variables were determined as Gleason scores, biopsy, PSA, SUV_max, and age. Prostate cancer treatments, which are dependent variables, were determined as hormone therapy(n = 30), radiotherapy(n = 28) and radiotherapy + hormone therapy(n = 30). Machine learning was carried out in the Python with SVM, RF, DT, ETC and XGBoost. Metrics such as accuracy, ROC curve, and AUC were used to evaluate the performance of multi-class predictions. The model with the highest number of successful predictions was the XGBoost. False negative rates for hormone therapy, radiotherapy, and radiotherapy + hormone therapy treatments were, respectively, 12.5, 33.3, and 0%. The accuracy values were computed as 0.61, 0.83, 0.83, 0.72 and 0.89 for SVM, RF, DT, ETC and XGBoost, respectively. The three features that had the greatest influence on the treatment model prediction for prostate cancer with XGBoost were biopsy, Gleason score (3 + 3), and PSA level, respectively. According to the AUC, ROC and accuracy, it was determined that the XGBoost was the model that made the best estimation of prostate cancer treatment. Among the variables biopsy, Gleason score, and PSA level are identified as key variables in prediction of treatment.