In this paper we obtain central limit theorems for generalized Pólya urn models with L ≥ 2 colors where one out of K different replacements (actions) is applied randomly at each step. Each possible action constitutes a row of the replacement matrix, which can be nonsquare and random. The actions are chosen following a probability distribution given by an arbitrary function of the proportions of the balls of different colors present in the urn. Moreover, under the same hypotheses it is proved that the covariance matrix of the asymptotic distribution is the solution of a Lyapunov equation, and a procedure is given to obtain the covariance matrix in an explicit form. Some applications of these results to random trees and adaptive designs in clinical trials are also presented.
