A function to find the optimal combination of first two splits, and the corresponding Q-between given the first split
find_second_split(xk, first.splits, y, vi, minbucket, minsplit)
moderator vector
possible first splits
effect size
sampling variance
the minimal number of studies in child nodes
the minimal number of studies in parent node
a list including all possible combinations of the triplet