抵消GLM因子R-标准误差

使用内置的

mtcars

数据集，我刚刚运行了以下代码：

my_mtcars <- mtcars
my_mtcars$cyl <- as.factor(my_mtcars$cyl)
my_mtcars$gear <- as.factor(my_mtcars$gear)

我想偏移

cyl

。如果我将偏移设置为与glm参数估计值相同，则应获得与

档位相同的模型

my_glm <- glm(mpg ~ cyl + gear, data = my_mtcars)
summary(my_glm)

my_mtcars$cyl_offset <- case_when(
  my_mtcars$cyl == 4 ~ 0,
  my_mtcars$cyl == 6 ~ -6.656,
  my_mtcars$cyl == 8 ~ -10.542
)

my_mtcars$cyl_offset通过使用偏移量，您正在修复cyl的参数估计值。由于完整模型中的参数估计值之间存在非零协方差，因此，如果修复cyl，您将减小齿轮参数估计值的方差。相反，如果向模型添加共线参数，则会增大方差，假设它与此相反。只是一个猜测，但是由于模型中不包括
cyl
（只包括它以抵消模型中的其他项），那么它就不会包括在协方差矩阵中，协方差矩阵用于计算标准误差。这真的更适合。。。
my_glm <- glm(mpg ~ offset(cyl_offset) + gear, data = my_mtcars)
summary(my_glm)