The incorporation of interstitial hydrogen in yttria was studied by means of ab initio calculations based on density-functional theory (DFT) and muonium spin polarization spectroscopy (μSR). The density-functional calculations, based on a semilocal functional within the GGA-PBE and a hybrid functional, uncovered multiple geometrical configurations for the neutral, H0, and the negatively charged, H , states of hydrogen, thus demonstrating the existence of metastable minimum-energy sites. It was observed that the low-energy configurations for H0 and H are similar: they prefer to relax in deep, interstitial sites, whereas the equilibrium configurations for the positively charged state, H , were bond-type configurations with the hydrogen forming a covalent O-H bond with an O anion. For all neutral and negative configurations, localized defect levels were found inside the gap. Overall, the results for the formation energies obtained by the two different functionals are qualitatively similar; an amphoteric behavior was found for hydrogen after considering the lowest-energy structures for each charge state. The calculated acceptor transition level, obtained by the hybrid functional and seen near midgap, is consistent with μSR data from literature. The results are consistent with the present μSR data, where the observed diamagnetic signal is attributed to a donor-like muonium at the oxygen-bonded configurations and the paramagnetic signal to an acceptor-like deep muonium at the interstitial sites.