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