A recent paper walks through what appears to be a similar analysis although possibly with a different methodology. Its plots and tables appear to have been produced in R; however the authors do not provide their code.
In principle, Equation 10 can be implemented in R; doing so is not trivial, however, and may require a sophisticated understanding of several different statistical tools.
Without a citation to the paper shown, it's not possible to make any more concrete suggestion than to run searches on rseek, an R-tuned Google front-end portal against the author(s)'s names and keywords. Someone may already have implemented the algorithm in R.