# NOT RUN {
data("mwp", package = "feisr")
feis.mod <- feis(lnw ~ marry + enrol | exp,
data = mwp, id = "id")
vcov(feis.mod)
all.equal(vcov(feis.mod), feis.mod$vcov) # FALSE: not equal, because vcov() unscaled
all.equal(vcov(feis.mod, scale = TRUE), feis.mod$vcov) # equal
# }
Run the code above in your browser using DataLab