We formulate p-adaptive finite element methods based on the technique called gradient recovery. The answers to two open questions on the design of transition elements and the existence of of approximate formulae for the error will be addressed.