Normal ordering of powers of the bosonic number operator can be used t
o define a discrete probability distribution associated with the numbe
r of elements of a random set. We represent these random sets by vecto
rs in a Hilbert space and obtain q-bosons.