Update estimate.R