updated PP code for fixed sigma.squared.beta