# \donttest{
## model as in Xinyi & Mokhtarian (2024)
f <-
## ordinal and continuous outcome
twing_status | vmd_ln ~
## selection model
edu_2 + edu_3 + hhincome_2 + hhincome_3 +
flex_work + work_fulltime + twing_feasibility +
att_proactivemode + att_procarowning +
att_wif + att_proteamwork +
att_tw_effective_teamwork + att_tw_enthusiasm + att_tw_location_flex |
## outcome model NTW
female + age_mean + age_mean_sq +
race_black + race_other +
vehicle + suburban + smalltown + rural +
work_fulltime +
att_prolargehouse + att_procarowning +
region_waa |
## outcome model NUTW
edu_2 + edu_3 + suburban + smalltown + rural +
work_fulltime +
att_prolargehouse + att_proactivemode + att_procarowning |
## outcome model UTW
female + hhincome_2 + hhincome_3 +
child + suburban + smalltown + rural +
att_procarowning +
region_waa
fit <- opsr(f, telework_data)
texreg::screenreg(fit, beside = TRUE, include.pseudoR2 = TRUE, include.R2 = TRUE)
# }
Run the code above in your browser using DataLab