I have trained a linear regression model based on some data. Now I have new data, and I need to find a way to calculate the CIs for each $Y_{new}$ sample. Meaning now I have n samples
in the features $X_{new}$ and the same n samples in the labels $Y_{new}$.
I was able to calculate the CIs for $E[Y_{new}|X_{new}]$ for each point in $Y_{new}$ (so overall n confidence intervals), using percentiles bootstrap. For each iteration of the bootstrap, I resample the data, fit a linear regression model on the resampled data, then predict $Y_{new}$ for $X_{new}$ and collect the predictions in an array. Eventually I found the confidence intervals.
What about $Y_{new}$ using bootstrap? just an explenation of this would be really helpful no need for the implementation in code. Thanks you!