r - 用 visreg() 说明反向转换的 GLS

Question

我有一个带有对数转换因变量的 GLS。使用 visreg() 进行可视化会给出错误代码：

 library(visreg); library(nlme)
 fit1 <- gls(log(Ozone) ~ Wind , data=na.omit(airquality))
 visreg(fit1,"Wind",trans=exp,ylab="Ozone")

有趣的是，同样适用于 LM：

library(visreg); library(nlme)
fit2 <- lm(log(Ozone) ~ Wind , data=na.omit(airquality))
visreg(fit2,"Wind",trans=exp,ylab="Ozone")