Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
This paper addresses the computational bottleneck of large-scale Cox regression by developing a novel stochastic minimization technique for LogSumExp objectives. It introduces a softplus surrogate for normalizer estimation, enabling unbiased single-sample gradients and achieving improved convergence rates ($O(T^{-1/2})$ averaged objective bound) for smooth convex objectives. For Cox regression specifically, by exploiting nested risk sets through grouped failures and shared auxiliary variables, the method achieves a mean-square rate of $T^{-4/5}$ relative to the full solution, matching its asymptotic distribution.
A single auxiliary variable per group of failures unlocks near-optimal $T^{-4/5}$ convergence for large-scale Cox regression, outperforming stochastic baselines.
Motivated by the computational challenges of large-scale Cox regression, we study stochastic minimization of LogSumExp objectives over large sets. Mini-batch normalizer estimates generally yield biased gradients. We instead use a softplus surrogate that introduces one auxiliary scalar per normalizer and admits unbiased single-sample gradients. For smooth convex LogSumExp objectives, we prove an $O(T^{-1/2})$ averaged objective bound, improving the previous $T^{-1/4}$ analysis. With a strongly convex regularizer on the original variable, we also obtain a last-iterate squared-error rate of $\widetilde{O}(T^{-1})$ without strong convexity in the auxiliary variables. For Cox regression, the normalizers are defined over nested risk sets. We exploit this structure by grouping neighboring failures and sharing one auxiliary variable per group. The resulting compressed objective admits uniform score and curvature bounds that control the errors from grouping and softplus approximation. Together with the general optimization result, these bounds give a mean-square rate of $T^{-4/5}$, up to logarithmic factors, relative to the full Cox solution. The compressed estimator also matches the full estimator's asymptotic distribution. Experiments on synthetic and real survival datasets with slowly decreasing risk sets show a favorable performance relative to stochastic baselines.