  We investigate the associated production of heavy quarkonia, with
angular-momentum quantum numbers ^{2S+1}L_J = ^1S_0, ^3S_1, ^1P_1, ^3P_J (J =
0, 1, 2), and photons, Z bosons, and W bosons in photon-photon, photon-hadron,
and hadron-hadron collisions within the factorization formalism of
nonrelativistic quantum chromodynamics providing all contributing partonic
cross sections in analytic form. In the case of photoproduction, we also
include the resolved-photon contributions. We present numerical results for the
processes involving J/psi and chi_{cJ} mesons appropriate for the Fermilab
Tevatron, CERN LHC, DESY TESLA, operated in the e^+ e^- and gamma gamma modes,
and DESY THERA.
