We compute the decay rates for the exclusive decays B+/- --> (eta',eta)(K+/-,K*(+/1)) and B-0 --> (eta'eta)(K-0,K*(0)) in a QCD-improved factorization framework by including the contribution from the process b --> sgg --> s(eta',eta) through the QCD anomaly. This method provides an alternative estimate of the contribution b --> sc (c) over bar --> s(eta',eta) to these decays as compared to the one using the intrinsic charm content of the eta' and eta mesons determined through the decays J/psi --> (eta,eta',eta(c))gamma . The advantage of computing the relevant matrix elements via the QCD anomaly governing the transition gg --> (eta',eta) is that there is no sign ambiguity in these contributions relative to the matrix elements from the rest of the operators in the weak effective Hamiltonian. Numerically, the QCD anomaly method and the one using the radiative decays J/psi --> (eta,eta',eta(c))gamma give similar branching ratios for the decays of interest here. The resulting branching ratios are compared with the CLEO data on B+/- --> eta'K+/- and B-0 --> eta'K-0 and predictions are made for the rest. (C) 1998 Elsevier Science B.V.