I am running beast on time series as discussed here Detect changes in the seasonal component using bfast . My code is working and plotting but I am rather getting the follow