Return to Question

improvement over the question

Source Link

edit approved Sep 16, 2018 at 10:50

This is the results of my anova(glm())anova(glm()) and the post-hoc analyses emmeans()emmeans() :

         Df Deviance Resid. Df Resid. Dev  Pr(>Chi)    
NULL                                 515     1336.6              
Type_product        3   32.544       512     1304.0 4.019e-07 ***
Exhaustion_product  9   92.167       503     1211.8 5.977e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1[/code]

gl=glm(Effort ~  Type_product +  Exhaustion_product, family=poisson , data=vect)
library("emmeans")
emmp <- emmeans( gl, pairwise ~ Type_product)
summary( emmp, infer=TRUE)

 contrast    estimate         SE  df    asymp.LCL   asymp.UCL z.ratio p.value
 c - f    -0.14084934 0.04221684 Inf -0.249305743 -0.03239295  -3.336  0.0047
 c - m     0.33882907 0.08050197 Inf  0.132016967  0.54564118   4.209  0.0002
 c - s     0.31167356 0.12274400 Inf -0.003659682  0.62700680   2.539  0.0541
 f - m     0.47967842 0.08339482 Inf  0.265434484  0.69392235   5.752  <.0001
 f - s     0.45252290 0.12076017 Inf  0.142286189  0.76275962   3.747  0.0010
 m - s    -0.02715551 0.14317625 Inf -0.394979861  0.34066883  -0.190  0.9976

There is significant effect of Effort and Type_product (χ2(3)=32.5, p<0.001). Post-hoc test report decreasing Effort of C comparing to F.

Firstly, I want to be sure that my GLM reporting is correct.

After that, I can't found any APA format reporting for post-hoc test (p-value, estimate?)

How to report these values ?

This is the results of my anova(glm()) and the post-hoc analyses emmeans()

         Df Deviance Resid. Df Resid. Dev  Pr(>Chi)    
NULL                                 515     1336.6              
Type_product        3   32.544       512     1304.0 4.019e-07 ***
Exhaustion_product  9   92.167       503     1211.8 5.977e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1[/code]

gl=glm(Effort ~  Type_product +  Exhaustion_product, family=poisson , data=vect)
library("emmeans")
emmp <- emmeans( gl, pairwise ~ Type_product)
summary( emmp, infer=TRUE)

 contrast    estimate         SE  df    asymp.LCL   asymp.UCL z.ratio p.value
 c - f    -0.14084934 0.04221684 Inf -0.249305743 -0.03239295  -3.336  0.0047
 c - m     0.33882907 0.08050197 Inf  0.132016967  0.54564118   4.209  0.0002
 c - s     0.31167356 0.12274400 Inf -0.003659682  0.62700680   2.539  0.0541
 f - m     0.47967842 0.08339482 Inf  0.265434484  0.69392235   5.752  <.0001
 f - s     0.45252290 0.12076017 Inf  0.142286189  0.76275962   3.747  0.0010
 m - s    -0.02715551 0.14317625 Inf -0.394979861  0.34066883  -0.190  0.9976

There is significant effect of Effort and Type_product (χ2(3)=32.5, p<0.001). Post-hoc test report decreasing Effort of C comparing to F.

Firstly, I want to be sure that my GLM reporting is correct. After that, I can't found any APA format reporting for post-hoc test (p-value, estimate?)

How to report these values ?

This is the results of my anova(glm()) and the post-hoc analyses emmeans() :

         Df Deviance Resid. Df Resid. Dev  Pr(>Chi)    
NULL                                 515     1336.6              
Type_product        3   32.544       512     1304.0 4.019e-07 ***
Exhaustion_product  9   92.167       503     1211.8 5.977e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1[/code]

gl=glm(Effort ~  Type_product +  Exhaustion_product, family=poisson , data=vect)
library("emmeans")
emmp <- emmeans( gl, pairwise ~ Type_product)
summary( emmp, infer=TRUE)

 contrast    estimate         SE  df    asymp.LCL   asymp.UCL z.ratio p.value
 c - f    -0.14084934 0.04221684 Inf -0.249305743 -0.03239295  -3.336  0.0047
 c - m     0.33882907 0.08050197 Inf  0.132016967  0.54564118   4.209  0.0002
 c - s     0.31167356 0.12274400 Inf -0.003659682  0.62700680   2.539  0.0541
 f - m     0.47967842 0.08339482 Inf  0.265434484  0.69392235   5.752  <.0001
 f - s     0.45252290 0.12076017 Inf  0.142286189  0.76275962   3.747  0.0010
 m - s    -0.02715551 0.14317625 Inf -0.394979861  0.34066883  -0.190  0.9976

There is significant effect of Effort and Type_product (χ2(3)=32.5, p<0.001). Post-hoc test report decreasing Effort of C comparing to F.

Firstly, I want to be sure that my GLM reporting is correct.

After that, I can't found any APA format reporting for post-hoc test (p-value, estimate?)

How to report these values ?

Source Link

asked Sep 16, 2018 at 10:47

ranell

Report GLM and Posthoc with emmeans in APA format

This is the results of my anova(glm()) and the post-hoc analyses emmeans()

         Df Deviance Resid. Df Resid. Dev  Pr(>Chi)    
NULL                                 515     1336.6              
Type_product        3   32.544       512     1304.0 4.019e-07 ***
Exhaustion_product  9   92.167       503     1211.8 5.977e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1[/code]

gl=glm(Effort ~  Type_product +  Exhaustion_product, family=poisson , data=vect)
library("emmeans")
emmp <- emmeans( gl, pairwise ~ Type_product)
summary( emmp, infer=TRUE)

 contrast    estimate         SE  df    asymp.LCL   asymp.UCL z.ratio p.value
 c - f    -0.14084934 0.04221684 Inf -0.249305743 -0.03239295  -3.336  0.0047
 c - m     0.33882907 0.08050197 Inf  0.132016967  0.54564118   4.209  0.0002
 c - s     0.31167356 0.12274400 Inf -0.003659682  0.62700680   2.539  0.0541
 f - m     0.47967842 0.08339482 Inf  0.265434484  0.69392235   5.752  <.0001
 f - s     0.45252290 0.12076017 Inf  0.142286189  0.76275962   3.747  0.0010
 m - s    -0.02715551 0.14317625 Inf -0.394979861  0.34066883  -0.190  0.9976

There is significant effect of Effort and Type_product (χ2(3)=32.5, p<0.001). Post-hoc test report decreasing Effort of C comparing to F.

Firstly, I want to be sure that my GLM reporting is correct. After that, I can't found any APA format reporting for post-hoc test (p-value, estimate?)

How to report these values ?