pip function

Calculates posterior inclusion probabilities (PIPs) for modifiers in HDLM & HDLMM