r - 运行线性回归时如何获得系数的协方差矩阵？

Question

我正在使用该lm()函数运行多元线性回归，我想获得估计系数的协方差矩阵。我怎样才能做到这一点？

这里lme 中系数的协方差（或相关）矩阵有一个函数的答案，lme()我想知道是否可以对lm(). 谢谢！