  We obtain shape derivative formulae for the first eigenvalue of the Robin
$p$-Laplace operator. This result is used to study the variation of the first
eigenvalue with respect to perturbations of the domain. In particular, we prove
that for large values of the boundary parameter, the first eigenvalue is
monotonic with respect to domain inclusion for smooth domains.
