如何使用predict获得标准偏差？

我正试图策划一个三方互动。首先，我使用

expand.grid

创建变量的不同组合。然后我使用

predict

获得预测值。它工作得很好，但是，我想得到标准偏差或置信区间。一旦我将选项

se.fit=T

包含在

predict

中，我就会得到一些我无法解释的结果，因为两个组合得到一个预测值，另外两个得到四个结果的范围。我不确定在这种情况下，

se.fit=T

选项是否错误，或者问题是否在我对结果的解释中。有人能帮忙吗

这是不带

se.fit

选项的代码：

set.seed(20170925) 
dat <- data.frame(time=gl(n = 2,k = 5000),
                  y= rnorm(10000, mean=1000, sd=400),
                  a=factor(sample(1:2, 10000, replace=TRUE)),
                  b=factor(sample(1:2, 10000, replace=TRUE)),
                  c=factor(rep(rep(1:3,times=c(2500,2500,5000)),2)))

m = lm(y ~ a*b*c, data=dat)
pred <- expand.grid(a=factor(1:2),b=factor(1:2), c=factor(1, levels=c(1,2,3)))
pred$y <- predict(m,pred)

---

predict

现在为您提供一个列表，您只需选择其中的元素，可能如下所示：

p
pred$y <- predict(m,pred, se.fit = T)