stepBF function

Identify the optimal number of shifts using Bayes factors