profile_phat_k = Profile(phats, i=k)