Within the framework of local analytic sector subtraction, we present the full analytic integration of double-real and real-virtual local infrared counterterms that enter NNLO QCD computations with any number of massless final-state partons. We show that a careful choice of phase-space mappings leads to simple analytic results, including non-singular terms, that can be obtained with conventional integration techniques.