fit@output$confidenceIntervals doesn't give estimate, just lbound and ubound