p.adjust in cuminc_pairs is not working. method="BH", method="bonferroni", and method="none" give equal results.