options(digits=4)
data(Ornstein)
mod<-lm(interlocks~assets+nation, data=Ornstein)
Var(mod)
## (Intercept) assets nationOTH nationUK nationUS
## (Intercept) 1.079e+00 -1.588e-05 -1.037e+00 -1.057e+00 -1.032e+00
## assets -1.588e-05 1.642e-09 1.155e-05 1.362e-05 1.109e-05
## nationOTH -1.037e+00 1.155e-05 7.019e+00 1.021e+00 1.003e+00
## nationUK -1.057e+00 1.362e-05 1.021e+00 7.405e+00 1.017e+00
## nationUS -1.032e+00 1.109e-05 1.003e+00 1.017e+00 2.128e+00
hccm(mod)
## (Intercept) assets nationOTH nationUK nationUS
## (Intercept) 1.4808521 -1.640e-05 -1.445e+00 -1.460e+00 -1.455e+00
## assets -0.0000164 2.425e-09 1.079e-05 1.324e-05 1.161e-05
## nationOTH -1.4451983 1.079e-05 6.743e+00 1.432e+00 1.430e+00
## nationUK -1.4600207 1.324e-05 1.432e+00 3.722e+00 1.440e+00
## nationUS -1.4550976 1.161e-05 1.430e+00 1.440e+00 1.840e+00
Run the code above in your browser using DataLab