  Multiple harmonic sums are iterated generalizations of harmonic sums.
Recently Dilcher has considered congruences involving q-analogs of these sums
in depth one. In this paper we shall study the homogeneous case for arbitrary
depth by using generating functions and shuffle relations of the q-analog of
multiple harmonic sums. At the end, we also consider some non-homogeneous
cases.
