We have used a path integral Monte Carlo technique to simulate positronium (Ps) in a cavity. The primitive propagator is used, with a pair of interacting chains representing the positron and electron. We calculate the energy and radial distribution function for Ps enclosed in a hard, spherical cavity, and the polarizability of the model Ps in the presence of an electrostatic field. We find that the positron distribution near the hard wall differs significantly from that for a single particle in a hard cavity. This leads to systematic deviations from predictions of free-volume models which treat Ps as an effective, single particle. A virial-type estimator is used to calculate the kinetic energy of the particle in the presence of hard walls. This estimator is found to be superior to a kinetic-type estimator given the interaction potentials, cavity sizes, and chain lengths considered in the current study.