我有一个具有对数变换因变量的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")