We present the results of the second-order Fourier coefficients (v2), which characterizes the azimuthal anisotropy, in p+p, p+Au, and d+Au collisions at sqrt(sNN) = 200 GeV. The v2 values are calculated via the two-particle correlation method. To investigate the role of initial geometry effects, the measured v2 is extracted by two different kinematic selections. The two results are compared and its ratio will be shown as a function of pT and charged-particle multiplicity to evaluate the validity of the flow factorization. The results are compared to AMPT calculations, which reproduces the flow factorization breaking in small collisions systems.