inst/tests/visreg2d-glm.R

In visreg: Visualization of Regression Models

data("birthwt",package="MASS")
birthwt$race <- factor(birthwt$race,labels=c("White","Black","Other"))
birthwt$smoke <- factor(birthwt$smoke,labels=c("Nonsmoker","Smoker"))

## Basic
fit <- glm(low~age+race+smoke+lwt,data=birthwt,family="binomial")
visreg2d(fit,"age","lwt")
visreg2d(fit,"age","lwt",scale="response")

## Transformation of X
fit <- glm(low~age+race+smoke+log(lwt),data=birthwt,family="binomial")
visreg2d(fit,"age","lwt")
visreg2d(fit,"age","lwt",type="contrast")
visreg2d(fit,"age","lwt",scale="response")

## Cond
visreg2d(fit,"age","lwt",scale="response",cond=list(smoke='Smoker'))

## Factors
visreg2d(fit,"age","race",cond=list(smoke='Smoker'), ylab="")
visreg2d(fit,"race","age",cond=list(smoke='Smoker'), xlab="")