add_p function

Add p-values