uci_matrix function

Uniform credible intervals in matrix-valued case