1. (important) update \pibeta and \pigamma using slice sampling; and keep tracking them for spike and slab prior
2. corresponding coef function and plot function