> print(median(Y[2,-burn.in]))