help seeking for prediction of fixed effect model

[gotodockforevildragon]

now i get a formula by FE MODEL, but i wanna use it to predict the upper and lower bound of 'Y'
however, it seems that the i could not predict if based on the orginial data by commands like 'get_prediction'
because this command is for OLS, not for the panel data or panel model
if i solve it by Stata, i just need to use "predict" and it will ouput the prediction directly
so are there any other commands for my demand QAQ ！！！

This is my code:

filtered_df['Y_winsorized'] = winsorize(filtered_df['Y'], limits=[0.2, 0.2])

filtered_df = filtered_df.set_index(["material","date"])
exog = sm.add_constant(filtered_df[['X1','X2']])
result = PanelOLS(filtered_df['Y_winsorized'], exog, entity_effects=True, time_effects=True)
result = result.fit()
print(result)

pred = result.get_prediction(exog)
conf_int = pred.conf_int()  
pred_int = pred.conf_int(observed=True)  

this is the warning and error:

Error:
AttributeError: 'PanelEffectsResults' object has no attribute 'get_prediction'

Larz60+ write Jun-18-2024, 08:03 AM:
Please post all code, output and errors (it it's entirety) between their respective tags. Refer to BBCode help topic on how to post. Use the "Preview Post" button to make sure the code is presented as you expect before hitting the "Post Reply/Thread" button.
I have added tags for you this time. Please use BBCode tags on future posts.