powered by
This function prevents the support of beta becoming greater than n_l/2. This heuristic approach prevents erratic values of BIC when in high-dimensions.
k.index(matrix, n_l)