我在恢复具有系数限制的线性模型的方差-协方差矩阵时遇到麻烦。不受限制的模型是
y = beta0 + beta1 * x1 + beta2 * x2 + beta3 * x3
受限模型具有等式约束beta1 = beta2。按照通常的惯例，我使用隔离函数 I（） 运行线性模型，如下所示：
受限模型 <- lm（y〜I（x1 + x2）+ x3，data = Data）
但是 vcov（restricted.model）给出了3 * 3矩阵，其列和行名称由[（intercept），I（x1 + x2），x3]给出，而不是4 * 4矩阵的列和行名称为[（intercept），x1，x2，x3]。有什么简单的方法可以在受限模型下恢复完整的方差-协方差矩阵？我的模型比这里的示例大得多，并且我想避免在填写缺失元素时出现任何错误。