In this paper, we study Lp-boundedness properties for higher order Littlewood-Paley g-functions in the Bessel setting. We use the Calderón-Zygmund theory in a homogeneous-type space (in the sense of Coifman and Weiss) ((0, ∞), d, γα), where d represents the usual metric on (0, ∞) and γα denotes the doubling measure on (0, ∞) with respect to d defined by dγα(x) = x2α+1dx, with α > −1/2.