I have about 8 features as my predictors in a logistic model I am trying to fit in python. One of the features is TotalAward (Financial Aid) and the second is NEED. I am attempting to predict the likelihood of student retention (not leaving/transferring). Since the feature (NEED) affects the relationship between TotalAward and Retention (y=1), I would like to interact it without interacting any other features in my (
X=df[[ 'x1', 'x2', etc]]). In other statistical software, like stata, you would simply do $X1*X2$. How can I do this in python?