It is by now well known that, at subleading power in scale ratios, factorization theorems for high-energy cross sections and decay amplitudes contain endpoint-divergent convolution integrals. The presence of endpoint divergences hints at a violation of simple scale separation. At the technical level, they indicate an unexpected failure of dimensional regularization and the MS subtraction scheme. In this paper we start a detailed discussion of factorization at subleading power within the framework of soft-collinear effective theory. As a concrete example, we factorize the decay amplitude for the radiative Higgs-boson decay h → γγ mediated by a b-quark loop, for which endpoint-divergent convolution integrals require both dimensional and rapidity regulators. We derive a factorization theorem for the decay amplitude in terms of bare Wilson coefficients and operator matrix elements. We show that endpoint divergences caused by rapidity divergences cancel to all orders of perturbation theory, while endpoint divergences that are regularized dimensionally can be removed by rearranging the terms in the factorization theorem. We use our result to resum the leading double-logarithmic corrections of order α n s ln 2n+2 (−M 2 h /m 2 b ) to the decay amplitude to all orders of perturbation theory.