Reduced order methods are powerful tools for the design and analysis of sophisticated systems, reducing computational costs and speeding up the development process. Among these reduced order methods, the Proper Generalized Decomposition is a well-established one, commonly used to deal with multi-dimensional problems that often suffer from the curse of dimensionality. Although the PGD method has been around for some time now, it still lacks mechanisms to assess the quality of the solutions obtained. This paper explores the dual error analysis in the scope of the PGD, using complementary solutions to compute error bounds and drive an adaptivity process. The energy of the error obtained from the dual analysis is used to determine the quality of the PGD approximations. We define a new adaptivity indicator based on the energy of the error and use it to drive parametric h- and p- adaptivity processes. The results are positive, with the indicator accurately capturing the parameter that will lead to lowest errors.