new old script 3vex_estim.R