library(MASS)
dat <- MASS::birthwt

##   race bwt
## 1    1  96
## 2    2  26
## 3    3  67

race 1=White 2=Black 3=Other