Takes the model at convergence and calculates the robust variance matrix accounting for correlated survival times
robust.var(model, data, cluster.name, n.legendre = 50)survPen object with robust variance Vr
survPen object, see survPen.fit for details
original dataset
name of cluster variable in data
number of nodes for Gauss-Legendre quadrature; default is 50