Thanks for your feedback, Brenda.

df_panel_with_dummies contains two extra columns: COUNTRY and YEAR which are not part of the LSDV model.

The following should work correctly:

X_col_names = list(df_dummies.columns) + X_var_names

lsdv_model_results.predict(df_panel_with_dummies[X_col_names])

'best