We study an extension of the stochastic SIS model generated under the assumption that the disease transmission coefficient in an infinitesimal time interval is \(\beta (t) + \sigma (t)\xi _t\) , where \(\beta (t)\) and \(\sigma (t)\) are bounded time-dependent functions and \(\xi _t\) is a white noise. We verify that there is a global positive solution to the stochastic differential equation (SIS-Itô-SDE) resultant for the infected and give the necessary conditions in their coefficients to achieve endemic persistence or disease extinction. Also, a numerical procedure to estimate the parameters of the respective SDE is proposed to make forecasts useful in the initial stages of an epidemic. The method begins with the discretization of infectives; our approach assumes that the SIS-Itô-SDE involves unknown transmission \(\beta (t)\) and diffusion \(\sigma (t)\) rates, which we estimate. Specifically, we employ the maximum likelihood estimator (MLE) to gauge the transmission rate’s value, ensuring its strong consistency. For estimating the diffusion parameter, we utilize the quadratic variation of the data. To validate our methodology, we simulate using various datasets. Subsequently, using the empirical distributions of these same estimators and the Milstein method, the number of infected is predicted. To evaluate the effectiveness of this method, by using the RMSE metric, the forecast accuracy of the number of infected individuals reported in two studies is compared with real historical data.