# NOT RUN {
library("Formula")
inst <- ~ sic3 + k_serv + inv + engsci + whitecol + skill + semskill + cropland +
pasture + forest + coal + petro + minerals + scrconc + bcrconc +
scrcomp + bcrcomp + meps +
kstock + puni + geog2 + tenure + klratio + bunion
tradeprotection <- dplyr::mutate(tradeprotection,
y = ntb / (1 + ntb),
x1 = exports / imports / elast,
x2 = cap * x1)
smbltest(Formula::as.Formula(y ~ x1 + x2 + labvar, inst), tradeprotection)
# }
Run the code above in your browser using DataLab