对merMod对象的更新给出了不同的拟合

Question

我试图用intercept（null模型）重新构建类merMod的完整模型。但是，使用update.merMod进行重新设定会给出与手动拟合空模型不同的答案，例如：

# Generate random data
set.seed(9)

dat <- data.frame(
  x = do.call(c, lapply(1:5, function(x) rnorm(100, x))),
  random = letters[1:5]
)

dat$y = rnbinom(500, mu = exp(dat$x), size = 1)

library(lme4)

# Get full model
full <- glmer.nb(y ~ x + (1 | random), dat)

# Write out intercept-only model by hand
null <- glmer.nb(y ~ 1 + (1 | random), dat)

# Update 
null2 <- update(full, . ~ 1 -. + (1 | random))

VarCorr(null)

VarCorr(null2)

知道为什么我可以使用update来获取相同的vcov矩阵吗？