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