Added the bit in tstat to do the p-value calculation in log space.
Matrix centering in kinship calculation