new r code for pairwise treatment comparisons