r - 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)

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

有谁知道该怎么做？

在此先感谢您的帮助！干杯，茉莉花