I have paired binary responses from same subjects under 9 different conditions. I conducted Cochran's Q which indicated that there significantly different pairs exist. I want to follow up but NOT with post-hoc comparisons: I am not interested in the 36 possible comparisons of these 9 conditions.
I have selected a smaller subset of pairwise comparisons to analyze using McNemar's test. I can manually apply Bonferroni correction by simply dividing my chosen alpha by the number of comparisons I'm running.
However, I worry that Bonferroni may be punishing me too much and I'm interested in also applying other corrections; perhaps Benjamini-Hochberg. How would I implement these alternative corrections (preferably in R). Packages like rcompanion correct for all pairwise comparisons, but I want to correct for a smaller subset that I will test.