R 如何提取特定值的方差-协方差矩阵？

我想提取变量b和c的方差协方差矩阵，并努力找到正确的命令。我的原始数据帧有100多个变量。因此，知道一个开脱的命令将是伟大的

给定数据：

                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

---

athat很快。非常感谢。你不知道下标越界在这种情况下意味着什么？我知道。所有这些都只是矩阵子集。您可以使用
vcov（test）[c（“x”、“z”）、c（“x”、“z”）]
重现示例中的错误
x
和
z
不是设计矩阵的一部分。解决了它！谢谢你的时间，时间很快。非常感谢。你不知道下标越界在这种情况下意味着什么？我知道。所有这些都只是矩阵子集。您可以使用
vcov（test）[c（“x”、“z”）、c（“x”、“z”）]
重现示例中的错误
x
和
z
不是设计矩阵的一部分。解决了它！谢谢你的时间，时间很快。非常感谢。你不知道下标越界在这种情况下意味着什么？我知道。所有这些都只是矩阵子集。您可以使用
vcov（test）[c（“x”、“z”）、c（“x”、“z”）]
重现示例中的错误
x
和
z
不是设计矩阵的一部分。解决了它！谢谢你的时间，时间很快。非常感谢。你不知道下标越界在这种情况下意味着什么？我知道。所有这些都只是矩阵子集。您可以使用
vcov（test）[c（“x”、“z”）、c（“x”、“z”）]
重现示例中的错误
x
和
z
不是设计矩阵的一部分。解决了它！谢谢你抽出时间
               vc.glmCube <- vcov(test[2:3]) # var-cov matrix

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