Update estimation.R