I'm working with panel data using the plm package in R to estimate a fixed effects model that includes interaction terms. I'm currently able to get my model results using plm and display them with stargazer, but I'm looking for a more refined way to present these results, specifically focusing on the interaction terms. Additionally, I need to calculate and display the marginal effects of these interaction terms, including their p-values, in the same table. I'm also interested in plotting these marginal effects to better understand their impact.
Here's an example of my plm model setup:
fe_model_interaction <- plm(life_satisfaction ~
data = data_analyse_mother,
index = c("pid", "syear"),
model = "within")
My questions are:
What is the best way to display the plm model results, especially for interaction terms, in a more detailed and customizable format than what stargazer offers? I'm looking for a method to include standard errors, coefficients, and p-values in a clean and professional-looking table.
How can I calculate and include the marginal effects of these interaction terms in the results table, along with their p-values?
What are the recommended approaches or packages in R for plotting the marginal effects of interaction terms in a fixed effects model estimated with plm?