hence the AGN effect can be ignored in this work.

@mos12 also reported that the galaxy color depends on the SFR uncertainties. The adopted fitting formula in this work is found to have a scatter of residual $\sim$ 0.19 for the star-forming galaxies and $\sim$ 0.45 for the quiescent galaxies. That is, the method may not yield precise SFR measurements for star-forming and quiescent galaxies. To estimate the uncertainty in our analysis that is introduced by this method, we have performed Monte Carlo simulations. Given a median SSFR for the quiescent and star-forming galaxies and a quiescent fraction, we construct a distribution of SSFR with double-Gaussian peaks to mimic distinct star-forming and quiescent populations. We then randomly perturb the SSFRs with a scatter 0.45 and 0.19 for the quiescent and star-forming population, respectively, and repeat the same procedure for ten thousand times to obtain the bias and deviation for the median SSFR of quiescent and star-forming galaxies and the quiescent fraction relative to our original inputs. Given different inputs based on our results in four different mass ranges, including the median SSFR for the quiescent and star-forming galaxies and the corresponding quiescent fraction, we find that the bias and scatter 