The Aegean region is geologically situated at the western end of the Gediz Graben system, influenced by the Western Anatolian Regime. In addition, the region is characterized by various active fault lines that can generate earthquake activity. Numerous earthquakes have been recorded in the region, causing significant material and moral damage from the past to the present. In this study, earthquake data from three different catalogs are examined. The non-clustered catalog is compiled for the years 1970 to 2020, including earthquakes with a moment magnitude (Mw) greater than 3.0. The monthly average magnitudes of earthquakes in the region are obtained and analyzed using ARIMA, singular spectrum analysis (SSA), and deep learning methods including convolutional neural network (CNN) and long short-term memory (LSTM), as these methods have not been compared for the region previously. Each method has a different benefit. ARIMA analyzes time series trends and seasonal patterns, while SSA focuses on decomposition and feature extraction. LSTM attempts to capture complex relationships using memory mechanisms, while CNN is powerful at pattern recognition and extracting important features. Thanks to this diversity, our study allows for more comprehensive and reliable forecasts of average earthquake magnitudes for the next 36 periods. The estimation capabilities and error rates of each method were analyzed based on earthquake magnitude data, and it was determined that the LSTM method provided the most effective and accurate predictions.