Bose-Einstein correlations in<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mrow><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo>+</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mrow><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">−</mml:mi></mml:mrow></mml:msup></mml:mrow></mml:math>collisions

Four data sets collected with the Mark II detector during its operation at the SPEAR and PEP ${e}^{+}$${e}^{\mathrm{\ensuremath{-}}}$ storage rings at SLAC are used to study the Bose-Einstein correlation between pairs and triplets of like-sign charged pions. The data sets represent four different regions of energy available for hadron production: the J/\ensuremath{\psi} at \ensuremath{\surd}s =3.095 GeV, typical hadronic energy of \ensuremath{\approxeq}5 GeV in two-photon events, and annihilation in the energy regions \ensuremath{\surd}s =4--7 GeV above the J/\ensuremath{\psi}, and \ensuremath{\surd}s =29 GeV. The Bose-Einstein correlation is studied as a function of ${Q}^{2}$, the four-momentum difference squared of the pair, and of ${Q}_{3}^{2}$, an analogous quantity defined for triplets. After corrections for Coulomb effects and pion misidentification, pair analyses indicate a nearly full Bose-Einstein enhancement (\ensuremath{\lambda}\ensuremath{\approxeq}1) in the J/\ensuremath{\psi} and the two-photon data, and about half the maximum value in the two higher-energy data sets. The pair analysis parameter r lies within a band of \ifmmode\pm\else\textpm\fi{}0.1 fm around 0.7 fm and is essentially the same for all four data sets. Pion triplet analyses give consistent results for the triplet parameters ${\ensuremath{\lambda}}_{3}$ and ${r}_{3}$. In an attempt to investigate the shape of the pion source, we also study the Bose-Einstein correlation in pion pairs using two-dimensional distributions in components of ${Q}^{2}$.