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