powered by
Uniform credible intervals in matrix-valued case
uci_matrix(fpsd.sample, alpha, uniform_among_components = F)
see (6.5) in Meier (2018)