vcov.NilsEstimate function

Covariance matrix for NILS category estimates