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.
integers can be tested with `assert_equal` (no floating point errors possible)