ENH: PHReg formula improvements 

This PR allows extra arguments of array type to be passed to the PHReg from_formula method as variable names.  The data are read out of the same data object that is used for the formula.
I don't feel particularly strongly about it but we have `missing='none'` as the default everywhere. The idea is not to do any unnecessary computations unless asked for.