We present an ab initio study of dopant-dopant interactions in beryllium-doped InGaAs. We consider defect formation energies of various interstitial and substitutional defects and their combinations. We find that all substitutional-substitutional interactions can be neglected. On the other hand, interactions involving an interstitial defect are significant. Specially, interstitial Be is stabilized by about 0.9/1.0 eV in the presence of one/two BeGa substitutionals. Ga interstitial is also substantially stabilized by Be interstitials. Two Be interstitials can form a metastable Be-Be-Ga complex with a dissociation energy of 0.26 eV/Be. Therefore, interstitial defects and defectdefect interactions should be considered in accurate models of Be doped InGaAs. We suggest that In and Ga should be treated as separate atoms and not lumped into a single effective group III element, as has been done before. We identified dopant-centred states which indicate the presence of other charge states at finite temperatures, specifically, the presence of Beint +1 (as opposed to Beint +2 at 0K).