Gaussian process approximations for multicolor Pólya urn models

Motivated by mathematical tissue growth modelling, we consider the problem of approximating the dynamics of multicolor Polya urn processes that start with large numbers of balls of different colors and run for a long time. Using strong approximation theorems for empirical and quantile processes, we establish Gaussian process approximations for the Polya urn processes. The approximating processes are sums of a multivariate Brownian motion process and an independent linear drift with a random Gaussian coefficient. The dominating term between the two depends on the ratio of the number of time steps n to the initial number of balls N in the urn. We also establish an upper bound of the form for the maximum deviation over the class of convex Borel sets of the step-n urn composition distribution from the approximating normal law.