  A sequence of rational numbers as a generalization of the sequence of
Bernoulli numbers is introduced. Sums of products involving the terms of this
generalized sequence are then obtained using an application of the Fa\`a di
Bruno's formula. These sums of products are analogous to the higher order
Bernoulli numbers and are used to develop the closed form expressions for the
sums of products involving the power sums $\displaystyle
\Psi_k(x,n):=\sum_{d|n}\mu(d)d^k S_k(\frac{x}{d}), n\in\mathbb{Z}^+$ which are
defined via the M\"obius function $\mu$ and the usual power sum $S_k(x)$ of a
real or complex variable $x.$ The power sum $S_k(x)$ is expressible in terms of
the well known Bernoulli polynomials by $\displaystyle
S_k(x):=\frac{B_{k+1}(x+1)-B_{k+1}(0)}{k+1}.$
