cigmales <- dplyr::mutate(cigmales,
age2 = age ^ 2, educ2 = educ ^ 2, educage = educ * age,
age3 = age ^ 3, educ3 = educ ^ 3)
expreg(cigarettes ~ habit + price + restaurant + income + age + age2 + educ + educ2 +
famsize + race | . - habit + reslgth + lagprice + age3 + educ3 + educage,
data = cigmales)
expreg(birthwt ~ cigarettes + parity + race + sex | parity + race + sex +
edmother + edfather + faminc + cigtax, data = birthwt)
Run the code above in your browser using DataLab