add exact computation of PR in R