  We seek to understand the impact on shape estimators obtained from circular
and elliptical shapelet models under two realistic conditions: (a) only a
limited number of shapelet modes is available for the model, and (b) the
intrinsic galactic shapes are not restricted to shapelet models.
  We create a set of simplistic simulations, in which the galactic shapes
follow a Sersic profile. By varying the Sersic index and applied shear, we
quantify the amount of bias on shear estimates which arises from insufficient
modeling. Additional complications due to PSF convolution, pixelation and pixel
noise are also discussed.
  Steep and highly elliptical galaxy shapes cannot be accurately modeled within
the circular shapelet basis system and are biased towards shallower and less
elongated shapes. This problem can be cured partially by allowing elliptical
basis functions, but for steep profiles elliptical shapelet models still depend
critically on accurate ellipticity priors. As a result, shear estimates are
typically biased low. Independently of the particular form of the estimator,
the bias depends on the true intrinsic galaxy morphology, but also on the size
and shape of the PSF.
  As long as the issues discussed here are not solved, the shapelet method
cannot provide weak-lensing measurements with an accuracy demanded by upcoming
missions and surveys, unless one can provide an accurate and reliable
calibration, specific for the dataset under investigation.
