The pharmaceutical industry is facing challenges due to various factors such as supply chain disruptions, changing consumer behavior, and regulatory changes. Accurate demand forecasting is essential to ensure an adequate supply of drugs. The goal of this work is to forecast paracetamol product demand. For this purpose, we propose a hybrid forecasting model combining two effective forecasting techniques: SARIMA (Seasonal AutoRegressive Integrated Moving Average) and ANFIS (Adaptive Neuro-Fuzzy Inference System). This proposal consists of nonlinear components of time series by ANFIS and adjusting the result by the mean of the residuals of the SARIMA to improve the accuracy and performance of ANFIS predictions. Before the prediction phase, we preprocess our data and detect the anomalies in our dataset with Locally Selective Combination in Parallel Outlier Ensembles (LSCP). Then, by treating these anomalies as missing values, they are imputed using the combination of Fuzzy-Possibilistic c-means (FCM) with support vector regression (SVR) and a genetic algorithm (GA). Finally, we evaluate the performance of the model and some known models based on MAPE. We choose the hybrid model SARIMA-ANFIS that provides the most accurate and reliable forecasting.