如何为特定值提取方差协方差矩阵？

Question

我想提取变量b和c的方差协方差矩阵，并且有一些努力找到正确的命令。我的原始数据框有100多个变量。因此要知道exctract的命令会很棒

鉴于数据：

                a<-rnorm(1000, mean = 0, sd = 1)
                b<-rnorm(1000, mean = 0, sd = 1)
                c<-rnorm(1000, mean = 0, sd = 1)
                d<-rbinom(1000, size = 1, prob = .5)
                e<-rbinom(1000, size = 1, prob = .5)
                f<-rbinom(1000, size = 1, prob = .5)
                data<-data.frame(a,b,c,d,e,f)

               test<-glm(a~b+c+d+e+f,data=data)

               pe.glmCube<-test$coefficients[2:3] # point estimates

我尝试使用方差矩阵。但这样做似乎毫无意义：

               vc.glmCube <- vcov(test[2:3]) # var-cov matrix

Answer 1

vcov(test)[c("b", "c"), c("b", "c")]
#              b             c
#b  1.083964e-03 -2.532682e-05
#c -2.532682e-05  9.779278e-04