evaluate threshold to get threshold distribution
Adding nonstationarity module