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