In USGS-R/WREG: USGS WREG v. 2.02

knitr::opts_chunk$set(echo = TRUE)

Summary text

print(wregOUT)

Figures

Fitted vs. Residuals

plot(wregOUT$fitted.values,wregOUT$residuals,
                        xlab="Fitted values",ylab="Residuals",
                        main="Fitted vs Residual")

Y vs Leverage

plot(wregOUT$Y,wregOUT$ResLevInf$Leverage,
                          xlab="Y",ylab="Leverage",
                          main="Y vs Leverage")
                     abline(h=wregOUT$LevLim,lty=2,col="red")
                     par(mar=c(0, 0, 0, 0))
                     # c(bottom, left, top, right)
                     plot.new()
                     legend('center',"Critical value", lty = 2,
                            col="red",ncol=1,bty ="n")

Y vs Influence

layout(rbind(1,2), heights=c(7,1))
                     plot(wregOUT$Y,wregOUT$ResLevInf$Influence,
                          xlab="Y",ylab="Influence",
                          main="Y vs Influence")
                     abline(h=wregOUT$InflLim,lty=2,col="red")
                     par(mar=c(0, 0, 0, 0))
                     # c(bottom, left, top, right)
                     plot.new()
                     legend('center',"Critical value", lty = 2,
                            col="red",ncol=1,bty ="n")

USGS-R/WREG documentation built on May 9, 2019, 6:48 p.m.