Renamed cp and pc and now using glm formulation

cp is now cpa and pc is now pcv. Last name change, hopefully, for good!
These functions are now using glm formulation, i.e. y ~ x1 + x2 + x3 …