Performing time series decomposition using the stl() function