My logistic regression model includes the following interaction:
x1*x2*time*time^2
Where x1 and x2 are both dummy variables, and time is a variable ranging from 0 to 25.
In a previous model that only contains a x1*x2 interaction, I was able to calculate the AME of the interaction term using the following code (I am using mi imputed data):
Is it also possible for me to calculate individual AME values to be presented in a regression table for the interaction terms that arise from my more complex model, namely:
x1*x2
x1*time
x2*time
x1*time^2
x2*time^2
x1*time*time^2
x2*time*time^2
If so, how should I go about doing it? (I'm using StataSE v16.01)
At present, I am presenting log-odds in the regression table to give an initial idea of the relationships, and then plotting the AMEs over time, but I have received feedback asking me to already present AMEs in the regression table as well, which is why I am asking.
Thank you!
x1*x2*time*time^2
Where x1 and x2 are both dummy variables, and time is a variable ranging from 0 to 25.
In a previous model that only contains a x1*x2 interaction, I was able to calculate the AME of the interaction term using the following code (I am using mi imputed data):
Code:
mimrgns x1, dydx(x2) predict(pr) pwcompare
x1*x2
x1*time
x2*time
x1*time^2
x2*time^2
x1*time*time^2
x2*time*time^2
If so, how should I go about doing it? (I'm using StataSE v16.01)
At present, I am presenting log-odds in the regression table to give an initial idea of the relationships, and then plotting the AMEs over time, but I have received feedback asking me to already present AMEs in the regression table as well, which is why I am asking.
Thank you!

Comment