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