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