  We perform a series of numerical experiments to study how the nonlinear
metallicity--color relations predicted by different stellar population models
affect the color distributions observed in extragalactic globular cluster
systems. % We present simulations in the $UBVRIJHK$ bandpasses based on five
different sets of simple stellar population (SSP) models. The presence of
photometric scatter in the colors is included as well. % We find that unimodal
metallicity distributions frequently ``project'' into bimodal color
distributions. The likelihood of this effect depends on both the mean and
dispersion of the metallicity distribution, as well as of course on the SSP
model used for the transformation. % Adopting the Teramo-SPoT SSP models for
reference, we find that optical--to--near-IR colors should be favored with
respect to other colors to avoid the bias effect in globular cluster color
distributions discussed by \citet{yoon06}. In particular, colors such as \vh\
or \vk are more robust against nonlinearity of the metallicity--color relation,
and an observed bimodal distribution in such colors is more likely to indicate
a true underlying bimodal metallicity distribution. Similar conclusions come
from the simulations based on different SSP models, although we also identify
exceptions to this result.
