Results of fitting a double (z < 3) and single Schechter (z ≥ 3) functional forms to the observed galaxy stellar mass function. Measurements are binned in redshift and inferred from both χ² minimization (red) and likelihood methods (evolving colors) taken at the median of each parameter posterior distribution. In both cases, Schechter functions are convolved with parameterized, redshift-dependent kernels in δ M (upper right subpanels) to account for Eddington bias (dashed lines) which is then removed in the corrected fit (solid lines). For completeness, the corrected maximum likelihood solutions are also indicated (dash-dot colored curves). In the lower panels, fitting residuals between the uncorrected kernel convolved fits and the data are shown by both the relative fractional difference (dashed curve) and the difference weighted by the uncertainty (square points). Corresponding parameters are shown in Table C.1.