Dirichlet–Neumann operators arise in many applications in the sciences, and this has inspired a number of studies on their analytical properties. In this paper we further investigate the analyticity properties of Dirichlet–Neumann operators as functions of the boundary shape. In particular, we study the size of the disc of convergence of their Taylor-series representation. For this we use a complexification technique which requires a novel reformulation of the problem, coupled with methods for systems of elliptic partial differential equations. Numerical results to illustrate our theoretical conclusions are presented.