We present an analysis of convergence of a quasi-regression Monte Carlo method proposed by Glasserman and Yu (2004). We show that the method surely converges to the true price of an American option even under multiple underlyings via polynomial chaos expansion and weaker conditions than those used in Glasserman and Yu (2004). Further, we show the number of simulation paths grows exponentially in the number of basis functions to obtain convergence in implementing the method. Finally, we propose a rate of convergence considering regularity of value functions.