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)a list including all possible combinations of the triplet
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