Needs["MultivariateStatistics`"]
n=1000
pdpd={1/6,1/6,1/6,1/6,1/6,1/6}
pdid={1/7,1/5,1/6,1/6.1/6,33/210}
u=RandomInteger[MultinomialDistribution[n,pdpd],1]
v=RandomInteger[MultinomialDistribution[n,pdid],1]
All of the above works.
But u and v appear more than once in subsequent calculations, and
RandomInteger reexecutes at each appearance. The reexecutions change
the values of u and v, and that is not what I want to happen.
Any advice will be appreciated.
John