如何从线性混合效应模型中获取提取的固定效应估计值的CI？

Question

我有一个线性混合效果模型，如下所示：

model.1 <- lmer(x ~ 0 + treatment + (1|block), data)

我从模型中提取了固定效应估计值：

data$FittedValues <- fixef(model.1)

我制作了拟合值的分布直方图，我需要知道拟合值的95％CI。我尝试了confint（），它给出了每个治疗的CI，但是我需要一个CI用于整套拟合值。我可以在拟合值上运行t.test，但我不认为这给了我正确的答案。

t.test(FittedValues, data = data, 
   alternative = 'two.sided',    
   conf.level = 0.95, 
   na.rm = TRUE)

我是统计数据和R的新手，但我搜索了很长时间，但无法找到答案。如果对于这个委员会来说这个问题太简单，请原谅。