  We compute the skewness of the matter distribution arising from non-linear
evolution and from non-Gaussian initial perturbations. We apply our result to a
very generic class of models with non-Gaussian initial conditions and we
estimate analytically the ratio between the skewness due to non-linear
clustering and the part due to the intrinsic non-Gaussianity of the models. We
finally extend our estimates to higher moments.
