compute_post_hoc_test_with_limma function

Compute pairwise post-hoc comparisons for a multi-level factor