“…For p = 2, we use the simulation configurations from Table 1 of Paul [7] with m 1 = 27, n i = 12, 7, 6, 6,7,8,10,7,8,6,11,7,8,9,2,7,9,7,11,10,4,8,10,12,8,7,8, and m 2 = 19, n i = 5, 11, 7, 9, 12, 8, 6, 7, 6, 4, 6, 9, 6, 7, 5, 9, 2, 6, 9. For p = 3, we use m 1 = 27, m 2 = 19 with the same cluster sizes as those for p = 2. In addition, we add m 3 = 32 and join the sizes of the control and the treatment groups from the data in Section 3 of Williams [13] to form n i = 13, 12, 9, 9, 8, 8, 13, 12, 10, 10, 9, 13, 5, 7, 10, 10, 12, 11, 10, 9, 11, 10, 10, 9, 9, 5, 9, 7, 10, 6, 10, 7. Paul and Islam [8] parameterized the correlation ρ as θ/(1 + θ) to express the beta-binomial distribution in a more compact form.…”