A function to compute an individual's approximate expected survival time using numerical integration. Note this appears to be unstable; the function is based on R's integrate function. Not intended for general use (yet!).
inputs for the function including the model matrix, frailties, fixed effects and the parameters of the baseline hazard derived from this model
the expected survival time for the individual, obtained by numerical integration of the density function.