I have developed a bunch of models to predict house prices. It is a 3 fold process:
I fit a gbm (first_model) and get the first prediction (first_pred),
there are some sub-models (simple lineer regression) for garage (sub_pred1), terrace(sub_pred2), energy(sub_pred3) and state of building(sub_pred4) to predict their impact. And get second_pred = first_pred * (1 + sub_pred1 + sub_pred2 + sub_pred3 + sub_pred4)
One last factor about quality is added to the second_pred. So main_pred = second_pred * (1 + quality) where quality in [0.1, -0.1]
Now I am interested in finding prediction intervals for the main_preds
I have been playing around conforming prediction through nonconformist library for python and I could directly use it if I had only 1 gbm model. But now I wonder what would be the best and doable way to come up with prediction intervals when having those different steps affecting the final prediction?
Thanks in advance for all helps and guidance!