quick-start.R
In visreg: Visualization of Regression Models

## ---- include=FALSE-----------------------------------------------------------
library(visreg)
knitr::opts_knit$set(aliases=c(h = 'fig.height', w = 'fig.width'))

## -----------------------------------------------------------------------------
fit <- lm(Ozone ~ Solar.R + Wind + Temp, data=airquality)

## ---- h=5, w=5----------------------------------------------------------------
visreg(fit, "Wind")

## -----------------------------------------------------------------------------
airquality$Heat <- cut(airquality$Temp, 3, labels=c("Cool","Mild","Hot"))
fit <- lm(Ozone ~ Solar.R + Wind*Heat, data=airquality)

## ---- h=4, w=9, out.width='100%'----------------------------------------------
visreg(fit, "Wind", by="Heat")

## ---- h=5, w=5----------------------------------------------------------------
visreg(fit, "Wind", by="Heat", overlay=TRUE)

## ---- h=5, w=6----------------------------------------------------------------
fit <- lm(Ozone ~ poly(Wind, 2)*poly(Temp, 2), data=airquality)
visreg2d(fit, "Wind", "Temp")