置信区间或gls {nlme}预测值的SE

Question

我正在运行多变量gls模型：

m <- gls(y ~ x + factor1 + factor2, cor = corPagel(1,phylogeny), weight= ~1/log(n))

我想绘制结果，我可以得到这样的预测值：

newdata <- data.frame(expand.grid( x = mean(x), factor1= unique(factor1), factor2 = unique(factor2)))

predvals <- predict(m,newdata)

虽然我想绘制的两个变量是因子，但我想获得这些预测值的置信区间或标准误差。

有谁知道怎么做？

提前感谢您的帮助！ 干杯， 茉莉