Create estimation.R