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 ?