Extension of weak orthogonality ͑WO͒ techniques developed by Szalewicz et al. and coupled pair equations derived by Jeziorski et al. to positronium compounds are presented. These methods enable us to calculate many-electronic positronium compounds with explicitly correlated functions. These methods were applied to positronium hydride ͑PsH͒ with Hylleraas-type functions ͑HTF's͒, and the total energies, the positron-electron two-photon annihilation rates, and the positronium binding energies were calculated. Extended coupled pair equations were solved by CEPA͑2͒-type, complete coupled pair ͑CCP͒ procedures, and independent pair approximation ͑IPA͒ of CCP. For the total energies, IPA, CEPA͑2͒, and CCP procedures gave Ϫ0.78899, Ϫ0.78238, and Ϫ0.77471 au, respectively. For the two-photon annihilation rates, the respective procedure gave 2.088, 2.064, and 1.972 ns Ϫ1 , respectively.