Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Improve report() for anova() outcome of regression models #53

Open
atnplab opened this issue Oct 3, 2019 · 0 comments
Open

Improve report() for anova() outcome of regression models #53

atnplab opened this issue Oct 3, 2019 · 0 comments
Labels
enhancement 💥 Implemented features can be improved or revised

Comments

@atnplab
Copy link

atnplab commented Oct 3, 2019

I noticed that report() of an anova() object does not provide correct/complete output for regression models (eg lm, lmer, glm).

In particular, it could be very helpful if report() could give an anova-like table of main and interaction effects of lmer or complex multiple lm() as coefficients parameters in these cases is sometimes hard to understand or some coefficients do not make sense empirically (i.e. expected difference for a numerical (age) coefficient of 0 between the levels of a categorical (use/non-use of electronic cigarette) coefficient for the VD (motivation to quit smoking)) and as reference for categorical comparison is taken without control on experimental or theoretical reasons in R.
Use of report for these models could be helpful for providing also effect size from sjstats and other relevant parameters of an anova() object to be reported in papers.

This is secondary but might help to report relevant information: report() does not work for anova() of model comparisons. Honestly, I do not know whether there are guidelines for model comparison reporting in psychology/neuroscience. However, report() could help for picking the right parameters (e.g. AIC, BIC or Chi2 and so on) and correct writing.

Lastly, would it be possible to get main and interaction effects for stan_lm/glm/lmer model with report()? the coefficient parameters output is the same as classic regression models and would be helpful to get the general effect of a predictor instead of its coefficient estimation.

How could we do it?
For instance, significance for fixed effects in lmer could be done with anova(model.lmer, dff = "Kenward-Roger) (Luke, 2017) or for other regression models as anova(model.lm) or anova(model.glm, test ="F").

I am not an expert in R programming and cannot help with actual coding, this is just a suggestion that I noticed during practical use of report() and R in general.

I appreciate your work!

Luke, S.G. Behav Res (2017) 49: 1494. https://doi.org/10.3758/s13428-016-0809-y

@strengejacke strengejacke added bug 🐛 Something isn't working enhancement 💥 Implemented features can be improved or revised labels Mar 19, 2020
@strengejacke strengejacke removed the bug 🐛 Something isn't working label Sep 22, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement 💥 Implemented features can be improved or revised
Projects
None yet
Development

No branches or pull requests

2 participants