Techniques are described for forecasting and dynamic routing of incoming routing service requests in a service environment that includes one or more call centers. Implementations apply a doubly stochastic modeling technique to modeling call volumes, call wait times, calling handling times, and/or other parameters in a service environment that includes one or more call centers that include multiple skill codes for the service representatives. By simultaneously modeling on different time scales, such as modeling both inter-day and intra-day correlations, implementations are able to provide an analysis that exploits more of the structure in the data compared to traditional techniques. Implementations can also integrate in other types of time series for other effects that may impact the call volume or other parameters in a service environment, such as exogenous and/or anomalous variables that are independent of other inter-day and/or intra-day trends exhibited in the data.