How to find the upper and lower CI values for new observation prediction using an ols model in Python? I want to plot those CI in a same plot

Question

Like in R we can set the Confidence interval value as we want with a single line:

model = lm(formula = assay~time + batch,data = AssayData) anova(model) summary(model) ne