Update estimateTheta.R