Plot prediction interval in r