Randomized algorithms and PAC bounds for inverse reinforcement learning in continuous spaces
Angeliki Kamoutsi, Peter Schmitt-Förster, Tobias Sutter, Volkan Cevher, John Lygeros
TL;DR
The paper tackles inverse reinforcement learning in continuous-state and continuous-action MDPs by formulating IRL as an infinite-dimensional linear program over occupancy measures, linking costs to discounted value via $V_c^{\pi}(\nu_0)=\langle \mu_{\nu_0}^{\pi},c\rangle$. It proves that the inverse feasibility set $\mathcal{C}(\pi_E)$ comprises costs for which the observed policy is optimal, characterized through duality with the operator $T_\gamma^*$ and a normalization constraint to exclude trivial rewards. The authors then develop a scalable, finite-dimensional approximation using linear function bases and apply the scenario approach to obtain $\varepsilon$-optimal solutions with probabilistic guarantees, plus sample-based bounds when only finite demonstrations and a generative model are available. These results provide formal guarantees for recovering cost functions in continuous IRL and have practical impact for robotics and safety-critical decision-making, especially where forward solves are expensive or intractable.
Abstract
This work studies discrete-time discounted Markov decision processes with continuous state and action spaces and addresses the inverse problem of inferring a cost function from observed optimal behavior. We first consider the case in which we have access to the entire expert policy and characterize the set of solutions to the inverse problem by using occupation measures, linear duality, and complementary slackness conditions. To avoid trivial solutions and ill-posedness, we introduce a natural linear normalization constraint. This results in an infinite-dimensional linear feasibility problem, prompting a thorough analysis of its properties. Next, we use linear function approximators and adopt a randomized approach, namely the scenario approach and related probabilistic feasibility guarantees, to derive epsilon-optimal solutions for the inverse problem. We further discuss the sample complexity for a desired approximation accuracy. Finally, we deal with the more realistic case where we only have access to a finite set of expert demonstrations and a generative model and provide bounds on the error made when working with samples.
