This paper is concerned with identification of linear parameter varying (LPV) systems in an input-output setting with Box-Jenkins (BJ) model structure. Classical linear time invariant prediction error method (PEM) is extended to the LPV PEM. Under the new LPV framework, identification of two types of input-output LPV models is considered: one is based on parameter interpolation and the other is based on model interpolation. The effectiveness of the proposed solution is validated by comparison with other existing LPV identification approaches through simulation examples and demonstrated by experiment studies. (C) 2011 Elsevier Ltd. All rights reserved.