Computes the variance decomposition by simulation

m2.mixt.vdec(model, nsim, stayer_share = 1, ydep = "y2")