added computation of mean f samples in bedpostx