Takes the model at convergence and calculates the variance matrix corrected for smoothing parameter uncertainty
cor.var(model)
model
survPen object, see survPen.fit for details
survPen.fit
survPen object with corrected variance Vc